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INTHODUDCION 


Iv] 2 de agosto de 2027 habra un eclipse total de Sol que sera vi- 
sible en Espafia, Marruecos, Libia, Argelia, Egipto, Arabia Saudi, 
Yemen y Somalia. El eclipse comenzar4 a las 8 horas, 41 minutos 
y 36 segundos (tiempo universal) y finalizara a las 11 horas, 27 
minutos y 41 segundos; su fase maxima tendra lugar a las 10 ho- 
ras y 6 minutos. 

Tales de Mileto predijo un eclipse el 28 de mayo de 585 a.C. 
No es que tengamos documentos histéricos de tanta precision, 
sino que, al igual que podemos calcular lo que ocurrira en el fu- 
{uro, también podemos calcular lo que ocurrié en el pasado. Los 
astr6nomos han determinado que en esa fecha hubo un eclipse 
total de Sol visible desde Asia Menor, y durante muchos afios no 
hubo otro de similares caracteristicas visible en la zona, por lo 
que debio de ser ese el que Tales predijo. 

La mecanica clasica es determinista. Rige el movimiento de 
planetas, satélites y proyectiles. Sus leyes fueron establecidas 
por Galileo y Newton, y desarrolladas durante los siglos xvm, xvm 
y xx por fisicos y matematicos. Conocidas las fuerzas que actian 
sobre un cuerpo, y dadas su posicion y su velocidad iniciales, las 
leyes de la mecanica clasica permiten calcular con toda preci- 
sién su trayectoria futura. Eso es el determinismo. 


i) detorminiamo absolute dentro de la fisica empezo a ser 
cuestionado a finales del siglo xx por tedricos como Ludwig 
Holwwnwane, que introdujeran el concepto de probabilidad para 
dar une inferpretaci6n mecAnica ala termodinémica. La concep- 
eién de Boltzmann wtilizaba la mecanica cl4sica para describir el 
movimiento de las moléculas, por ejemplo, en los choques entre 
dow moléculas en un gas. Pero aunque estos choques siguieran 
una dindimica estrictamente determinista, la descripcién global 
de un cuerpo macroscépico no dependfa del movimiento deta- 
llado de sus moléculas, sino de promedios estadisticos sobre un 
gran numero de ellas. Y ahi entraba la teoria de probabilidades. 
£1 determinismo sufrié una revisién atin mas profunda con el 
advenimiento de la mecdénica cudntica y, en particular, del prin- 
wiplo de incertidumbre de Heisenberg. Segin Heisenberg, en 
mecanica cudntica no podemos predecir con absoluta certeza 
la posicién de una particula, solo dar la probabilidad de que se 
éncuentre en determinado punto del espacio. Podemos hacer 
predicciones, y de hecho las predicciones de la mecAnica cudn- 
len se encuentran entre las mas precisas nunca realizadas, pero 
no determinar la trayectoria de una particula con absoluta pre- 
cisién. En realidad, el concepto de trayectoria de una particula 
carece de sentido en mecanica cuantica. 

Exceptuando las cuestiones estadisticas suscitadas por Boltz- 
mann, la mecdnica newtoniana, que rige el movimiento de los 
cuerpos macroscépicos y, por tanto, la fisica de nuestro dia a 
dia, parecfa ajena a estos ataques al determinismo. El] primero 
que vislumbr6é que también la mecanica clésica podfa no ser tan 
determinista como se daba por cierto fue el gran matematico y 
Nsico francés Henri Poincaré. En el curso de sus investigacio- 
nes sobre el problema de los tres cuerpos sometidos a su mutua 
atraccion gravitatoria, qué le vallé un premlo convocado por el 
rey Oscar Il de Suecia, Poincaré encontro que hab{a situaciones 
en lag que no era posible dar una solucton cerrada a las trayec- 
torlas de los planetas, Ademis, viskumbré el primer ejemplo de 
trayectoria cadlica, 

Pero no fue haste medindos del piglo xx cuando lo que hoy 
conocemos coma cnos determinisie hive au trrupcién definitiva. 
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Y lo fiige ie o elon alnullinennente en varlog pales y en 
varias disciplinas, Misieos y malemAticos fueron los artifices: de 
fate Lereera revolucion dé Ja fiaien del siglo xx, tras la relatividad 
y la mecanica cudntica, pero también participaron en ella bidlo- 
fos, Quimicos y economistas, y la aparicién de los primeros orde- 
fadores fue una de las claves de esta irrupcién. Los ordenadores 
permitian abordar problemas que los métodos tradicionales de 
ln fisica y las matemAticas no habian podido resolver. 

Il caos determinista se presenta cuando un sistema sometido 
ii leyes o reglas perfectamente determinadas se comporta de ma- 
here, errAtica y, abarentemente, aleatoria. E] comportamiento di- 
fiimico de un sistema cadtico determinista tiene caracteristicas 
diferentes alas de un sistema puramente aleatorio, y a lo largo 
de la segunda mitad del siglo xx los cientificos fueron hallando 
y definiendo estas caracteristicas. La mas importante de todas 
elias es la sensibilidad a las condiciones iniciales: dos estados 
inicialmente muy préximos pueden dar lugar a estados futuros 
completamente separados. Aunque la sensibilidad a las condi- 
ciones iniciales fue prevista por Poincaré, fue el meteorélogo es- 
taclounidense Edward Lorenz el primero en desarrollar un siste- 
iti, cadtico determinista que posefa esta propiedad. Para ilustrar 
fa sensibilidad a las condiciones iniciales, Edward Lorenz acufié 
lz expresion «efecto mariposa», que hoy en dia ha trascendido el 
campo de las matematicas. 

Ademas de la sensibilidad a las condiciones iniciales, los.sis- 
lemas cadticos deterministas exhiben otras propiedades carac- 
leristicas como son la recurrencia, la autosimilaridad y la fracta- 
lidad. La recurrencia se refiere a que el sistema vuelve.una y otra 
vez a estados muy similares al de partida. La autosimilaridad y la 
lractalidad se explicaran a:lo largo de esta obra.. 

El sistema solar y la atmésfera son sistemas cadticos, aunque 
lo son en diferentes .escalas de tiempo. Ambos estan entre los 
primeros sistemas cadticos. estudiados. El sistema solar fue el 
objeto de estudio en que Poincaré hallé por primera vez algunos 
de los elementos del caos, Los ordenadores permitieron, a fina- 
les del siglo pasado, estudiar el movimiento del sistema solar en 
detalle y calcular su futuro a millones de afios vista. Se revelé 
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enitonees que el fitiro del sistema, en ew eseala de Liempo, es 
Here + iipreadecthie, La winoslern tambien ea cadticn y Lo- 
rong (He de ios primeros elontiticos en iratar de cuantificar este 
heels, recurrlendo tambien al caleulo por ordenador. 

NG os necesario que un sistema dependa de un gran mimero 
dé variables para que su comporiamiento pueda ser cadtico. Este 
ee uno de loa descubrinientos mas relevantes de la teoria del 
eno, in el caso de sisternas cuya evolucién depende de forma 
eonlinun del tiempo bastan, de hecho, solo tres variables para 
que ¢l Caos sea posible. Esto permite mostrar las caracteristicas 
osenciales del caos determinista en sistemas relativamente sim- 
ples, como um péndulo o una pelota que bota sobre una mesa que 
viliva, Jon uno de los problemas estudiados por Poincaré, se tra- 
faba de un planeta que orbita entre otros dos. Las tres variables 
involucradas son las dos coordenadas del planeta en el plano de 
In Orbila y su velocidad. El bidlogo tedérico Robert May estudié un 
modelo aun mas simple, con una sola variable. El contexto era la 
Uinfmica de poblaciones y en su modelo la tnica variable era el 
numero de individuos de una especie dada. Una sencilla funcién 
Matematica establecia la poblacién de una generacién conocido 
él numero de individuos de la generacién anterior. May encontr6 
que, bajo determinadas condiciones, el numero de individuos de 
la poblacion podia variar erraticamente de una generacion a otra, 
6s decir, se comportaba de manera cadtica. El fisico Mitchell Fei- 
genbaum estudié el mismo modelo y demostr6 que representaba 
uma evolucién hacia el caos que estaba caracterizada por canti- 
dades universales, en el sentido de que no dependian de la forma 
detallada de la funci6n matematica utilizada en el modelo. 

La evoluci6n en el tiempo de un sistema cualquiera puede re- 
presentarse graficamente. Por ejemplo, podriamos representar 
en un plano todas las orbitas posibles de un planeta en torno 
al Sol, dependiendo de su posicién y velocidad en un momento 
dado. En ese caso podriamos rellenar ¢] plano con todas las tra- 
yectorias posibles. Pues bien, cuando 8¢ representa graficamen- 
te la evolucion de un sistema cadtico aparece un objeto geomé- 
trico muy peculiar: el atractor extrano, Un atractor extrafio es un 
objeto que no Hena completamente ¢l volumen 6 el area que ocu- 
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ph Ne dite que un ataelor oxtranio os (etal Las ties que To 
fecorren, lin ieayeetorins popibles de aiycona que representa, se 
eepnein Y ee cunen lave, Las propledades matemaéticus de jos 
ALICLOPeS OXfPAnos son, Come Veromos, cleramente singtulares:. 

Low alraciores: extrafios Aparecen en elstemas que disipan 
energia y que precisan de un aporle continuo de energia desde el 
estenor. Por ejemplo, el movimiento almostérico toma su ener- 
Win fel calentamiento solar, y la viscosidad del aire se encarga 
ile disipar la energfa en los torbellinos o en el rozamiento con 
sl aiielo, Pero el caos también existe en sistemas que no disipan 
energia, como el propio sistema solar. Ks a una prodigiosa gene- 
fieion de matematicos de la Rusia soviética a quienes debemos 
iH buena parte de la teorfa del caos en los sistemas de este tipo, 
on los que la energia total se conserva, y que reciben el nombre 
wenenco de sistemas conservativos. 

Cuando la teoria del caos determinista empez6 a tomar for- 
ini, Ja mecanica cudntica ya era una teorfa madura. Surgié muy 
pronto la pregunta de como influia el caos en mecdnica cuantica. 
Dicho de otro modo, jqué ocurre en los sistemas cudnticos cuyo 
eorrespondiente sistema clasico es cadtico? La sensibilidad a las 
®ondiciones iniciales, una de las claves del caos determinista, 
it puede estar presente en la mecdnica cudntica, porque, como 
hemos dicho, el concepto de trayectoria no tiene cabida. De al- 
una forma, la mecanica cu4ntica limita o suaviza el caos. Pero 
Ins consecuencias del caos si se dejan sentir en algunas carac- 
leristicas, como la distribucién de los niveles de energia en los 
wislemas cuadnticos cadticos. 

Los conceptos matematicos elaborados para entender el caos 
determinista han encontrado aplicaci6n en otros campos de la 
ciencia, mas alla de la fisica y las matematicas. La biologia, la qui- 
mica e incluso la economia han hallado en las matematicas del 
ci08 herramientas para estudiar sus propios sistemas. Como ha 
demostrado la teoria del caos, un sistema que se comporta de ma- 
nera irregular no tiene por qué ser extremadamente complicado, 
(uizé baste seguir la evolucién de solo unas pocas variables para 
describir ese comportamiento. Esta idea de partida ha motivado 
numerosos estudios en todas estas disciplinas con diversa fortu- 
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no Se han encontrade veaeconesn quimlicas que sluen las leyes 
del cos delerminigtn y su Comportaniento complejo se puede 
oyplicar a partir de unas poons variables, También se han cons- 
truido circultos clectronieos endticos y 8¢ usa él caos para deter- 
minadas tareas en clertos microdispositivos. Pero la existencia 
del caos determinista y su interés en biologia ha quedado patente 
solo on alunos casos. Y mucho mas discutible es su aplicabilidad 
# interés en economia, aunque algunas herramientas matemati- 
ons desarrolladas en el contexto del caos determinista han sido 
Incorporadas plenamente por los economistas. 

Caos es una palabra antigua, de origen griego. Significa «confu- 
#idn, desorden». Su antonimo original es cosmos, también de ori- 
fen griego, que significa «orden». Asi, segin la Biblia, reinaba 
@l Caos antes de que Dios pusiera orden en las cosas, creando el 
cosmos, Pero el uso, el tiempo y, a veces, también la ciencia, van 
modiflcando el sentido de las palabras, ampliandolo, matiz4ndo- 
lo, La @xpresion caos determinista es, en cierta forma, una con- 
iyadiccién en si misma, pero ilustra un tipo de comportamiento 
muy atractivo, equidistante del caos absoluto y del determinismo 
abaeeoluto, Justo donde suceden las cosas interesantes. Es un com- 
portamiento irregular, pero que sigue reglas definidas. Parcial- 
mente predecible, parcialmente impredecible, pero con normas 
y estructuras subyacentes. Coro la vida misma. 
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OAPITULO 1 


La atmoésfera y ei sistema solar: 
Caos en Ia tierra y en el cielo 


iin un sistema cadtico las magnitudes fisicas 
varian de forma erratica en el tiempo. En el 
caos determinista este comportamiento erratico 
es debido a reglas perfectamente definidas, 
determinadas. Dos ejemplos caracteristicos 

de caos determinista son el sistema solar 

y el tiempo atmosférico. 


Cuenta una antigua leyenda china que habia una region Nlama- 
ila (Qi cuyos habitantes se preguntaban si algin dia el cielo cae- 
ria sobre sus cabezas. Les preocupaba que el Sol y las estrellas 
ilejaran de moverse y se precipitaran sobre la Tierra. Esa preo- 
#upacion se extendid por toda la region hasta tal punto que sus 
habitantes empezaron a olvidarse de todo, incluso de las acti- 
vidades mas necesarias, como son dormir y comer. Los qiitas 
se pasaban el dia mirando al cielo, temerosos de que les cayera 
encima. Esta leyenda se usa habitualmente en China para ad- 
vertir de los peligros de las preocupaciones sin sentido y los 
qiitas son vistos, por lo tanto, como un ejemplo de imbecilidad. 

Sin embargo, nuestra civilizacidn debe mucho a los grandes 
personajes que a lo largo de Ja historia se han planteado, de 
una u otra forma, la misma pregunta que los qiitas. Podemos 
recordar aqui a Tomas de Aquino (ca. 1224-1274), el famoso 
traile y fil6sofo medieval. Para él era Dios, el primer motor, 
quien habia puesto en movimiento los objetos celestes, y este 
movimiento era eterno e imperturbable. Si el cielo no se nos 
cafa encima era por obra y gracia de Dios, que asi lo habia 
dispuesto. 
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1) flaice inglés sane Newton establecié au teoria de la grave: 
dad a finales del siglo xvn, La teoria de Newton unificaba en una 
sola ley fisiea dos fenémenos aparentemente distintos: la caida 
de un cuerpo én la Tierra y el movimiento orbital de planetas y 
satélites. La fuerza gravitatoria era la responsable de ambos mo- 
vimmientos. Segin Newton, la Luna «cafa» realmente a la Tierra, 
atraida por la misma fuerza que hace caer a una piedra, pero el 
movimiento de la Luna en direcci6n paralela a la superficie de la 
‘Tierra y la esfericidad de esta hacen que esta caida se transforme 
en un dar vueltas sin fin. Newton fue capaz de demostrar que la 
ley de la gravedad por él propuesta producia las orbitas elipticas 
que Kepler habia observado en los planetas, pero no decia nada 
acerca del porqué de sus movimientos transversales. Si la Luna 
no estuviera animada de un movimiento transversal caeria «a 
plomo» sobre nosotros, tal y como tem{an los qiitas. 

En este punto Newton, hombre profundamente religioso, re- 
currié a Dios, como siglos antes habia hecho Tomas de Aquino. 
Dios, el creador, hnbia animado los planetas y sus satélites con 
gus movimientos iniciales, lo que hacia posible que la Tierra y 
todos log planetas orbitaran alrededor del Sol, y la Luna alre- 
dedor de la Terra, Pero habia algo inquietante: la Tierra no era 
atrafda solo por ¢l Sol, sino también por el resto de los plane- 
tas. Como la masa del Sol es muchisimo mas grande que la de 
cualquier otro planeta, su fuerza es la que manda, y los planetas 
describen sim drbilie alrededor de él. Pero estas 6rbitas no son 
elipses inmutsbles, emo cabria esperar si cada planeta estu- 
viera solo don ¢] 86) én el universo, sino que estan perturbadas 
por las fuerzie de wireccion de los otros planetas. Estas fuerzas 
son muy Pequehas omparadas con la del Sol, casi impercepti- 
bles, pers Newton eamprendié que su accién continuada a lo 
largo de 168 siglo podta desbaratar el engranaje y dar al traste 

con todo. A Malia de aie explicacién mejor, Newton, hijo de una 
Spook ajo a 1 item de evolucién, acudié de nuevo a Dios, al 
que Olds 1 apel se vigilante del buen orden césmico: era 
Dine oi aie, We Vee OH Cuando, reajustaba las 6rbitas para que 
jae pequefine fires entre los distintos cuerpos no rompieran 
el len database, 
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Pero ,cdmo adquirieron li Lina y loa planetas su moviniente 
iraneversal? Que sepamos, fue Immanuel Mant @l primero que 
Witrodujo, un siglo deapudés de Newton, la hipdteais de que el sis 
forma solar habia evolucionado a partir de una nebulesa. En la 
#poea de Kant, el avance en la construccién de telescopios ha- 
bia permitido deseubrir varias nebulosas: objetos no estelares 
6 pspecto difurninado. Estas nebulosas inspiraron al fildsofo 
la hipdtesis, que, en lo esencial, se sigue considerando correcta 
hoy en dfa, de que el sistema solar se originé a partir de una 
fnmensa nube de gas y polvo: la nebulosa solar. Los planetas y 
¢] Sol heredaron sus movimientos de traslacién y rotacién del 
fliro que la nebulosa primordial efectuaba sobre si misma. A medi- 
(la que la nebulosa se contraia y se formaban los planetas, estos 
iban aumentando su velocidad de traslacién, de forma andloga 
fi la aceleracién que experimenta el agua en el torbellino que se 
produce cuando escapa por el desagiie del lavabo. Una de las vir- 
tudes de esta hipdtesis es que explica por qué todos los planetas 
arbitan en el mismo sentido: es el mismo sentido en que giraba la 
nebulosa sobre s{ misma. ;Menos trabajo para Dios! 


ZES INESTABLE EL SISTEMA SOLAR? 


Pierre-Simon de Laplace, astr6nomo y matematico francés de la 
época de la Revolucion, estudié el movimiento de los planetas 
ul detalle, intentando averiguar si las perturbaciones que unos 
@jercian sobre otros, esas que tanto inquietaban a Newton, ter- 
Minarian por desbaratar el orden cdésmico o si, por el contrario, 
fis Orbitas eran estables a largo plazo. Laplace ideé métodos 
de calculo aproximados para encontrar la érbita de un planeta 
perturbado por la presencia de otro y descubrié que los movi- 
Mmientos medios de los planetas eran constantes. Un caso que 
interesaba especialmente alos astrénomos era el de Juipiter y Sa- 
turno. Segtin las medidas de la época, Jupiter se aceleraba poco 
a poco, describiendo su érbita cada vez mas répidamente (hay 
que aclarar que estamos hablando de incrementos mintsculos 
én su velocidad orbital). A su vez, Saturno iba desacelerdndose. 
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Lapliee deiiosted que ealae veriaciolies eran eonsecuencla de 
su abraecion mutue y, lo mie iaportante, que el comporlimniento 
era perigdion: cad 450 anos La situacion se revertin y era Jupiter 
él qué a Trenaba y Satumo el que se aceleraba, A los 900 atos 
ambos volvian a su posicion inital, 

KE] eAleule realizado con Jipiler y Saturno hizo creer a La- 
place que las perturbaciones entre planetas no amenazaban 
It Ostabilidad del sistema solar, y asi lo publicé en su estudio 
que presenté al mismisimo Napoleén con el titulo de Tratado 
ie mecanica celeste, Napoleon era un amante de las ciencias, 
Micionado a las matemiéticas, e impuls6 la creacion de varios 
tentros de ensetianza superior en Francia. Segtin se cuenta, re- 
procho a Laplace que no hablara de Dios en su libro, a lo que 
ésle contest6: «No he necesitado esa hipétesis». Laplace habia 
librado a Dios del trabajo newtoniano de reajustar las 6rbitas 
cada cierto tiempo. 

Pero las cosas no eran tan faciles, y a mediados del siglo xx 
la estabilidad del sistema solar de nuevo fue puesta en cuestion. 
Ii] astr6nomo francés Urbain Le Verrier, famoso por predecir la 
existencia del planeta Neptuno, amplié los caélculos de Laplace 
y se percat6 de que habia perturbaciones mAs pequefias que las 
consideradas por este que tal vez invalidaban los c4lculos. Le 
Verrier se pregunt6 si estas pequefias perturbaciones, en un prin- 
cipio menores, podian causar a largo plazo efectos devastado- 
res, jlora licito no tenerlas en cuenta? 

Asi, el problema fue uno de los temas del premio mateméati- 
co que el rey Oscar II de Suecia convocé con motivo de su 60.° 
cumpleafios. El] matematico sueco Gésta Mittag-Leffler habia 
convencido al rey, muy interesado en las matematicas, para que 
convocara el premio y lo dotara generosamente, de forma 
que atrajera a los mejores matematicos de toda Europa. El pre- 
mio consistia en una medalla de oro conmemorativa y dos mil 
quinientas coronas suecas. E] matematico francés Henri Poinca- 
ré gano el premio con una memoria titulada «Sobre el problema 
de los tres cuerpos y las ecuaciones de la dinamica». Volveremos 
mas adelante a Poincaré y a su memoria ganadora del premio 
del rey Oscar porque esta forma parte de la historia del caos 
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determinista, pero por ahora digamos solo que una de las con 
Chisiones de Polvouré fue que ol problema era, en elerto modo, 
lesoluble, fn Stras pilabris, ala pregunta «{8e nos caeré el 
Hlelo encima algim dia?», ln respuesta de Polnearé era: «No po- 
demoa saberlo con certeza», 

Considerar e? sistema solar coro un todo es un problema 
formidable, aunque lo restrinjamos solo alos planetas y sus sa- 
{elites més grandes. Habria que resolver la trayectoria de dece- 
fia de cuerpos en mutua atraccién gravitatoria. De modo que 
Poincaré no se planted en su trabajo el problema completo, sino 
{itia versién reducida de un problema més sencillo: e] llamado 
problema de los tres cuerpos. El problema de los tres cuerpos 
fe en encontrar las trayectorias de tres cuerpos masivos 
sometidos a mutua atraccion gravitatoria. Pero incluso este pro- 
blema en su forma general, sin poner restricciones a la masa y 
ft las trayectorias de los cuerpos, es demasiado dificil. Asf que 
Poincaré se planted a fondo un problema mas sencillo: el proble- 
ma. de los tres cuerpos restringido. 

insu memoria, Poincaré escribia: 


Imaginemos dos cuerpos, el primero de gran masa, el se- 
gundo de masa finita, pero muy pequefia, y supongamos 
que estos dos cuerpos describen en torno a su centro de 
gravedad comun una circunferencia en un movimiento uni- 
forme. Consideremos ademas un tercer cuerpo, de masa 
infinitamente pequefia, de forma que su movimiento se ve 
afectado por la atraccién de los dos primeros cuerpos, pero 
él no puede afectar la érbita de estos dos primeros cuerpos. 
Limitémonos ademas al caso en que este tercer cuerpo se 
mueve en el plano de las circunferencias descritas por los 
dos primeros cuerpos. Este es el caso de un pequefio pla- 
neta que se mueve bajo la influencia de Jipiter y del Sol 
cuando se desprecia la excentricidad de Jupiter y la incli- 
nacién de las 6rbitas. Este es el caso también de la Luna 
moviéndose bajo la influencia del Sol y de la Tierra cuando 
se desprecia la excentricidad de la Tierra y la inclinacién de 
la érbita de la Luna respecto.de la ecliptica. 
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Tombien seria el ease, por ajenpto, de inn de las ives Apolo 
ena Vvinje de la Terri a le Lana, 

Los modemos ordenadores permilen resolver numeéricamente 
l problema con mucha precision, Los ordenadores no propor 
Monan una formula mhatematica de la trayectoria del planeta, 
fino que, ulilzando algoritmos especialmente disehados, dan 
los valores muméricos de la posicion del planeta en sucesivos 
instantes de tiempo previamente escogidos. Poincaré no tenia 
ordenador, pero si una gran perspicacia y una profunda intui- 
ion matematica, Uno de sus muchos descubrimientos fue darse 
cuenta de que habia dos tipos de soluciones. Por un lado, aque- 
llas que podriamos Hamar regulares: trayectorias normales que 
admiten bien una solucién explicita o bien aproximaciones sufi- 
ientemente buenas. Pero, por otro lado, habia trayectorias que 
No sefulan esa pauta. Y esas trayectorias no regulares impedian 
resolver de manera general el problema de los tres cuerpos y, 
por extension, el del sistema solar como un todo. Hoy en dia 
lunariamos cadticas a estas soluciones. En el caso del sistema 
solar la clave del problema estaba, como Poincaré explicaba en 
su memoria, en las resonancias. 


Las resonancias perturban las orbitas planetarias 


La palabra resonancia describe un gran nimero de fendmenos 
en campos muy diversos, alguno de los cuales incluso ha entra- 
do a formar parte del habla cotidiana, como, por ejemplo, la re- 
sonancia magnética nuclear. La resonancia se produce en siste- 
mas que presentan algiin tipo de comportamiento periddico. Las 
resonancias mecdénicas son las mas sencillas y se dan entre dos 
cuerpos que tienen la misma frecuencia natural de oscilacién. 
Pensemos en una guitarra bien afinada en la que pulsamos una 
cuerda. Esta oscila con una frecuencia que viene determinada 
por el material del que esta hecha, su grosor y su longitud. Si 
acercamos esta guitarra a otra guitarra que esté bien afinada, 
observaremos que esta ultima, que no hemos tocado, también 
suena levemente. El sonido de la cuerda que hemos hecho vi- 
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hior in exellade oseilaeiones eh ln euerdea de Ih guitare que noe 
hemos locwio, Pir ae este fonomeno oeurra, lhe dos euerdis 
Jeben estar bien alinsdas, es decir, vibrar, de manera natural, 
exholamnente Ala misma (recuencia, 

iro fenoémeno resonante al que estamos habituados en nues- 
lr Vida diaria es la sintonizacion de una emisora de radio. Las 
emisoras de radio emiten ondas electromagnéticas que se pro- 
pagan por el aire hasta nuestros receptores. Cada emisora tiene 
iia frecuencia caracteristica de emision, asignada por las autori- 
dades competentes. Inn el caso de la FM estas frecuencias estan 
en el rango de los cientos de megahercios. Pues bien, cuando gi- 
famos e] dial de nuestro receptor estamos ajustando la frecuen- 
ela de oscilacién de un circuito eléctrico interior para que oscile 
exactamente ala misma frecuencia de la emisora que queremos 
oir, (cl circuito entra asf en resonancia con la onda entrante y pro- 
duce una sefal eléctrica que es amplificada, demodulada y trans- 
milida alos altavoces. 

Volvamos al sistema solar. Las resonancias entre dos cuerpos 
#0 dan cada vez que sus periodos orbitales, esto es, el tiempo que 
lardan en completar una vuelta alrededor del Sol, estan en una 
razon que sea un naimero racional, por ejemplo 1:2, 2:3, 3:5, etc. 
Supongamos, para fijar las ideas, que tenemos dos planetas ima- 
winarios, el planeta A y el planeta B, cuyos periodos orbitales son, 
respectivamente, 40 y 60 afios. La raz6n entre ambos periodos es 
ile 2:3. Si partimos de una posicién en la que los dos planetas es- 
lan alineados con el Sol, como en el primer dibujo de la figura 1, al 
cabo de 120 afios el planeta A habra dado tres vueltas completas, 
y el planeta B dos vueltas completas, con lo que ambos se encon- 
(raran de nuevo en la misma posici6n relativa. Supongamos que el 
planeta B es muy masivo y el A es mas bien pequefio. El resultado 
de esta situacién es que cada 120 afios el planeta B tirara leve- 
mente del A alejandolo del Sol. El efecto puede ser muy pequenio, 
porque la atraccion solar es mucho mas importante, pero se pro- 
duce periddicamente, cada 120 afios, con lo que tiene un efecto 
acumulativo que terminara por sacar a A de su orbita. 

Una de las consecuencias mas notables de la existencia de 
las resonancias gravitatorias en el sistema solar es la existen- 
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la del eintaron de nateroldes, Adenia, las resonaiclns (ambiers 
determina la deposleion de los asteroides dentro del cintur6n. 
Venimos como ocurre esto. 

Las distancias dé los planetas al Sol siguen aproximadamente 
una regia, Conocida como ley de Titius-Bode. La figura 2 mues- 
ira la distancia de los planetas al Sol en funcién de su ntimero 
de orden tomando a Mercurio como el primero. Se observa que 
estas distancias siguen una pauta clara, pero entre Marte y Ju- 
piter parece haber un hueco que no ha sido ocupado por ningun 
planeta. 

iste hueco est4 ocupado por un gran nimero de pequefios 
euerpos rocosos, los asteroides. Estos ocupan una zona en la 
que hay varias érbitas que estan en resonancia con la érbita de 
Jupiter. A la distancia de 3,28 UA, la érbita de un astro tiene un 
periodo que es justo la mitad del de Jupiter. Esto significa que 
un objeto que estuviera orbitando a esa distancia del Sol encon- 
traria a Jupiter en la misma posicién cada dos vueltas (resonan- 
cia 2:1), En esta circunstancia la atraccion gravitatoria del pla- 
neta, el segundo objeto mAs masivo del sistema solar, daria un 
pequeno tirén del hipotético astro cada vez que se encontraran, 
y como ja posicion relativa seria siempre la misma, el efecto de 
log tirones se irfa acumulando, hasta sacar al objeto de esa érbi- 
ia y enviarlo a una 6rbita no resonante. 

La grafica de la figura 3 muestra la distribuci6n de los asteroi- 
des en funcién de la distancia al Sol. Vemos que en la posicién de 
la resonancia 2:1 no hay casi asteroides y que hay muy pocos en 
las 6rbitas mas cercanas. La maxima concentracién de asteroi- 
des se produce en puntos alejados de érbitas resonantes. A las 
distancias a las que hay otras resonancias la concentracién de 
asteroides es casi nula. Ello ocurre, por ejemplo, en la distancia 
2,5 UA, en la que el periodo de revolucién alrededor del Sol es un 
tercio del de Jupiter, es decir, un objeto en esa 6rbita encuentra 
a Jipiter en la misma posicién cada tres vueltas alrededor del 
Sol (resonancia 3:1). La existencia de 6rbitas resonantes entre 
Marte y Jupiter ha impedido la formacién de un gran planeta en 
esa area. Es una zona fuertemente perturbada por la presencia 
del mayor planeta del sistema. 
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Posicidn tras 60 afios Posicion tras 120 afios 


Ios planetas en drbitas resonantes con una razén 2:3 entre sus periodos. Se alinean cada 
120 afios. 
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Posicidn de los planetas del sistema solar en funcidn de su numero de arden, Se observa un 
hueco sin planeta entre Marte y Jupiter, justo donde se encuentra el cinturén de astercides. 
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Puede perenne qué ln otielencin de orbitas resoninites es an 
Tenomjens Muy PAeliculir, que sfeeta solo & elerlas posiciomes y 
que no deberia tener inflienets en el stems solar como in todo, 
Pero Polneare se dio cuenta de que el problema de las resonan- 
elas lenin earsecnencias profiundas cuando se trataba de calcu 
iar el comportamiento dé un planeta a muy largo plazo. Como 
hemos dicho, siempre que el cociente de los periodos de dos 
planeias sea. un mimero racional tendremos una resonancia. Las 
resonancias asociadas a numeros como 1:2 0 2:3 son muy fuertes 
y lenen consecuencias evidentes. Pero resonancias como 17:53 
también se haran notar cuando el planeta de periodo mas corto 
haya dado varias veces 53 vueltas. Es decir, es cuestidn de tiem- 
po, de que pase suficiente tiempo. Y el tema de la estabilidad del 
siatema solar también es cuesti6n de tiempo: queremos saber 
qué va a pasar con él en el futuro lejano. Poincaré demostré que, 
puesto que feran cuales fueran los periodos de dos planetas su 
cociente siempre iba a estar cerca de algun naGmero racional, era 
inposible, en general, encontrar una solucion regular, por muy 
altos que fueran los ntimeros que formaran la fracci6n. Dicho en 
lenguaje téenico: Poincaré demostr6 que el desarrollo en serie 
dé las soluciones alas ecuaciones de movimiento no era unifor- 
memente convergente. 

Polncaré respondio ala pregunta de Le Verrier de forma nega- 
tiva: en general, no se pueden despreciar las resonancias peque- 
fas én un estudio a tiempos muy largos. El matematico francés 
replanted completamente los métodos de la mecanica. Demostré 
que no es posible encontrar una soluci6n analitica, una soluci6n 
expresada en términos de funciones matematicas conocidas, 
que represente el movimiento planetario en un intervalo infinito 
de tiempo, También demostré que la consideracién de términos 
cada vez mas pequefios en la interaccién entre los planetas, la 
suma del efecto de mas y mas resonancias, no tenia por qué con- 
ducir a un valor finito. Esto tiene que ver con el problema mate- 
malico de sumar infinitos términos muy pequenos: el resultado 
puede ser un valor finito o infinito. Poincaré demostré que no 
si@mpre la suma de todos los efectos que afectan ala 6rbita de un 
planeta da una 6rbita definida. 
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La consecuencia fisica de estas dificultades matematicas es el 
caos. Determinadas é6rbitas pueden comportarse de forma erra- 
ica en el tiempo, sufriendo variaciones que son, en apariencia, 
indistinguibles de un movimiento aleatorio. Podemos servirnos 
de nuestros modernos ordenadores para explorar.algunas de 
las consecuencias de Ja no regularidad del movimiento de los 
planetas. En la figura 4 se muestran dos trayectorias posibles 
del cuerpo de masa despreciable en el problema restringido de 
los tres cuerpos, el que abord6é Poincaré en su memoria del pre- 
mio del rey de Suecia. Este ejemplo muestra una de las carac- 
teristicas esenciales de los sistemas ca6ticos: la sensibilidad a 
las condiciones iniciales. Las dos trayectorias parten de pun- 
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Hos trayectorias del cuerpo dé masa despreciable en el problema restringldo de los tres cuerpos. 

Lue trayectoriag difloren muy poco on au posicidn iniclal, Las pequefias diferencias en la posicidn inicial 

ao {raducon an trayectorias muy diferentes, Una de las trayectorias rodea el planeta de la izquierda y queda 
aruitando varias vaces al dé la dorecha. La otra trayectoria vuelve al planeta de la izquierda tras una breve 
axcuralén por al de ia derecha. 


a 


tos muy préximos, pero, al cabo de varias vueltas alrededor de 
uno y otro planeta, se separan hasta dar lugar a dos trayectorias 
c omplet ‘amente diferentes. 

Ademas, el problema de los tres cuerpos muestra otra de las 
aefias de identidad del caos determinista: la recurrencia. El cuer- 
po pequefio esta atrapado gravitatoriamente por los dos planetas 
mayores. Aunque su trayectoria sea erratica no se pierde en el 
inflnito, sino que vuelve una y otra vez a dar la vuelta alrede- 
dor de uno u otro de los centros de atraccion. Sensibilidad a las 
condiciones iniciales y recurrencia son dos de los ingredientes 
esenciales del caos determinista. 
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Pl sistema solar 68 oadtion 


Ml adveniniento de los ordenadores a mediados del siglo xx sig- 
Hilles una revolucion en todas las disciplinas cientificas, Gracias 
i allow fue posible resolver de forma numérica sistemas de ecua- 
Hones que hasta entonces habfan sido inmanejables, ya fuera por 
mi complejidad o, sitnplemente, por el ntimero de ecuaciones. L6- 
ficumente los astrénomos, como todos los demas cientificos, co- 
jenzaron a utilizar la nueva herramienta para abordar problemas 
que hasta entonces habian sido inabordables. El comportamiento 
del sistema solar como un todo era uno de esos problemas. 

itt Bureau des Longitudes se creé durante la Revolucion fran- 
fosa, el 7 mesidor del afio III (25 de junio de 1795). Su objetivo, 
fal declarado por la Convencion Nacional, era arrebatar el con- 
trol de los mares a los ingleses. Para ello era necesario mejorar 
fiolablemente el cAlculo de las longitudes en alta mar, un calculo 
de vital importancia para la navegacién segura y eficiente. Entre 
fos fundadores del organismo estaban Lagrange y Laplace. Un 
Algo después de su creacién, Henri Poincaré tuvo el honor de 
presidirlo. 

Desde tiempos inmemoriales el ser humano se ha servido de la 
fistronomia para orientarse, tanto en tierra como en el mar. Las 
@atrellas siempre han guiado a los marinos, y no pocos avances 
clentificos estan ligados a la navegacion. Si el cdlculo de la latitud 
(le un lugar a partir de observaciones astronémicas es relativa- 
mente facil (basta observar la altura de la estrella Polar sobre el 
horizonte o medir la altura del Sol a mediodia), conocer la longi- 
fides un problema mucho mas complejo. Su calculo se obtiene 
midiendo la diferencia horaria entre un lugar conocido (el puerto 
de partida, por ejemplo) y la posicién en la que el navio se en- 
cuentra en un momento determinado. El problema no se resolvid 
completamente hasta que el inglés John Harrison inventé a me- 
diados del siglo xvmi su cronémetro marino, que funcionaba con 
precision.incluso en las tormentas mas severas, como demostr6é 
él capitan Cook en uno de sus viajes alrededor del mundo. 

El conocimiento preciso de las efemérides astrondémicas, las 
coordenadas de los planetas y las estrellas en el cielo para cada 
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Soluci6n exacta * 


Solucion a la ecuacidn 
diferencial dy/dt= y. La linea 
continua representa la solucién 
exacta. Los circulos, una 
solucién numérica. 


LAATMOSFERA Y EL SISTEMA SOLAR: CAOS EN LA TIERRA Y EN Bl. GIELO 


Son 


LA ATMOSFERA Y EL SISTEMA SOLAR: CAOS EN LA TIERRA Y EN EL CIELO 


Mo 


dia, tarnibién permite deduelr la longinad del lugar de observacion, 
Y el Bureau dea Longi(udes fue creado eon Ia express mision de 
olaborar tadlas 16 suflelentemente precisas como para permitir 
alos marinos franceses caleniar su posicién a partir de la obser- 
vacion de la posicicn de los astros en el cielo, Desde entonces la 
Oficina ha asumido otras tareas relacionadas, como son la orga- 
nizacion de expediciones cientificas y la participacién en las con- 
ferencias internacionales que establecen los patrones de tiempo. 

Un astronomo del Bureau des Longitudes, Jacques Laskar, em- 
prendié en los afios noventa del siglo pasado la formidable tarea 
dé calcular las efemérides astronédmicas «para siempre». Laskar 
no pretendia, evidentemente, proporcionar tablas de navegacién 
para los marineros del siglo cinco mil. Su trabajo estaba en reali- 
dad en la estela de los de Lagrange, Laplace y Poincaré: responder 
de una vez por todas al tema que tanto preocupaba a los qiitas. 

Lo primero que a uno se le ocurre si quiere calcular las trayec- 
(orias de todos Jos planetas del sistema solar con ayuda de un 
ordenador es poner en la ecuacién de cada planeta la fuerza de 
fravedad que, segtn la ley de Newton, ejercen el Sol y los demas 
planetas sobre é]. Es decir, meter todas las ecuaciones sin apro- 
ximaciones y sin simplificaciones. Un grupo de cientificos del 
Jet Propulsion Laboratory, en Estados Unidos, hicieron eso en la 
década de 1980 y «solo» pudieron calcular el futuro del sistema 
para los préximos 44 siglos. Laskar queria mAs y, siguiendo los 
méLodos de aproximacion ideados por sus ilustres predecesores 
(Laplace, Poincaré y De Verrier, entre otros), pergefié un modelo 
que tomaba en cuenta los ocho planetas principales incluyendo 
los efectos conocidos mas importantes y despreciando términos 
irrelevantes. Las ecuaciones aproximadas de Laskar contenian 
150000 términos polinomiales y ocupaban cerca de ;1000 folios! 
Eso le permitié ir mucho mas all4 que sus colegas estadouniden- 
ses y pudo calcular con fiabilidad las trayectorias de los planetas 
hasta dentro de 200 millones de afios. Con estas ecuaciones Las- 
kar demostr6 que e] sistema solar es cadtico. ,Qué queremos de- 
cir cuando decimos cadtico? Que las trayectorias, en esas escalas 
de tiempo, se muestran irregulares, como si fueran aleatorias. En 
particular, si calculamos las trayectorias de un mismo planeta, por 
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Qenplo la herr, pariende de dos poslelanes iy aliilares, al 
eabo de LOO millones de itos lis trayectorias serin complelamen 
fe diferentes, In estns eorcunstancias no lene aentido eulerar las 
efemerides: mas alla de lox 100 millones de afios, porque esta sen- 
fibilidad a las condiciones iniciales las hace inexactas. Pero eso 
HO significa que no podamos sacar conclusiones, Lo que los fisicos 
¥ malematicos hacen en estos casos es calcular todo un conjunto 
de trayectorias partiendo de condiciones iniciales diferentes. 
Para calcular la trayectoria de un planeta que se mueve en 
#1 espacio tridimensional necesitamos dar seis nimeros inicial- 
thente: tres para especificar su posicién en el espacio y otros tres 
para su velocidad. Estos conjuntos de ntimeros forman un espa- 
elo abstracto que en fisica matematica recibe el nombre de espa- 
o1o de las fases. Fl espacio de las fases tiene dimensién 2 para un 
movimiento en una linea, dimensién 4 para un movimiento en el 
plano y dimensién 6 para un movimiento en el espacio. 
ara ilustrar el concepto consideremos un cuerpo que esté 
res(ringido a moverse sobre una l{nea, por ejemplo un skater que 
se lira por una pista estrecha de skateboard. En este caso, basta 
Conocer su posicién y su velocidad inicial para saber cuanto tax 
dara en recorrer la pista. En el espacio de las fases el movimien- 
to del skater queda especificado si representamos la velocidad, 
v, en cada instante en funcion de la posicién, x, como vemos en 
la figura 5. El movimiento de vaivén queda reflejado en el espa- 
cio de las fases en una curva en forma de elipse: la velocidad 
es maxima cuando el patinador pasa por el centro de la pista 
(v=0), su punto mas bajo, y cero cuando Jlega a uno de los extre- 
mos (% maximo o minimo), donde la altura es maxima. Si el ska- 
ter consigue mas altura en su ejercicio, la elipse correspondiente 
iene una amplitud mayor. Si el movimiento se produce en una 
pista muy ancha, de forma que el patinador pueda desplazarse 
por su superficie, necesitaremos dos nimeros para especificar 
Su posicion inicial y dos nimeros para especificar su velocidad. 
Isl espacio de las fases tiene ahora dimension 4. Y si la pista tiene 
rampas y monticulos que permitan al skater saltar y hacer acro- 
bacias en el aire, necesitamos dos nimeros mas para especificar 
su altura y su velocidad vertical. El espacio de las fases tiene 
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Chater \iniridoae por unt plsta y espacio de las fases. En la trayectoria interior se deja caer desde una altura mas 
hile y alenn2n ink menor velosldad maxima. 


entonces seis dimensiones, como el espacio de las fases de cual- 
quler particula que se mueva en el espacio al que nos referiamos 
on el parrafo anterior. 

hi en lugar de un cuerpo tenemos N cuerpos, el espacio de las 
fases tiene dimensién 6N, para un movimiento tridimensional. 
EI espacio de las fases del sistema solar tiene, por tanto, una 
dimension muy graride, y para realizar una integracién de sus 
ecuaciones, esto es, para resolverlas, hay que dar un numero de 
valores iniciales igual a seis veces e] numero de planetas y saté- 
lites que incluyamos en el calculo. 

Lo que hizo Laskar fue explorar el espacio de las fases, es de- 
cir, estudiar qué zonas del espacio de las fases eran accesibles 
desde los valores actuales de las posiciones y velocidades de los 
planetas o valores cercanos a ellos. La conclusién a la que Ille- 
g6 fue que los planetas interiores (Mercurio, Venus, la Tierra y 
Marte) estan sometidos a variaciones erraticas de sus elementos 
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orbitales, Por ejemplo, la exeentrieidad de Marte, que refleja lo 
Alargacds que es la elipae que deseribe, varia cadticamente entre 
0, drbita cireulay, y 0,1 (un valor cereano al actual), La excentrl- 
eidad de la ‘Tera varia entre 0 y 0,06 (el valor actual es 0,017). 
La causa principal dé estas varlaclones errdticas son las reso- 
MiNielas entre los movinuentos orbitales de la Tierra y Marte. La 
parte exterior del sistema solar, de Jupiter hacia fuera, tiene una 
dinfimica menos cadtica. 

Como todos hemos experimentado, la capacidad de célculo 
de los ordenadores avanza sin cesar y hoy en dia es posible Ile- 
war Jos caleulos mucho més lejos que en en la década de 1990, 
Recientemente, Laskar y sus colegas del Bureau des Longitudes 
han utilizado supercomputadoras para calcular la evoluci6n del 
sistema solar, sin aproximaciones, hasta dentro de 5000 millo- 
nes de afios. Mas allé de este intervalo no tiene sentido hacer 
efilculos, puesto que ese es el tiempo de actividad que se espera 
para nuestro Sol. A partir de entonces seré la evolucién del Sol 
eomo estrella la que marque el destino de todos sus planetas., 

Estos calculos han confirmado todo lo dicho anteriormente y 
muestran, ademas, las siguientes conclusiones: 


-— La mayoria de las érbitas accesibles desde la situacién ac- 
tual no prevén colisiones entre planetas. En el 1% de los 
casos la érbita de Mercurio podria deformarse lo suficiente 
como para que el planeta fuera engullido por el Sol o cho- 
cara con Venus. En un escenario muy improbable, la excen- 
tricidad de Marte podria incrementarse lo suficiente como 
para pasar cerca de la Tierra o incluso colisionar con ella. 
Ello desestabilizaria las 6rbitas de los planetas terrestres y 
podria provocar una colisién entre la Tierra y Venus. 


— Es posible, por fin, responder a la pregunta que Newton se 
hacia hace mas de 300 afios: ;pueden las érbitas planetarias 
desestabilizarse lo suficiente como para que dos planetas 
choquen entre si o alguno sea eyectado del sistema? La res- 
puesta es que, con un horizonte de 5000 millones de afios, 
estos sucesos catastréficos tienen una probabilidad del 1%. 


LA ATMOSFERA Y EL SISTEMA SOLAR; CAOS EN LA TIERRA Y EN EL CIELO 


33 


‘Tmpoes hay que procoupare denuwdide por estas concluslo 


nee Joe preos miles de files 60 upirain beta franquilos. 


LSE PUEDE PREDECIR EL TIEMPO ATMOSFERICO? 


In torne a 19.1 Lewis Fry Richardson comenzo a barruntar la 
idea de predecir el iernpo no de la manera en que se habia reali- 
yado hasta entonees, sino calculandoto. fin su época, y durante 
muchas décadas mas tarde, la prediccién del tiempo atmosféri- 
co se basaba en algunos conceptos generales y en la experien- 
cia del meteordlogo. Los meteorélogos recopilaban datos de 
las distinias estaciones meteorolégicas y dibujaban con ellos, 
manualmente, los Correspondientes mapas de presion, vientos, 
temperaturas y cualquier otra variable que consideraran de inte- 
rés. A partir de estos mapas trataban de inferir como evoluciona- 
ria In utmésfera, 0 buscaban situaciones similares en el pasado 
y hacian su prediccién baséndose en lo que habia ocurrido otras 
VOCCUN 

Parn Richardson esta situacién era poco satisfactoria. Ponia 
camo ejemplo lo que ocurre en astronomia: los astronomos cal- 
culan las efemérides, las predicen a partir del conocimiento de- 
iillado de las érbitas de los planetas y su movimiento. En sus 
propias palabras, escritas en el prdlogo de su obra Prediccién 
del tiempo mediante procesos numéricos: 


Pero el almanaque nautico, esa maravilla de prediccion pre- 
cisa, no se basa en el principio de que la historia de la as- 
tronomia se repite una y otra vez. Sc puede asegurar con 
bastante certeza que una particular disposicion de estrellas, 
planetas y satélites nunca se repite dos veces. {Por qué fba- 
mos a esperar entonces que el mapa meteoroldgico de hoy 
estuviera exactamente recogido en un catalogo de mapas 
anteriores? 


En contraposicién al método habitual, Richardson proponia 
«calcular» el tiempo atmosférico: resolver las ecuaciones fisicas 
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que deter ta evolaeion de la 
velooliid del viento, ln presion o la 


feniperaturn, Pero estas ¢cuneiones 


son enormemente complejas. Por 


“jemplo, la lemperatura del aire en 
tin punto y un instante dados depen: 


Jal vez algun dia en el borroso 
luluiro Sea posible hacer los 
calculos mas rapido que el 
avance del tiempo atmosférico 
y aun coste menor que lo que 


de de Ja temperatura en un instante se ahorre la humanidad debido 
anterior y de la temperatura en los q la informacion conseguida. 


puntos circundantes. También de- Pero eso es un suefio 
pende de la densidad del aire y de su 


velocidad en los puntos de alrededor. 
Todo ello se convierte en un conjunto 
ile ecuaciones diferenciales que hay que resolver. 

Para aplicar su técnica de calculo, Richardson reemplaz6 las 
disiribuciones continuas en el espacio y el tiempo por una red 
discreta. Dividié la parte de Europa que su modelo describia en 
una malla de 25 cuadrados, cada uno de ellos de 200 km de lado. 
Verticalmente dividié la atmésfera en cinco capas. En total, te- 
Nia, por tanto, 125 celdas tridimensionales. Dentro de cada una 
de esas celdas, las variables se suponia que tomaban un valor 
uniforme. Por ejemplo, se consideraba que la temperatura era la 
misma en todos los puntos de una celda dada. 

Richardson también tenia que discretizar el tiempo y para ello 
considero intervalos de seis horas. En sus ecuaciones reempla- 
zaba las derivadas temporales de las variables atmosféricas por 
simples diferencias. Por ejemplo, la cantidad d7/dt, la derivada 
de la temperatura con respecto al tiempo, era reemplazada por 
la siguiente expresion: 


donde T,,,, es la temperatura en el instante que se quiere cal- 
cular, 7’, es la temperatura conocida en el instante t y At es el 
intervalo de tiempo que transcurre entre la situacion conocida y 


la que se quiere calcular. 
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Partiends dé low datos Gheervados en ung parte de Europa el 
20 de mayo de 1910, Richardson pretendia ealeufar la prosian y 
la velocidad del viento sels horas mis tarde en un punto escogi- 
do del centro de Buropa. Evidentemente, no se trataba de prede- 
cir el tiempo ese dfa, porque, entre otras cosas, ese dia ya habia 
pasado, sino de poner de manitiesto que el tiempo atmosférico 
era «calculable». Richardson empezé a trabajar en serio en su 
proyecto en 1918, en el Observatorio Eskdalemuir, en Escocia. 

Richardson pertenecta a la comunidad cuaéquera, una sec- 
ta cristiana que rechaza la guerra, y cuando estall6 la Primera 
Guerra Mundial se declaré objetor de conciencia, por lo que 
no ingres6 en el ejército. En cambio, se alisto voluntario como 
conductor de ambulancias y fue destinado a Francia. Fue en los 
aiios 1916 y 1917 cuando Richardson formulé el problema de 
forma completa. Aprovechaba para ello el tiempo muerto que 
transcurria entre servicio y servicio de su ambulancia. Segtin sus 
propias palabras: «Mi oficina era un monton de paja en un frio 
barracon militar». Con todos los datos disponibles y todas las 
ecuaciones preparadas, en seis semanas pudo hacer el calculo 
completo para el punto que habia escogido como ejemplo. Un 
papel, un lapiz, unas plantillas hechas por él mismo y unas tablas 
dé logaritmos eran todo su equipo. 

Cuando Richardson terminé sus cuentas el resultado fue un 
rotundo fracaso. El valor que obtuvo para la presién estaba muy 
alejado del que se habia observado aquel dia. La diferencia no 
era atribuible a un error en sus cAlculos, sino que tenia una ra- 
z6n mas profunda y que desvelaremos a lo largo de este libro a 
su debido tiempo. Baste decir por ahora que la dificultad tiene 
que ver con los valores iniciales utilizados, algo que el propio 
Richardson ya intuy6. 

Richardson recogié todos los detalles de su aventura aritmé- 
tica en la obra mencionada. Aunque no se habia visto culminada 
con éxito, su experiencia arrojaba interesantes lecciones para el 
futuro. El, en cualquier caso, se mostraba optimista. Segtin sus 
cuentas, 64000 personas operando simultaneamente podrian cal- 
cular el tiempo atmosférico de todo el planeta lo suficientemen- 
te rapido como para predecir los acontecimientos atmosféricos 
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antes de que tuviernn lugar, Hielardson sefabe eon una inmensa 
sala len de personas caleulands y remitiendo sus informes a 
Halo otlenistas jefes, los cuales distripulrian por (elégrafo el 
(lempe predicho a los observatorios de todo el mundo. 


La era de los ordenadores 


16 que Richardson no podia imaginar era que las 64000 perso- 
fia (ue él presuponia necesarias no iban a hacer falta: un in- 
fenio automatico, el ordenador, haria los calculos por ellas. Al 
(ryminar la Segunda Guerra Mundial, el matematico John von 
Neumann, miembro permanente del Instituto de Estudios Avan- 
gados de Princeton, en Estados Unidos, se interesé por el cal- 
fulo automatico. Von Neumann dirigié el equipo que realiz6 la 
primera resolucién numérica por ordenador de las ecuaciones 
del tiempo atmosférico. Von Neumann y Jule Charney utilizaron 
para ello el primer computador electrénico de propésito general 
(le la historia, el ENIAC (Electronic Numerical Integrator and 
Computer). El ENIAC habia sido desarrollado por cientificos 
dt la Universidad de Pensilvania a finales de la Segunda Guerra 
Mundial y fue financiado por el ejército de Estados Unidos con la 
intencién de utilizarlo en los calculos de artilleria. Von Neumann 
trabajaba entonces en Los Alamos en el Proyecto Manhattan, el 
proyecto de disefio y construccion de la primera bomba atomi- 
ca, Von Neumann supo de la construccién del computador y co- 
laboré activamente con sus creadores. De hecho, los primeros 
cdlculos que se realizaron en el ENIAC tuvieron que ver con el 
Froyecto Manhattan y fueron propuestos por él. 

Von Neumann y sus colaboradores publicaron sus primeros 
resultados sobre la prediccién del tiempo atmosférico en 1950. 
Dividian Estados Unidos en una malla cuadrada de 15x18 y el 
lado de cada una de las celdas resultantes tenia 736 km. En estas 
condiciones no cabia esperar una prediccién detallada. Ademas, 
con esta discretizacién, el ENIAC tardaba unas 24 horas en cal- 
cular la evolucion de la atmésfera a 24 horas vista, por lo que no 
hubiera servido como predictor. Pero los cientificos esperaban 
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True invite oT 


faciociones meteoralégicas fiables. 

modelo Incorpora no solo la 

spe elevantes, Gomo la forma de 
per iss oaraotersioas ¢ eseficiales de fa 


UNOS ( 


Jo), qué podriamos comparar con las 60 operaciones. 
@f SUS primeros trabajos, Los dos supercom- 

y son alitvientades de manera independiente con las” 
in frente a fallos eléctricos y de otro tino. De esta forma se 
WE proparciona cada 12 horas predicciones para 
roporciona predicciones, con menor resolucion, a 


SRAY-1, el primer superorcenador adquirido por el ECMWE, en 1977. 
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qe Un neve demputador, que ya OatAbnen eenatrucelon én el 
intariisime instiiute de Betudios Avanzados cor una arquilec- 
jure interna disefiada por Von Neumann, pudiern realizar los 
cAleulos en menos de media hora y con una malla mas tupida. 
Los resultacios obtenidos no predecian el tiempo correctamente, 
pero sf algunas caracteristicas generales observadas. Id] suefio 
de Richardson estaba mas cerca. 


Edward Lorenz y una atmdsfera de juguete 


Y aguf entra en escena el meteordlogo estadounidense Edward 
Lorenz. Lorenz habia servido como meteordélogo de campo del 
ejército del aire y, al terminar la guerra, se incorporé al Departa- 
mento de Meteorologfa del Instituto Tecnolégico de Massachu- 
setts, el afamado MIT. Lorenz adquirié allf una sdélida formaci6n 
en meteorologia dindmica, la que trata de entender la dinamica 
atmosférica a partir de las ecuaciones que rigen los procesos at- 
mos{éricos, y empezé a interesarse por la prediccién numérica 
tal como la realizaban Charney y Von Neumann. 

in esos primeros afios de los ordenadores, algunos meteor6- 
logos se hab{an inclinado por utilizarlos con un enfoque estadis- 
tico, mas parecido a la prediccién sinéptica habitual: situaciones 
semejantes daran lugar a comportamientos, semejantes. Los es- 
tadisticos utilizaban los ordenadores para manejar gran canti- 
dad de datos y tratar de hacer pronésticos a partir del compor- 
tamiento pasado de la atmésfera. Era, digamoslo asi, una suerte 
de extrapolacion. 

En torno a 1956 Lorenz aaa una forma de discernir si-el en- 
foque estadistico y el dindmico podian ser equivalentes. Partiria 
de un sistema de ecuaciones simplificado, pero cuyas soluciones 
tuvieran suficiente variabilidad. Esta variabilidad era necesaria 
porque las soluciones estacionarias, periddicas o muy. repetiti- 
vas son facilmente predecibles. Y utilizaria un ordenador para 
resolver las ecuaciones y generar una solucién numérica sufi- 
cientemente larga, que permitiera valorar la.capacidad de pre- 
diccién a medio y largo plazo. Tratando la solucién como si fuera 
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ne Coleccion dé datos meteorologicas reales, Lorenz utilizaria 
procedimientos ostadisiicos estandar park predeciy el compor- 
iimlento de la serie numérica, la secuencia de valores obteni- 
dos, Sila prediccién estadistica se correspond{a con la solucién 
numérica verdadera, la validez del método estadistico quedarfa 
ostablecida. 

Por sugerencia de un colega, Lorenz propuso al departamento 
la adquisici6én de un ordenador para su proyecto. A mediados 
Ue la década de 1950 los ordenadores no eran, ni mucho menos, 
algo corriente. El ordenador que adquiri6 Lorenz era un Royal 
McBee LGP-30 con una memoria interna de unos 16 kilobytes y 
que realizaba unas 60 multiplicaciones por segundo. La primera 
tarea a la que tuvo que dedicarse Lorenz fue la de programar. 
Muchas veces se nos escapa la gran cantidad de horas de traba- 
jo abnegado y rutinario que conileva a menudo la investigaci6én 
cientifica. En una época en que la programacisn informatica es- 
taba dando sus primeros pasos, pocas personas sabian progra- 
mar un ordenador y la tarea era compleja e intrincada. Lorenz 
paso largos meses aprendiendo a utilizar su ordenador antes de 
poder hacer ningun célculo efectivo con él. 

Inspirado por las ecuaciones reales de la atmésfera, Lorenz 
ideo un sistema de unas pocas ecuaciones, doce en la versién 
final, que constitufan una especie de atmédsfera de juguete: ha- 
bia calentamiento solar, ciclones y anticiclones, vientos. del este 
y del oeste, y poco mas, Pero, de nuevo tras muchos ensayos 
fallidos y muchas horas de trabajo, consiguié un rango de para- 
metros en los que su atmésfera mostraba el comportamiento es- 
perado. La grafica de la figura 6 muestra una recreacién de cémo 
evolucionaba una de las variables de la atmésfera de Lorenz du- 
rante cinco «meses». El comportamiento no es completamente 
periédico, pero tampoco completamente aleatorio. 

El ordenador arrojaba sus datos por medio de una impreso- 
ra, todavia estaba lejos la invencién de los terminales grdaficos. 
Lorenz y sus colegas de departamento se agrupaban a veces en 
torno a la impresora para ver salir los resultados. Incluso ha- 
clan pequefias apuestas sobre el comportamiento de la siguien- 
te ristra de numeros, igual que hacen a veces los meteordélogos 


LA ATMOSFERA Y EL SISTEMA SOLAR: CAOS EN LA TIERRAY EN EL CIELO 


41 


eon el (lenge del dia aigutente, Pronto se perentaron ce alu 
nos Comporiaiionios peculiaras, de patrones que se reprodu 
elan con mide feeweneia que otyes, [én clerte mode jugaben a 
sor meteordlogos de aquella «almoésfera numercns, Lorenz pudo 
comprobar en su atmosfera imaginaria que los procedimientos 
caladistieos habituales eran incapaces de predecir lo que iba a 
ocurrir, 8] comportamiento de su atmésfera mimetizaba el de la 
almosfera real. 

Un buen dia Lorenz decidié repetir algunos caélculos, con el 
lin de examinar con detalle una de las variables. Detuvo el or- 
(denador, tecleo una linea de nimeros que habia salido por la 
impresora un rato antes y volvi6 a ponerlo en marcha. Se fue a 
tomar un calé y volvid al cabo de una hora (debia de ser un café 
iimericano), Cuando regres6 comprob6 que el ordenador habia 
simulado unos dos meses de tiempo meteoroldgico, pero los nt- 
ievos que Salian por la inypresora no tenian nada que ver con los 
que habla calculado antes. 

Lo primero que le vino a la cabeza era que el ordenador se ha- 
bia estropeado. Pero antes de llamar a nadie repasé los numeros 
que habian estado saliendo por la impresora durante su ausen- 
cia. Mnseguida se dio cuenta de que los nuevos valores repetian 
los anteriores al principio, pero pronto empezaban a diferir. Las 
diferencias aumentaban duplicandose cada cierto tiempo, hasta 
(ue las cifras, al final, diferian completamente. La grafica de la 
figura 7 muestra lo que Lorenz observé. 

Lorenz comprendioé lo que habia ocurrido: los nimeros que él 
habia tecleado en la computadora no eran los originales exactos. 
Su maquina calculaba con seis cifras significativas, pero la impre- 
sora solo imprimia tres, para ahorrar espacio. Por lo tanto, los nt- 
meros introducidos diferian de los originales a partir de la cuarta 
cifra. Lorenz habia descubierto la sensibilidad a las condiciones 
iniciales, el rasgo mAs caracteristico del caos determinista. 

Si una atmosfera tan simplificada ya mostraba este compor- 
tamiento, qué cabria esperar de la atmésfera verdadera? Las 
magnitudes fisicas en un observatorio meteorolégico rara vez 
se miden con una precision de tres cifras. Y la distribucién de 
observatorios es irregular y obliga a hacer interpolaciones que, 
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Recreacion de los resultados del modelo simplificado de atmdsfera utilizado por Lorenz. La variable indica la 
latitud de ios vientos de poniente mas fuertes durante cinco «meses». 


FIG. 7 


Grafica obtenida por Lorenz en su modelo simplificado de la atmésfera, en la que se aprecia la sensibilidad a las 
condiciones iniciales. 


en algunos casos, tienen mucho margen de error. Por lo tanto, 
Lorenz hab{fa descubierto que la predicci6n meteoroldégica a lar- 
go plazo era una quimera. 

En la actualidad los modelos numéricos son la herramienta 
basica con la que se realizan las predicciones meteoroldégicas. 
Estas son bastante flables a dos o tres dias vista, incluso a una 
semana, pero los modelos son incapaces de proporcionar un 
pronéstico aceptable con mas de diez dias de antelacidn. 
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es ied Cuda ih diferentes 
_ de caos existen? Trataremos de responder a estos interrogantes 
~ en los siguientes capitulos. 


in los afios setenta y ochenta del siglo pasado 

fue surgiendo entre cientificos de varias 

disciplinas la conciencia de que la aleatoriedad 

aparente que mostraban algunos fen6émenos 

se podia explicar mediante leyes relativamente 
encillas. 
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Entre los afios 1950 y 2015, la poblacié6n mundial pasé de 2500 
millones a 7370 millones de personas. Dicha evolucién se reco- 
en la grafica de la figura 1. La linea continua trazada sobre la 
prafica corresponde a un incremento exponencial de un 9% cada 
cinco afios. 
Esta evolucidon es el resultado de que la poblacién, cada cinco 
_afios, se multiplica por un factor constante, lo que se refleja en la 
siguiente ley matematica: 


donde x, es la poblacién en un determinado afio, x, es la po- 

blacién cinco afios mas tarde y & es una constante, que viene de- 

terminada por la tasa de crecimiento. En este caso, o vale apro- 
ximadamente 1,09. Si la poblacién en 1960 era de 3000 millones 
de habitantes, en 1965 se acercaba a 1,09-3000=3270 millones. 

m Al multiplicarse por un factor constante cada vez, los va- 

lores de la poblacion cada cinco afios siguen una progresién 
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jleometrion y fe iepliran de fork exponeneial Male erect 
iene exponciew de la poblacian human fue predicho por 
@) penendar ingles Thomas Malihus on 1798, Malthiis aiuguraba 
que la poblaeian, en aueencia de obelaculos, ereria seein una 
prosresion peametrien, de ln forma ce, @", ce’, @... mientras 
que los recursos numentarian ¢n progresion arilmética, a la 
manner de la sucesion t, 2, 4, 4,,.. Como la progresién geomé- 
Ines ereee mucho mas deprisa que la arilmeética, era solo cues- 
ion de lempo que la poblacién humana empezara a toparse 
con difieultades sino se lomaban medidas para controlar el 
erecimuento. 

in ausencia de factores limitantes, Malthus predecia que la 
pablacion se doblaria cada 25 afios. En la grafica de la figura 1 
VeINOSs que eso no ocurre en la actualidad (se necesitan cerca 
de 40 anos para duplicar la poblacién mundial), pero lo que si 
ex cierto, en el intervalo analizado en esta grafica, es que el au- 
mento de la poblacién sigue una progresidn geométrica, con un 


miuvemento casi constante del 9% cada cinco afios. Es decir, se. 


nmple la ley de Malthus pero el coeficiente @ es mas pequefio 
de lo que éf calceulé. 

La ley de Malthus es un caso particular de lo que en mate- 
miticns se conoce como un sistema dindmico discreto. En ge- 
neal, un sistema dindmico es cualquier sistema cuyo estado 
eanibia en el tiempo. La palabra sistema la debemos entender 
jit en un sentido completamente general: puede tratarse de 
un conjunto de cuerpos celestes, el cuerpo humano, un ecosis- 
tema, la cotizacion de las acciones en la bolsa, etcétera. Para 
estudiar la evolucion de un sistema dinamico se precisa una ley 
que nos diga c6mo cambia con el tiempo. El sistema es discreto 
si esta ley proporciona el estado del sistema en intervalos fijos 
de tiempo, por ejemplo, el valor de las acciones cada dia, o el de 
una poblacion cada generacion o cada cierto numero de afios, 
como en la figura 1. Si la ley proporciona el estado del sistema 
de instante en instante, para cualquier tiempo, se dice que el 
sistema es continuo. 

Un sistema dinamico discreto viene definido por una variable 
x, que representa el estado del sistema (el numero de individuos 
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con un incremento de! 9% cada cinco afos. 


dle la poblacién mundial cada cinco afios, en el ejemplo ante- 
rior), y una ecuacion que dice como evoluciona el sistema: 


x, =F) n=l, 2, By eee 


La ley de evolucién viene representada por f(z), una funcién 
matematica. Lo que expresa la ecuacion anterior es que en cada 
iteraci6n obtenemos el nuevo valor de x a partir del valor an- 
terior. Para un valor inicial x,, la aplicacion de f(x) (primera 
iteracion) produce x, =f(%,). Una segunda aplicacion, esta vez a 
“,, produce x,=f(x,) (segunda iteracién); haciéndolo una vez 
mas obtenemos x, = f(%,). De esta forma, la aplicaci6n reiterada 
de f(x) produce una sucesi6on de valores %,, %,, Xp)... a partir 
de un valor inicial x, dado. Por ejemplo, podemos construir un 


{QUE ES EL CAOS? 


20 


20 


Evolucin de la poblacién humana del planeta entre 1950 y 2015. El aumento sigue una progresién geométrica, 
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sleternn dindmioo seneilla duplleando cada pase el valer de 
© Le finelsn sere entonces 7(a)e ee, y vn Vator inletal de ja 
variable av, « 2 produciria la sucesién de valores 2, 4, 8, 16, 82,... 
Dando valores inicliles diferentes obtendriamos aucesiones di- 
ferentes, Cada una correspondiente a una evolucién posible del 
sistema, La figura 2 representa la funcion f(2) @ 2x, asf como la 
evolucién a la que da lugar su aplicacién reiterada a partir del 
valor inicial 2, 

En la teorfa de Malthus f(a) es una funcién lineal: f@v) = az, 
y su representacion en una grafica es una linea recta similar a la 
de la figura 2, aunque menos inclinada, porque el factor de creci- 
miento es menor que 2. Una funcién lineal tiene dos propiedades 
fundamentales: 


1, St multiplicamos la variable por un valor determinado el resul- 
indo queda multiplicado por el mismo valor: 


Stk x)= k-f(x). 


#2. La Suma de dos valores da como resultado la suma de aplicar la 
funcién por separado a cada uno de ellos: 


S(%, + %,) = M(H) +f). 


In lo que nos ocupa en este libro, un sistema lineal nunca 
puede producir caos. En particular, la evolucién propuesta por 
Malthus solo da lugar a tres comportamientos perfectamente 
predecibles. Si el factor « es mayor que 1, la poblacion crece 
de manera exponencial, como en la figura 1, y a muy largo plazo 
se haria infinitamente grande. Si el factor @ es menor que 1, la 
poblacién disminuye en cada iteracién de manera continuada; a 
muy largo plazo la poblaci6n tiende a cero. Por tltimo, si el valor 
a es igual a 1, la poblacién no cambia. Un mundo lineal es un 
mundo perfectamente predecible, aburrido, cabria decir. Para 
que pasen cosas mas interesantes, el comportamiento caético 
entre ellas, las leyes que rigen los procesos naturales deben ser 
no lineales. 
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if 


0,4 0,6 


Representacién grafica de f(x) = 2xy la sucesién de valores a que da lugar su aplicacién reiterada: 2, 4, 8,... 


ROBERT MAY ESTUDIA LA DINAMICA DE POBLACIONES 


En 1976 un cientifico multidisciplinar, Robert May, estudié con 
detalle una modificacién de la ley de Malthus que representaba 
un modelo mas realista de la dinamica de poblaciones. May ha- 
bia estudiado Fisica en su Australia natal. Tras doctorarse en 
['{sica Tedrica pasé un tiempo en Harvard como profesor de Ma- 
temAética Aplicada, y en 1971 pas6 un afio en el Instituto de Es- 
tudios Avanzados de Princeton. Desde su fundacion, el objetivo 
del instituto ha sido el de proporcionar a los cientificos un am- 
biente relajado donde trabajar sin obligaciones, sin la presién de 
obtener resultados o completar proyectos en un plazo dado, una 
especie de «monasterio cientifico». En Princeton los cientificos 
pueden embarcarse en proyectos que no Ileven a ninguna parte, 
sin temor a ser expulsados. Alli May entablé relacién con los 
bidlogos de la universidad y empezo a interesarse por los proble- 
mas de la ecologia teérica. May ofrece un claro ejemplo de hasta 
qué punto el caos determinista es un campo multidisciplinar. El 
sistema dindénico que él estudié ya hab{fa sido considerado por 
otros bidlogos como modelo simplificado para la dinamica de 
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poblaciones, Pero, donde otras vieron lrrewularidades y eoripor 
Fimenitos eifioe de ciniowe interés, é)vio olge que mereeta da 
pena caning en Cetalle, 

Come el propio Malthus habin eomprendids, el crecimiento 
exponencial no podia mantenerse indelinidamente, La falta de 
recursos limitaria necesariamenté él crecimiento. Para tomar en 
eventa este hecho, el modelo estudiado por May considera- 
ba que el coeliciente @ debia depender de la poblacién, x, de 
modo que, a mayor poblacién, menor tasa de crecimiento. May 
adoplo una forma en la que o, la tasa de crecimiento, decrece 
eumndo crece la poblacién, 4%, de una manera sencilla: 


Is decir, si la poblacién, x, es una fraccién muy pequefia del 
valor Constante «, & es esencialmente constante y su valor es 
muy proximo a7, que representa la tasa de crecimiento mas alta 
posible. Por ejemplo, si x es una centésima parte de x,,, el mame- 
ro entre paréntesis vale 0,99 y 0=0,997. Pero a medida que cre- 
ce la poblacion, la tasa de crecimiento « va disminuyendo. Si la 
poblacion es la mitad del valor x, &S decir, =0,5x,, la tasa de 
crecimiento a disminuye a la mitad: r/2. Por ultimo, la tasa 
(le crecimiento se hace cero para x=. La constante %,, Yepre- 
senta, pues, el maximo valor de poblacién que el entorno puede 
admitir, ya que alcanzado este valor la poblacién se extinguiria 
en el siguiente paso. 

Para que la ecuacion tenga un aspecto mds manejable, vamos 
a redefinir la variable x de forma que represente no ya el nime- 
ro de individuos, sino la fraccién entre el niimero de individuos 
existente y el total admisible, x,,, De esta forma, en adelante, 
“= 1 representa una poblaci6én con el nimero maximo de indivi- 
duos, %,,; £=0,5 representa una poblacién con la mitad del valor 
x, etcétera. Con este criterio, la ley de May se escribe asi: 


©, =rd-2%)a,. 


n+ 
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Thy que entender esta ley coma un modelo no solo de Ia po 
Hlaeion himane mundial, aloe de cualquier eapecie en un ecosis 
tema dado, Las diversas eractones corresponden a la poblacion 
Sada cierto intervalo [jo de liempo, ya sean semanas, meses o 
floes, pero también pueden referirse a sucesivas eneraciones de 
individuos. Si la poblacién es muy pequeria, se supone que hay 
aificiente alimmento en el entorno y no hay depredadores especifi- 
os que se hayan desarrollado. La poblacién crece entonces con 
una tasa constante, como predijo Malthus. Pero, si la poblacién 
erece en exceso, empiezan a escasear los alimentos y proliferan 
Ins depredadores especificos, reduciendo el crecimiento. En caso 
ile superpoblacion, representada en la formula anterior por el va- 
lor w= 1, la poblacién alcanza el valor maximo posible, se produ- 
ce la catastrofe, y la especie desaparece en la siguiente iteracion. 

La funcién f(v)=r(1—x)a, que define el sistema dinamico 
asociado, es ahora una funcidén no lineal. Esta funcién tiene dos 
(érminos: f(a) =rz—ra*. El primer término es como el de la ley 
de Malthus y es lineal; el segundo es el que introduce la no linea- 
lidad, ya que el cuadrado de la suma de dos nimeros no es la 
suma de sus cuadrados y, por tanto, f(%, +x) #f(@,) +f(@,). 

Esta funci6n esta representada en la figura 3. La aplicaci6n 
recibe el nombre de aplicacién logtstica. 


Representacion grafica 
ly de f(x} = r(1—x)x 
1 para r= 2. j 
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May doseubrid que este sence modelo da lugar & eomporta, 
Mmiontos muy varlados dependiondo de in law de ereeinients 7, 
Vamos a explora, 

Siel valor de? ee inferior a 1, la poblaeion dlaaninuye en cada 
ifeyaeion y, al cabe del Gempo, la eepecte ae extingue, Este Com. 
portamento 6s ¢! que vemos en la iigura 4, en la que hemos re- 
presentado 1a evalucién de la poblacion para r= 0,5. lin este caso 
ln poblacion se redues mas oO menos ala mitad en cada paso. En 
(lev ilernctones casi ha desaparecido. 

Sires mayor gue L, la poblacién comienza creciendo. Pero 
él derming —e que hemos afadido en el paréntesis evita que la 
poblacion erezea indefinidamente, Como pasaba. en el modelo de 
Malthus, Para valores de r entre 1 y 3, la poblaci6n acaba por 
ostabilizarse en un valor fijo. Esto lo vemos en la figura 5, don- 
dle hemos representado la evolucién de la poblacién para r= 1,5. 
tion este parametro la peblacién se estabiliza para un valor de 
w proximo a 0,33 independientemente del valor inicial de x que 
hayamos utilizado. 

Pero ahi no se agotan las posibilidades. Si seguimos aumen- 
{ando 7, aparecen comportamientos mas variados. Por ejemplo, 
pura 7 = 3,2 la poblacién acaba por oscilar entre dos valores de- 
linidos, [in Ja figura 6 hemos calculado la evolucion del sistema 
para este valor empezando con x=0,1. Inicialmente la pobla- 
Mon crece y después de un par de iteraciones empieza a osci- 
lar entre dos valores que al principio estan bastante cercanos y 
luego se separan. Al final se alcanza un régimen de oscilacion 
entre dos valores definidos, repitiéndose el mismo valor cada 
dos iteraciones. 

Incrementando un poco mas r se llega a una situacion en la 
que la poblacién pasa por un ciclo de cuatro valores diferentes, que 
van repitiéndose sucesivamente, tal y como vemos en la grafica 
de Ia figura 7. 

Las sucesivas transiciones a ciclos de repeticién cada vez mas 
altos se producen cada vez mas pronto. Y a partir de 7=3,57 la 
dinamica se vuelve cadtica. 

El comportamiento observado parece indistinguible del con- 
junto de datos que arrojaria un proceso completamente aleato- 
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0,26 | 
0,24 
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4 6 8 10 


Evolucién de la poblacién segun el modelo de 
May para r= 0,5. Al cabo de varias iteraciones 


la poblacion desaparece. se alcanza un valor estable. 


rio, como puede ser tirar un dado. Por ejemplo, los valores de 
la centésima iteracién y siguientes son: 2,,, = 0,766; %,,, = 0,659; 


W 99 = 9,827; X93 = 0,526; Xo, = 0,917... Sin embargo, estos datos 


no se corresponden con los que obtendriamos en un proceso 
puramente aleatorio, en el que los valores obtenidos no estan 
correlacionados unos con otros, no dependen unos de otros, 
sino que son totalmente independientes. Esto no ocurre con 
las series caéticas como la de Ja figura 8, porque los datos 
obedecen a una ley determinista y siempre hay una estructura 
subyacente. En este caso es facil comprobarlo (figura 9): basta 
representar x, ,, frente ax, para recuperar parte de la parabola 
de la figura 3. 

.Cémo podemos recopilar de forma compacta toda esta varie- 
dad de comportarmentos? La forma habitual de hacerlo es repre- 
sentar en el eje de abscisas los valores de 7, y eneleje de ordenadas 
los valores de x, que se obtienen para un numero suficientemente 
alto de n, de forma que se haya eliminado la evoluci6n transitoria 
de las primeras iteraciones. Al hacerlo se obtiene el diagrama de 
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May para r= 1,5. Al cabo de varias iteraciones 
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-Evolucién de la 

poblacién segiin 

el modelo de May 

= para r= 3,678. 

ee eeeeoe0eeseecee8ee ee e La poblaci6n varia 

ee ° erraticamente. 
ee 

e 


Representacién 
dex, , frente 
a X, para los valores 
del modelo ce 

May obtenidos 


con r= 3,678. 


ZQUE ES EL CA0S? LQUE ES EL CAOS? 


No solo an investigacion, alo 
on ol dia a dia, an politica y 
goonamla, nog irla mejor a todos 
al mas gente entendiera que 
Alstemas no lineales simples 

no necesariamente presentan 
propladades dindmicas simples. 


In fiona 10. Mere dingenma, que May 
ya viatumbyro, se eaGnoee como tite 
game de Feigenbaum, en honor de 
Mitchell Peigenbaum, matematico y 
fisico esladounidense del que habla- 
remos mas adelante. 

La interpretacién de este diagra- 
ma es la siguiente. Para valores de 7 


ie, 


entre 0 y 1, todas las posibles evolu- 
ciones de la poblacién terminan en 
la extinecién: «=0. Este valor repre- 
senta un punto de equilibrio del sistema dinaémico y en este 
raingo de valores es un punto estable del sistema. Entre r igual 
alyrigual a3, la poblacién tiende a estabilizarse en torno a un 
valor distinto de cero. Este valor crece a medida que crece 7’. Si 
se aumenta ry por encima de 3, la poblacién ya no se estabiliza 
en un nivel concreto, sino que oscila entre dos valores. En r 
igual a 3,449 se produce un nuevo cambio y la poblacion em- 
pleza a oscilar entre cuatro valores diferentes. Los puntos en 
los que se produce el cambio de comportamiento, esto es, r= 1, 
re}, r=3,449, etcétera, se conocen como puntos de bifurca- 
cron. Al aumentar 7, las bifurcaciones se producen cada vez 
mas pronto y el sistema va pasando por ciclos de 2, 4, 8, 16,... 
valores, El] namero de puntos en el ciclo tiende a infinito para 
un valor concreto de 7, cercano a 3,57. A partir de ese punto la 
dinamica se vuelve cadtica. 

lista transicién al caos a través de ciclos cuyos periodos se 
van doblando (2, 4, 8,...) se conoce como transicién al caos por 
desdoblamiento de periodo o cascada subarmoénica. La aplica- 
cién logistica es quiza el ejemplo matematico mas simple que 
presenta este comportamiento. Pero este modo de transici6n al 
caos se ha encontrado, desde los estudios de May, en sistemas 
muy diversos, ya sean naturales o artificiales. En particular, se 
observo experimentalmente por primera vez en el movimiento 
de un fluido calentado desde abajo, como veremos mas adelante. 

En lo que se refiere a la dinamica de poblaciones, el modelo 
es muy simplificado, pero deja un mensaje radical: un comporta- 


Ropernt May 
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Diagrama de Feigenbaum. A partir de un valor de rcercano a 3,57 la dinamica se vuelve 
cadtica. 


miento complejo no se debe necesariamente a una multiplicidad 
de causas. Por ejemplo, para r=4 la poblacién varia de iteracién 
en iteraciOn barriendo todos los valores posibles, de 0 a 1. Puede 
darse el caso de que una iteracién produzca un valor tan cercano 
a0 que la poblacion se extinga sin mas. No hacen falta catastrofes, 
epidemias, invasiones externas o meteoritos: la propia dinamica 
interna puede llevar a la extincidn. 


Autosimilaridad 
Ya hemos explicado que los sistemas cadticos deterministas 
se caracterizan por la sensibilidad a las condiciones iniciales, 


por la irregularidad en el comportamiento y por la no lineali- 
dad de las ecuaciones que los modelan. El diagrama de Feigen- 
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ee ——————————— 


Denn now Huet ane propiedad mviy de este tips de aston: a 
itoalinilenidaed, 

La tigure 1) representa toe mplivelon de la sone alta det dis 
grwnia de da Sgure 10, Le linea inlelal ne correspende a un punto 
estable de la poblacion, sino que se corresponde con los valores 
altos del cielo que se produce tras Ja segunda bifureacion. Ve- 
THOS que el aspecto de este dlagrama es similar al del diagrama 
anterion Si, a su vez, hiciéramos un geom de este diagrama en 
torneo & re 3,56 obtendriames otra vez una figura similar. Solo 
ceambiarian en la grafica kes escalas de ry de a, cuyo rango de 
verlacion seria cade vez mas estrecho. Esta caracteristica de una 
figura o diagrama, que presenta la misma estructura a escalas 
onda vez mds pequefias, es lo que se conoce como autosimilari- 
dad, y esta presente de una u otra forma en todos los sistemas 
CAOLICOS, 

May escribie un articulo para la revista cientifica Nature en 
1976 en el que resumia todos sus hallazgos, los mas importantes 
de los cuales lemes descrito aqui. El articulo terminaba con una 
linimada a los responsables del sistema educativo para que toma- 


ae ee eee. Fg 


"34 345 35 3,55 36 369 37 


feem el dlagrama de Feigenbaum para valores de r comprendidos entre 3,35 y 3,7 y valores de x entre 


O75 909. 


(QUE ES EL CAOS? 


ran conclencgia de le iportaneia de kano Uinealidad en nuestra 
entorne, May proponia que el estudio de la aplicaeidn logiaien, 
que se puede hacer con una calenladora de bolsillo, se incluyera 
on la asignatura de Matematicas a una edad temprana, ya que no 
oe mais dificil que las derivadas o las integrales. Pasados casi 40 
ahos, la propuesta de May no parece haber encontrado aun el 
eco que él hubiera deseado. 


EL CAOS EN UN POCO DE Liquipo 


Las primeras décadas del siglo xx son, en fisica, los afios de la fi- 
sica at6mica y nuclear. Los grandes tedéricos y los mejores fisicos 
experimentales pugnaban por desentrafiar las leyes del atomo y 
establecer la naturaleza de la radiactividad. Pero no todo ej mun- 
tlo se dedicaba a estas disciplinas. Una rama, que habia sido muy 
importante durante el desarrollo de la fisica y la tecnologia en el 
siglo xrx, siguid llamando la atencion de algunos de los mejores: 
la fisica de los fluidos. Uno de ellos fue Henri Bénard, un fisico 
francés cuyas habilidades experimentales le han hecho pasar a 
la historia por el descubrimiento de las celdas de conveccién en 
un liquido que se calienta desde abajo. 

Henri Bénard nacié en un pueblo de Normandia, en Francia, 
el 25 de octubre de 1874. En 1894 super6é el examen de ingre- 
so en la Escuela Normal Superior de Paris. La Escuela Normal 
era, y sigue siendo, uno de los centros de educaci6n superior 
mas prestigiosos de Francia, y fue fundada en la época de la 
Revolucion. Una idea del nivel de exigencia de los ex4menes 
de acecso y de la escuela la da el hecho de que Bénard fue 
uno de los 17 que superaron el examen de entre 307 candidatos 
llegados de toda Francia. Compafieros de promocion de Bénard 
en la Escuela Normal fueron dos de los mas grandes cientificos 
franceses del siglo xx: el fisico Paul Langevin y el matematico 
Henri Lebesgue. 

Al iniciar su tesis doctoral Bénard se interesé por un tema de 
moda: la recepcién de ondas electromagnéticas. Investigando 
nuevos materiales para la deteccién de las ondas observ6, por 
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Cosieidad, In foriaclin de Duis polyonwes an in baa de 
porating Nendida nla que se le habia afacide polve de pratita, 
Normaimente la densidad de un fuido diaminuye con la lempera: 
lair) cuanto mas caliente, menos denso, Eso hace que él aire ca- 
hente tienda a ascender, y lo mismo ocurre deniro de un liqui- 
do; i Ja temperatura no es uniforme las zonas mas calientes, y 
por tanto menos densas, tienden a ascender, poniendo al liquido 
én movimiento, il liquido caliente transporta Consigo el calor, 
iransferencia que recibe el nombre de conveccidn. En el expe- 
rimento de Bénard las particulas de grafito habfan servido de 
irazadores, revelando las estructuras convectivas que la parafina 
lfquida formaba al moverse impulsada por el calor aportado para 
fundirla, Las particulas de grafito se movian con el liquido y su 
movimiento reflejaba el de aquel, de modo parecido al humo del 
labaco cuando nos ayuda a ver los remolinos que se forman en 
¢l aire. Benard comprendié que las estructuras formadas mere- 
cian ser objeto de investigacién por si mismas e inicié su estudio 
sistematico, 

{in sus experimentos Bénard depositaba una fina capa de Ii- 
(juido sobre una placa metalica que calentaba por medio de va- 
por de agua. El éxito de las observaciones dependia de muchos 
factores, siendo los mas importantes la constancia y la uniformi- 
dad de la temperatura en toda la placa. La mas leve variaci6n de 
{emperatura o el mas minimo defecto de la placa arruinaban el 
experimento. El uso de vapor de agua como elemento calefactor 
llevé a Bénard a usar sustancias que son sélidas a temperatura 
ambiente y funden en torno a 50 grados. Entre ellas destacaba 
un aceite conocido como espermaceti, que almacenan en la ca- 
beza algunos cetaceos. En esa época, en la que la industria ba- 
llenera estaba en su apogeo, el espermaceti se usaba para hacer 
velas y también como ungtiento medicinal, por lo que era una 
sustancia relativamente facil de conseguir. La figura 12 mues- 
tra lo que hoy en dia se conocen como celdas de Bénard, uno de 
los patrones de conveccién que Bénard observ6 en wna delgada 
capa de espermaceti. Al calentar el liquido por debajo, las zonas 
mas cercanas a la placa estan mas calientes que las capas mas 
superficiales. El liquido caliente tiende a ascender y el mas frio, 
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La imagen de la izquierda muestra unas celdas de Bénard en espermaceti vistas desde arriba. El esquema de 
la derecha ilustra el movimiento de conveccion dentro de una de las celdas, girada y vista en seccion vertical. 


El I(quido asciende por el centro y baja por los bordes de cada celda. 


a descender. Si el calentamiento es muy uniforme y la capa de 
liquido muy delgada, el flujo se produce de forma ordenada. En 
el patr6én hexagonal de la figura, Bénard pudo determinar que el 
liquido ascendia por el centro de los hexagonos y descendia por 
los bordes. 

Que un sistema inanimado produjera semejante patrén hexa- 
gonal de manera espontanea podia parecer, a primera vista, sor- 
prendente. Sin embargo, hay que tener en cuenta que sea cual 
sea el patron debe llenar toda la superficie del liquido, y si no hay 
zonas diferenciadas, esto es, si la superficie es totalmente lisa, 
de espesor perfectamente constante y tiene temperatura unifor- 
me, el problema se reduce al de rellenar un plano con un patron 
regular. Los matematicos saben desde muy antiguo, ya lo sabfan 
los artesanos andalusies que decoraron las paredes de la Alham- 
bra de Granada, que solo hay tres formas de llenar una superficie 
con un motivo poligonal regular: usando cuadrados, triangulos y 
hexdégonos. Asi que el patrén observado por Bénard era una de 
las escasas posibilidades. 
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Tlenrl Deétard Tue un extriordindrio experimentadar Sus ca 
pan de espermaoel tenian apenne un Milimetrs de prosor y lh 
placa sobre la que reposaban debia estar en poxicion perfecta- 
mente horizontal para que el capesor no fuera mayor en unas 
7On08 que en otras. Para hacer visiblos las celdas Benard uso la 
reflexion producida en la superficie, deformada como un espejo 
concave, en él centro de las celdas. Fue capaz, de esta forma, 
de mediy diferencias de espesor del orden de una micra en ca- 
pas de un milfmetro que estaban a unos L100 °C de temperatura. 
Usando una técnica éptica que aprovechaba la luz reflejada en la 
superficie del liquido y en la pulida superficie metalica del fondo, 
consiguid medir diferencias de temperatura del orden de 0,1 °C 
en él interior de las celdas de conveccidén. 

Si el siglo xix vio nacer la fotografia, el cambio de siglo trajo el 
cinematografo. En la cultura moderna tanto la fotografia como 
el cine suelen estar asociados al arte y al espectaculo, pero esta 
concepcién ignora la enorme importancia que el desarrollo de 
ambas técnicas ha tenido en la ciencia. Henri Bénard fue uno 
dé los primeros cientfficos en utilizar el cinematdégrafo como téc- 
nica experimental. Grab6é sistematicamente sus experimentos y 
ello le permitié observar el movimiento del liquido dentro de las 
ealructuras convectivas. 

Después de pasar unos afios en Lyon y Burdeos, en 1922 Bé- 
nard fue nombrado profesor en la Facultad de Ciencias de la 
Sorbona, en Paris, ciudad en la que residié hasta su muerte en 
1939, Dirigid numerosas tesis doctorales y form6 a toda una ge- 
neracion de cientificos franceses que dedicaron su carrera a la 
mecanica de fluidos, todos ellos distinguidos por una sélida for- 
macion experimental. 

Los experimentos de Bénard llamaron la atencidn del fisico 
inglés John William Strutt, lord Rayleigh, que en 1916 propuso 
una explicacion tedrica de las observaciones del francés. Des- 
de entonces, el problema de un liquido que se calienta desde 
abajo es conocido como problema de Rayleigh-Bénard. Si re- 
presentamos la velocidad maxima de movimiento del lfquido 
en funcién de la diferencia de temperatura entre el fondo y la 
superficie, obtenemos un diagrama parecido al de la figura 13. 
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Valookiad dé éGonvacclan on 
funclén de ja difaroncla de 
tomporatura en el problema 
de Rayleigh-Bénard. 


Por debajo de la diferencia de temperatura critica AT, el If 
quido esta en reposo y la velocidad es cero. A partir de AT, el 
lfquido se pone en movimiento con una velocidad cada vez ma- 
yor. El valor de AT, representa un punto de bifurcacién, y el dia- 
grama nos recuerda inmediatamente a la primera parte del dia- 
grama de Feigenbaum de la figura 10, a la primera bifurcacion, 
en torno a r= 1, cuando la poblacién pasa de ser cero a crecer 
poco a poco. Qué cabe esperar si aumentamos la temperatura 
un poco mas? 


Helio en una caja pequefia 


En 1977 otro gran experimentador, también francés y también 
normatien, esto es, exalumno de la Escuela Normal, llev6 los 
experimentos sobre la conveccién de Rayleigh-Bénard a un gra- 
do de precisién y sofisticaci6n desconocidos hasta entonces. Si 
Bénard habia usado vapor de agua para controlar la tempera- 
tura, y eso lo hab{a llevado a usar sustancias que son solidas a 
temperatura ambiente, Albert Libchaber buscé su liquido entre 
las sustancias gaseosas a temperatura ambiente: el helio. 
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Libehaber ncié on Marie en LM, en ef geno de tina familia 
Judin de origen polaco, Cuando los navies invadieron Mraneia sus 
padres lo enviaron, Junto con su hermano mayor, Marcel, al sur 
Ue Mranela, con la esperanza de que alll tuvieran mas posibili- 
Jades de sobrevivir, Se hicieron pasar por huérfanos catélicos 
provententes de Alsacia y consiguicron la ayuda de algunas per- 
sonas, enire ellas un jefe de la policia local de Petain. Aunque 
aus padres sobrevivieron a la guerra, su madre perdié a toda su 
familia, y su padre ala mayor parte de la suya. 

Tras graduarse en Matematicas en la Universidad de Paris en 
1966 y como ingeniero de telecomunicaciones en la Escuela Na- 
cional Superior de Telecomunicaciones en 1958, Libchaber obtu- 
vo una beca Fullbright, que le llev6 a trabajar en Illinois con el ya 
premio Nobel de Fisica John Bardeen, una de las grandes figuras 
de la fisica del siglo xx y uno de los pocos cientificos que ha ga- 
nado dos veces el premio Nobel. De vuelta en Paris fue llamado 
a filas y enviado al Sahara para luchar en la guerra contra los 
independentistas argelinos. Enrolado en el escuadrén de armas 
alomicas, participO6 en uno de los primeros ensayos nucleares 
franceses, realizados en el desierto en torno a 1960. En 1965, 
Libehaber termin6 su doctorado en la Escuela Normal Superior. 

La forma de enfocar la ciencia por parte de Libchaber siempre 
fue muy original. Segun él mismo decia en una entrevista: «Soy 
judio y francés. Francés es ser racional, matematico. Judio, una 
vision mistica del estudio y el aprendizaje». Para él los cientffi- 
cos y los conocedores del Talmud hacen algo similar: «Hay un 
mensaje en clave, que hay que descifrar». 

Desde su estancia en Illinois, Libchaber se interes6 por la fisi- 
ca de las bajas temperaturas, lo que le llevé al estudio del helio 
superfluido. El helio liquido, estado en el que se encuentra este 
elemento por debajo de 4,22 K (4,22 grados por encima del cero 
absoluto), presenta propiedades sorprendentes si desciende la 
temperatura por debajo de 2,17 K. En particular se anula su vis- 
cosidad, por lo que fluye sin friccién. Ese estado se conoce como 
superfluido. 

En torno a 1977 Libchaber decidi6é usar sus conocimientos de 
criogenia para estudiar el viejo problema de Rayleigh-Bénard en 
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Gétta diferencia de jarfocreting critica por debajo de , 


he pudo leer el artioule de Rayleigh hasta ¢ ous acabé ta guerra, cuando 5s Rvisigh y ya 
muerto, por lo gue no pudo discutir'con él ios resultados. 5] Caso es que Bénard tomo” | 
Josas medidas tanto de! tamaria de sus celdas de conveccién como de la diferencia . 
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hetie Tiqiite on estade nord, co eupertiide, a tres yradoa por . ’ 


enol del cero absohito. {Que veritatis podia toner este iqnide 
frente 2 ioe uandos por Henared? Aun sin ser superttiuide, el helio 
liquide presenla una viscosidad muy inferior ala de la mayoria 
de los liquidos, y habia sole unos experimentos preliminares, 
ronlizados por Guenter Ahlers, de la conveccion de Rayleigh-Beé- 
hard en estas condiciones. 

Libchaber uso wna caja de acero noxidable que, en su inte- 
nor, tenia del orden de un milimetro de ancho y de profundo y 
menos de un milimetro de alto. En una caja tan pequefia solo era 
de esperar que cupiera una celda de conveccién, con el liquido 
ascendiendo por el centro y bajando por las paredes, formando 
dos «rollos» paralelos. La pequefia caja de acero tenia una entra- 
da aun mas pequena para el helio liquido y estaba encerrada en 
una camara de vacio sumergida, asu vez, en un bano criogénico. 
La parte inferior estaba en contacto con una placa de cobre cuya 
fomperatiura era controlada por una delicada resistencia eléctri- 
cu, La parte superior estaba hecha de un cristal de zafiro. Una 
sonda en miniatura monitorizaba la temperatura en un punto de 
a parte superior de la celda. La temperatura de esta sonda local 
s0ivia Como marcador del movimiento del liquido. Un rotulador 
registraba sobre papel continuo como evolucionaba con el tiem- 
po la temperatura de la sonda. Para que nos hagamos una idea 
de lo delicado del experimento montado por Libchaber, baste 
decir que la conveccién del helio dentro de la caja se iniciaba 
cuando Ja diferencia de temperatura entre la parte inferior y la 
superior era de unas dos milésimas de Kelvin, y la temperatura 
se controlaba con una precisién de hasta un microkelvin (una 
millonésima de kelvin). 

Tras casi dos afios de minucioso trabajo liegé el momento de 
poner en marcha el experimento. Todo el equipo trabajaba dia 
y noche y la temperatura de la sonda encargada de monitorizar 
la convecci6n arrojaba sus datos sin parar. Al principio, una vez 
establecida la conveccién, la sonda marcaba una temperatura 
constante, lo que queria decir que el helio se movia dentro de su 
celda a una velocidad que no cambiaba con el tiempo y de acuer- 
do con un patr6én establecido. Al ir aumentando la diferencia de 


{QUE ES EL CAOS? 


lenperaturn, Libehaber eorproba que riuy pronto, justo por 
encima del umbral en él que el helo empezaba fo moverse, el 
rohilador que mareaba la evolucion de la temperatura se des- 
plazaba de un lado a otro periddicamente, oscilando cada dos 
aepundos. La interprelacion que cabia hacer de este hecho era 
que la columna central ascendente de liquido se bamboleaba a 
un lado y aotro de la celda. Para una diferencia de temperaturas 
todavia mayor aparecia otro periodo, ahora de cuatro segundos, 
superpuesto al anterior. La temperatura de la sonda presentaba 
maximos cada dos segundos, como antes, pero estos maAximos 
no eran iguales, sino que uno era mas alto que el otro, repitién- 
lose la secuencia entera cada cuatro segundos. 

Libchaber encontr6 asf una secuencia continua de bifurca- 
ciones cada vez mas cercanas unas de otras. Cada bifurcacion 
aportaba una periodicidad cuyo periodo duplicaba al anterior. 
A partir de un cierto valor de la diferencia de temperaturas, el 
comportamiento de la sonda y, por tanto, el del liquido, se volvia 
caotico. Albert Libchaber no sabia entonces que ese compor 
tamiento que é] hab{a encontrado en el helio liquido encerrado 
en una pequefia caja era el mismo que Robert May habia visto en 
su aplicacién logistica. Libchaber habia descubierto el primer 
sistema experimental que mostraba un comportamiento cadtico 
determinista. 


UNIVERSALIDAD 


Mitchell Feigenbaum nacié en Nueva York en 1944, hijo de inmi- 
grantes judios. En 1970 se doctoré en Fisica de Particulas en el 
Instituto Tecnologico de Massachusetts, el MIT. Tras pasar por 
otro de los grandes centros universitarios de Estados Unidos, la 
Universidad de Cornell, recalé en el Laboratorio Nacional de Los 
Alamos. Su objetivo alli era, en principio, estudiar la turbulencia 
en los fluidos. En particular, comprender cémo se producia la 
transicién de un movimiento laminar a uno turbulento. 

En un flujo laminar el fluido se mueve de manera ordenada 
y suave. Cada particula sigue una linea bien definida, llamada 
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Hi milamo Feigenbaum atribuye au descnubrimienta al heehe de 
que use ina ealeulndora de batalla para obtener loa muicealvos 
valores de 7. Otros investigadores que, como May, se habian inte- 
pesado por la aplicacion log{stica, usaban ordenadores para rea- 
lize las iteraciones de manera auto- 


lined de corriente, ¥ las lineas de corrionte diseiirren pirilelan 
ies © OLS, aii mezelinde, Come Formanda Wiis, y sin eam 
blar con ¢] tlempo. Por el contrario, ¢l movimiento furbulento 
se cormeteriza por la irregularidad, Las particulas del Muido si, 
fuen trayectorias cambiantes y su velocidad fluctua continua- ; 
mente, Ademas, en un flujo turbulento abundan los torbellinos. matica. El hecho de tener que buscar De alguna forma lo maravilloso 

A velocidades suficientemente bajas los flujos suelen ser lami- los puntos de bifureacién uno por de la Tierra es que hay cosas 
nares, como el humo justo al salir de la punta de un cigarrillo. uno, to que llevaba su tiempo, obliga- bellas en ella, cosas maravillosas 
Pero, al aumentar la velocidad, los flujos se vuelven turbulentos, rie a Feigenbaum ° tratar de adivinar y atractivas, y gracias a tu trabajo 
“omo también vemos en el humo de un cigarrillo un momento donde se produciria la proxima bifur- ; 
después. Esta transicién de flujo laminar a turbulento era lo que cacién y eso le llevé a descubrir que uleres llegar a comprenderla. 
Mitcuent Feigensaum 
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(nteresaba a Feigenbaum en Los Alamos. 

Hemos visto que el comportamiento de la poblacién en el mapa 
logiatica usado por May pasaba por sucesivas bifurcaciones en 
las que el periodo se iba duplicando, hasta que, finalmente, la 
evolucion del numero de individuos de la poblacién se volvia 
castica. Este mismo comportamiento de desdoblamientos de pe- 
riodos sucesivos fue observado experimentalmente por Libcha- 
ber en su caja de helio. Estando en Los Alamos, y con el fin de 
reflexionar sobre la forma en que la turbulencia aparecia en los 
fluidos, Feigenbaum decidié tomar la aplicacién logistica como 
modelo, Jugando con una calculadora de bolsillo, descubrié que 
las sucesivas bifurcaciones obedecian a una ley sencilla. Para 
verlo, designamos con 7, los valores en los que se produce la 
bifurcacién con periodo 2”, y definimos el siguiente cociente: 


T44—7, 
sks n+l n 

6, = : 
Tria ~ Vast 


Al aumentar el valor de n, los valores de 6, tienden a un nu- 
mero fijo 6=4,6692016... Es decir, las distancias sucesivas entre 
dos bifurcaciones van disminuyendo en progresién casi geomé- 
trica, por lo que, al final, los valores de 7, convergen a un valor 
fijo 7, =3,5699456... Este valor de 7 corresponde a un periodo de 
repeticion infinito, es decir, la secuencia de valores no se repite 
nunca. Representa el limite a partir del cual se establece el caos. 
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se sucedian en una progresién casi 
geométrica. 

Habiendo encontrado una regla 
que determinaba los puntos de bifurcacién de la aplicacién 1lo- 
pisitica, la usada por May, Feigenbaum decidié hacer lo mismo 
para otra funcidén sencilla, Probé con otra funcién que no conte- 
iia un término cuadratico como la de May, sino una funcion tri- 
gonométrica. El nuevo sistema presentaba la misma cascada de 
bifurcaciones. Usando el mismo cociente para calcular los suce- 
sivos valores de ce se llevé una sorpresa mayiiscula: convergian 
al mismo valor § = 4,6692016... 

Feigenbaum pronto comprendi6 que el valor del ntimero 6 no 
dependia de la forma particular de la funcién que escogiera. Bas- 
taba con que la funci6n tuviera un maximo similar al que obser- 
vamos en la representaci6n de la figura 3 (pag. 53), para que las 
distancias entre las sucesivas bifurcaciones convergieran a este 
numero. 

Como hemos dicho, Feigenbaum habia estudiado fisica de 
particulas en el MIT. Estaba familiarizado con las técnicas de re- 
normalizacién que se aplican en electrodinamica cuantica para 
eludir ciertas cantidades infinitas que aparecen en la teoria. Fei- 
genbaum uso técnicas similares para demosirar la existencia de 
6 y calcular su valor. La idea fundamental del procedimiento que 
uso reside en la autosimilaridad que vimos reflejada en los dia- 
gramas de las figuras 10 y 11 (pags. 59 y 60): si hacemos un zoom 
en el diagrama de la figura 10, obtenemos una figura similar a la 
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original, Se trata, pues, de eneontivar low valares de eweala ade- 
CHAOS (0 Oe TeV de init [pe a oben, 

Loa valores (jos del diagrama de Veigenbaumn corresponden a 
eoluciones de in eounclon we f(a) Por ejemplo, & 762 le corres: 
ponde en el diigirama de la figura 10 el valor a 0,5, que 68 solu 
elon de x ery (1<a) para ese valor de 7, como es facil comprobar. 
Cuanda se produce la bifureacion que Neva al primer desdobla- 
mento de periodo, el comportamiento que vemos en la figura 6, 
log puntos del ciclo de dos valores que se observa son solucién 
dé ve f(fCe)), porque tras dos aplicaciones de la funcién vol- 
vemos al mismo valor. Por ejemplo, para r=3,2 los dos puntos 
del diagrama de la figura 10 corresponden a a=0,518 y #=0,8, 
de forma que f(0,8)= 0,518 y £(0,513)=0,8. Luego se tiene que 


JUF00;8)) 0,8 vy SC/(0,513)) = 0,518. Tras la siguiente bifurcacién, 


loa valores de] ciclo de orden 4 son soluci6n de x=fCf f(a), 
y ast sucesivamente. Las funciones compuestas, como f(f(%)), 


JUS (PD) y siguientes, tienen una forma fastidiosa, como el 


lector puede imaginar, pero Feigenbaum estudi6é su forma apro- 
ximada, alrededor de los puntos solucién y encontré que todas se 
podian superponer si se escalaban con un factor adecuado. Ese 
factor permitia calcular el valor de 5. 

Al darse cuenta de que la constante 6 no dependia de la for- 
ma concreta de la funcién que utilizaba, Feigenbaum compren- 
diG que su descubrimiento era algo profundo. Suele ocurrir 
en fisica que los detalles cualitativos de muchos fenémenos 
son similares, pero los valores cuantitativos son diferentes. 
Por ejemplo, las transiciones de fase de sdlido a liquido o de 
liquide a vapor de las sustancias quimicas se producen siempre 
de manera similar y se describen con la misma teoria termodi- 
namica, pero la temperatura a la que se producen, los calores 
latentes de fusién o ebullicién, o la cantidad de calor que la 
sustancia intercambia con el exterior durante el proceso, son 
diferentes en cada caso. En la cascada de Feigenbaum ocurre 
justamente lo contrario: funciones diferentes, que describen 
fenédmenos muy diversos y cualitativamente diferentes, llevan 
a una cascada sucesiva de bifurcaciones que dan lugar a una 
misma cantidad, 5. 
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Cuando Alber Libehaber aupe de ln tear de Metenbaain cab 
iid el cocienie enive tos sucesivos valores de bifureacion de la 
diferencia de Gampernwnras entre las placas superior ¢ inferior de 
fu eaja de helio, Dentro de las ineertidumbres experimentales, 
los sucesivos coclentes se aproximaban al valor predicho por 
leigenbaum. 

é,Por qué un modelo tan simple como el de May y Feigenbaum 
és relevante para algo tan complejo como un fluido en movi- 
miento? Recordemos que el nimero de grados de libertad de un 
sistema es el nimero minimo de nimeros que es necesario dar 
para especificar su estado. Pues bien, en un medio fluido los gra- 
dos de libertad disponibles son, en principio, infinitos, ya que 
hay que especificar la velocidad y la temperatura de cada punto. 
Pero la caja de Libchaber es tan pequefia, y la diferencia de tem- 
peraturas tan exigua, que el helio se mueve dentro de forma bas- 
tante ordenada. El liquido esta forzado a moverse en dos rollos, 
luego las particulas del fluido se mueven de forma coordinada, 
y de los infinitos grados de libertad existentes en principio, solo 
unos pocos son relevantes. Una vez que el sistema se desestabi- 
liza, silo hace por una bifurcacion por desdoblamiento de perio- 
do, la aplicacién logistica capta lo esencial del proceso. 

Feigenbaum titul6 sus articulos usando expresiones como 
«comportamiento universal», «universalidad cuantitativa» 0 «pro- 
piedades métricas universales». {Qué queria expresar con el uso 
del adjetivo universal? Lo que Feigenbaum pretendia resaltar con 
el uso de esta palabra era que el campo de aplicacién de su teoria 
era enorme, y que el escalado de las sucesivas bifurcaciones seria 
siempre el mismo en todos los fenédmenos que presentaran una 
cascada subarmonica, siempre y cuando la funcién subyacente 
presentara un maximo simple. 

No se puede objetar nada a Feigenbaum por usar el térmi- 
no universal en este contexto. En fisica se utiliza la expresién 
constanite universal para referirse a las constantes fundamen- 
tales que aparecen en las leyes fisicas basicas. Por ejemplo, en 
la ley de la gravedad de Newton aparece una constante que se 
denomina constante de gravitacién universal. Esta constan- 
te es universal en el sentido de que es la misma para todos los 
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Inversion 


Escalado x 2,5 


tia calcular §. 
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cuerpo, Independientemente de si comporiclon o de enagquier 
otra variable, y os ti Toikma en todos los hiagares y en todos los 
Hiempos, En ese mismo sentido, los valores de eseala encontra- 
dos por Feigenbaum son también aniversales, Sin embargo, esta 
denominacion ha Nevado a algunos escritores seudocientificos y 
esoléricos a ver en el desecubrimiento de Feigenbaum un conte- 
nido mistico del que esta claramente exento. 

La cascadasubarmonica, 0 de Feigenbaum, o de desdoblamien- 
Lo de periodo, que todos esos nombres recibe, es solo una de las 
rulas conocidas hacia el caos. Fue la primera en ser descrita en 
detalle y, quiza, la mejor entendida de todas. Con el tiempo, fisi- 
cos y matematicos han encontrado y descrito otras rutas hacia el 
vHos que Comparten Con la de Feigenbaum varias caracteristicas, 
siendo Ja universalidad, en el sentido ya explicado, una de ellas. 

En el primer capitulo vimos dos de las sefias de identidad del 
eaos determinista: un comportamiento irregular y la sensibili- 
dad a las condiciones iniciales. Esta segunda propiedad signi- 
flea que, partiendo de dos estados muy préximos, un sistema 
cadlico puede evolucionar a estados futuros muy diferentes. 
En este capitulo hemos afiadido dos caracteristicas nuevas: la 
autosimilaridad y la universalidad. Por autosimilaridad enten- 
dlemos que los sistemas cadéticos producen estructuras, como el 
cliagrama de Feigenbaum, que se asemejan a si mismas a distin- 
las escalas de observacion. Por universalidad entendemos el 
hecho de que sistemas muy diversos llegan al comportamiento 
cadtico a través de las mismas rutas. 
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~Qué tiene de extrafio 
un atractor extrano? 


El] comportamiento cadtico en sistemas 
deterministas gira en torno a un concepto 
matematico descubierto en la segunda 
mitad del siglo xx: el atractor extrano, 

un objeto geométrico de caracteristicas 
muy especiales. 


Los objetos fractales estan intimamente ligados a las matemati- 
cas del caos. En este capitulo una maquina imaginaria nos va a 
ayudar a entender qué es un fractal. Mas adelante veremos qué 
relacion tienen los fractales con los sistemas dinamicos. 
Llamaremos a nuestra maquina «fotocopiadora reductora 
multiple». Dada una imagen cualquiera, el triangulo de la figura 
la, por ejemplo, la fotocopiadora reductora multiple la reduce 
primero a la mitad, de ahi el calificativo reductora, y luego la co- 
pia tres veces, y de ahi que la llamemos miitiple. Las tres copias 
se disponen formando un triangulo equilatero, como muestra la 
figura lb. Nuestra maquina genera asi una figura de salida, que 
consiste en tres copias reducidas de la figura de entrada original. 
Utilicemos ahora la fotocopiadora reductora multiple reitera- 
damente, metiendo en ella como figura de entrada la figura de 
salida obtenida en el paso anterior. Si introducimos la figura ob- 
tenida en la primera aplicacién, obtendremos la figura lc, y si 
introducimos la figura 1c obtendremos la figura 1d. Repetir el 
mismo procedimiento seis veces produce la figura 2. El triangulo 
de Sierpinski, llamado asi en honor del matematico que lo imagi- 
no, es el conjunto que se obtiene repitiendo este proceso un nt- 
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o dv veces, RE oe Es 
de Git. porque ello supondrla tener unk precision infinita 
enol dibujo, lo Hgura 2 es ye une buena apreximacion a él, Como 


Landis veees ocurte en malemations, e] triingulo de Sterpinski es 


be Abetraceion, un objeto ideal, resultade de tin proceso infini- 
to, que en la praction solo es realizable de manera aproximada. 
COAntas fvis iteractones hagamos, mas nos acercamos al objeto 
matemAalico abstracto. 

MI) trifngulo de Sierpinski es un conjure fractal y tiene lodas 
ins caracteristicns de este tipo de objetos. Entre ellas destacan 
dow: la mutosimilaridad y la dimensién fractal, La autosimilaridad 
significa que al ampliar un fragmento de la imagen se nos revela 
una imagen igual a Ja original. Este concepto ya lo encontramos 
en el ceaptitule anterior, al referirnos al diagrama de Feigenbaum. 
Debido a la autosimilaridad, un fractal es una figura geométrica 
que io se revela mas simple cuando la vemos en escalas mas y 
ils pequefias cada vez, Por esta razon, un fractal es necesaria- 
mente un abjeto de aspecte quebrado, rasgo del que deriva su 
nombre, 

La segunda caracteristica de los fractales que queremos desta- 
car aqui es su dimension. El triadngulo de Sierpinski no es asimi- 
lable # na curva convencional y tampoco es una superficie. Por 
etio, su dimensi6n no es uno, como corresponde a una curva, 
iii dos, Gorno corresponderia a una superficie. En general, a un 
objeto. fractal no se le pueden atribuir las medidas habituales de 
jn geometria, como longitudes, areas o volimenes. En particular, 
el tridngulo de Sierpinski parece tener una longitud infinita, pero 
su Area es cero. Asi, se dice que la dimension de los fractales no 
es un namero entero, sino que es fraccionaria. Dimension frac- 
tales la expresion usada. La dimension atribuida al triangulo de 
Sierpinski es de 1,5849. 

La forma en que hemos construido el tridngulo de Sierpins- 
ki nos muestra dos caracteristicas adicionales de los objetos 
fractales que estan en el origen de su relacién con el caos: los 
conjuntos invariantes y el concepto de informacidén. La «fotoco- 
piadora reductora multiple» es una operacién que transforma un 
trozo del plano, el folio que fotocepiamos, en e] mismo trozo del 
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Construcelan de un triangulo de Sierpinski con una «fotocopiadora reductora miltiple». 


Wy, Ve" “a 


Tridngulo de Sierpinski (aproximacién obtenida tras seis aplicaciones sucesivas de la «fotocopiadora reductora 
multiple»). 
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s, como se muestra en la figura, pero sin reducirlo, La 
i =21, y un rea total a,=3a,, Si volvemos a-copiar - 
1a forma, el lado de la nueva figura mide /,=4/, y._ 


la longitud del lado del objeto y muttiplica por 3 su 
tendremos que la longitud del lado es /,=2"/, 
nia dlefiniclén usual, la dimension del triangulo de 


“ [log a, 93 : 
rar ees 
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Sierpinski 
1? en el limite nee. Tomando fogaritmos a uno y otro lado 


plano, otro folio de igual tamaie, Bn elerto sertide es un slater 
dindmico: par cade imagen inleia, la apiicacion reiterada de la 
Tolocopiadora genera una secuencia de imagenes, Esa secuencia 
va convergiendo a una imagen final, el triéngulo de Sierpinski, a 
1 que se acerca cada. vez mas. De hecho tal vez el lector se haya 
a cuenta de que Ia imagen de partida es irrelevante: da igual 

empecemos las iteraciones con un circulo, un cuadrado, un 
(riangulo, 0 una foto de la reina de Inglaterra, al cabo de muchas 
iteraciones ese detalle inicial se va difuminando y la figura del 
tridngulo va emergiendo poco a poco. El tridangulo de Sierpinski 
8 el patrén al que tienden todas las imagenes tras la aplicacion 
reiterada de nuestra fotocopiadora. Y si empezaéramos directa- 
mente con él las iteraciones sucesivas, lo dejarian invariante. El 
conjunto de Sierpinski es, por tanto, un «punto fijo» de nuestra 
aplicacion del folio sobre el folio. Y es a la vez un atractor, todas 
las imagenes tienden a él tras sucesivas iteraciones. 

En cuanto a la informacion, probablemente el lector también 
se habra dado cuenta de lo que ocurre: la aplicacién sucesiva de 
la fotocopiadora reductora multiple produce la pérdida progre- 
siva de informacién. Poco a poco vamos perdiendo los detalles 
de la imagen original, que, en el infinito, desaparecen totalmente. 
Esta pérdida de informacion esta relacionada con el hecho de 
que la fotocopiadora reductora miltiple es una transformacién 
contractiva: dados dos puntos cualesquiera del plano, la contrac- 
cién los acerca en la siguiente iteraci6n. 


LA GEOMETRIA FRACTAL DE LA NATURALEZA 


La palabra fractal fue inventada por Benoit Mandelbrot en 1975, 
inspirandose en el latin fractus, para designar la geometria de 
objetos similares al triangulo de Sierpinski. Mandelbrot nacié en 
Varsovia en 1924 en el seno de una familia judia. La inminente 
persecucién nazi llevé a su familia a buscar refugio en Paris en 
1936. Antes de que los nazis llegaran a la capital francesa, huye- 
ron de nuevo hacia el sur y se instalaron en un pueblo cercano a 
Limoges. Al acabar la guerra, de vuelta en Paris, el joven Benoit 
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apr (tion sc seed Ete Meeps nSecunla 
Poliiéenten, & ingresé en ewta ailtina, 
S Algunce aio antes, Honrl Polnearé, el gran matematico y fisi- 
86 francés nacido er) Nancy en 1854, mantuvo un debate abierto 
eon otros matemAticos sobre los fundamentos de las mateméati- 
as, HS] aleman David Milbert era el mas insigne de sus oponen- 
(és, Hilbert abogaba por fundamentar las matematicas sobre una 
s¢ puramente formal: construir todo el edificio de las matema- 
S partiendo de unos pocos postulados de los cuales se de- 
ducian todos los teoremas, usando para ello las reglas inobjeta- 
bles de la légica. A esta concepcion Poincaré oponia una vision 
én la que la intuicién era imprescindible. De Poincaré son estas 
palabras, pronunciadas en el Congreso Internacional de Mate- 
maticos celebrado en Paris en 1900: «De este modo la légica y la 
intuicién tienen cada una su propio papel necesario. Ambas son 
indispensables. La l6gica, que es la nica que puede dar certezas, 
es el instrumento de la demostraci6n. La intuici6n es el instru- 
mento de la invencién». 
Pero Poincaré era un hombre que trabajaba en solitario y casi 
no tuvo discipulos. Por ello no creé escuela. A pesar de que tuvo 
en vida un gran prestigio en su pais, a su muerte las matemati- 
on cas francesas no siguieron su estela. Al otro lado de la frontera, 
UCesIVOS 8 cmentos. Asi 5, al : David Hilbert si tuvo innumerables colaboradores y su influen- 
16 puntos. Sin embargo, su longi- ; ; 
de longitud 1/3”, por lo que la cia se extendid, desde Alemania, a toda Europa. Cuando Benoit 
names ‘en cade svn 82/9 Este valor tends Mandelbrot inicié sus estudios en la Escuela Politécnica, la mis- 
4 4, : ma donde se habia formado Poincaré mas de medio siglo antes, 
ae las concepciones filosdéficas de los mateméticos franceses esta- 
ban mas cerca de Hilbert que del célebre matematico de Nancy. 
Y Mandelbrot, un joven especialmente dotado para el dibujo y 
con una singular visi6n espacial, encontraba poco satisfactoria 
esa forma de concebir las matematicas. 
Un ejemplo del ambiente de las matematicas francesas en esa 
época lo da la definici6n del nimero x que aparece en un tratado 
de Nicolas Bourbaki, nombre imaginario con el que firmaba sus 
publicaciones un grupo influyente de matematicos. Bourbaki de- 
fine asi el numero 7m: «La funcion e (a) admite en todo punto de R 
una derivada igual a 2x, donde z es una constante positiva». No 
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ca rocesarioa ontender este definicion para triaginar el prado de 
jbsleeciOn al que se hebia ogado, Uno de los nimeros mAs eo. 
nocidos de todas las malematicas, ligado hisioricanente a la geo- 
metria, aparece desprovisto de cualquier connolacion intuitiva. 

Tal vez para huir de este ambiente, poco después de terminar 
sus estudios en Parfs Mandelbrot emigr6é tna vez mas, en esta 
ocasion a Estados Unidos, donde se asent6. Alli se integré en el 
centro de investigacién de IBM Thomas J. Watson Research Cen- 
ter, donde empezo a interesarse por todo tipo de problemas, des- 
de las subidas y bajadas del precio del algod6n hasta la frecuen- 
¢ia con que se producen los fallos en los circuitos electronicos. 

iin 1982 Mandelbrot public6 La geometria fractal de la na- 
huevaleza, una revisién expandida de un primer ensayo sobre el 
mismo Lema escrito en 1977. El texto reunia las incursiones del 
autor en diversos campos. En él hablaba de economia, de ga- 
laxias, de turbulencia, de la longitud de la costa de Inglaterra, de 
arle, de montafias y de muchas mas cosas. Mandelbrot sostenia 
que las curvas continuas y diferenciables tan queridas para los 
fisicos, aquellas que se pueden dibujar sin levantar el lapiz del 
papel, no eran la herramienta adecuada para describir muchos 
(le Jos patrones y formas geométricas que encontramos en la na- 
luraleza. En su libro Mandelbrot introdujo conceptos que hoy 
lia Son esenciales en el campo del caos determinista: las nocio- 
nes de conjunto fractal, dimension fractal y autosimilaridad. 

De entre los muchos ejemplos que describe en su obra, la lon- 
pitud de la costa de Inglaterra se ha convertido quiza en el para- 
digma por excelencia de objeto fractal. gCuanto mide la costa de 
Inglaterra? Podemos empezar tomando un mapa con una escala 
determinada, por ejemplo 1:2250000, y medir con una regla la 
longitud de la costa. Obtendriamos un determinado valor. Si re- 
petimos la misma operaci6n con un mapa con una escala mayor, 
1: 100000, por ejemplo, apareceran detalles y pliegues en la costa 
que no eran perceptibles en el primer mapa, con lo que la medi- 
da de la longitud total se incrementara. El caso es que, a medida 
que aumentamos el detalle con el que describimos la costa, 
su longitud aumenta, y parece que no hay un limite, por lo que 
la longitud real seria infinita. La situaci6n contrasta con la de 
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ine curva de In veometriy sonvenelona, Por ejemplo, podemoas 
Thacribly en ann elreunmferenela un poligono de 7 lados, cada uno 
de longltud d 

La Jongitud total de la cireunferencia la podriamos aproxi- 
mar entonces por” veces la longitud de los lados, esto es, nl. 
A medida que aumentamos el valor de n, el poligono se parece 
euda vez mas a la circunferencia, y la longitud total medida por 
este procedimiento se acerca mas y mas a la longitud conoci- 
da de la circunferencia, 27k, donde R es el radio (a Bourbaki 
ésia definicion de « como el cociente entre la longitud de una 
tircunferencia y su diametro seguramente le pareceria poco ri- 
purosa). La longitud de las curvas sencillas de la geometria se 
puede determinar sin ambigtiedad por este procedimiento. Pero 
este falla para un objeto como la costa de Inglaterra: la suma de 
los Segmentos en que hemos dividido el contorno crece indefini- 
damente a medida que estos son mas pequefios. 

Este hecho ya habia sido observado por algunos cientificos. En 
particular, por el mismo Lewis Richardson que encontramos en 
el capitulo primero, el que intent6 predecir el tiempo atmosférico 
con papel y lapiz. Su conclusién fue que las costas de los paises, 
sus fronteras y otras caracteristicas similares no tienen una lon- 
gitud definida, y que la longitud atribuible depende de la regla con 
que la midamos. Pero Benoit Mandelbrot supo reinterpretar estos 
datos dando un salto cualitativo. La costa de Inglaterra no tiene 
una longitud definida porque es algo m4s que una curva, aunque 
menos que una superficie: la costa es fractal. Y Mandelbrot se 
aplico a tratar de determinar, aunque fuera aproximadamente, la 
dimension fractal de las costas de varios paises. Esta dimensién 
result6 ser proxima a 1,5 en varios casos. Ademas, Mandelbrot 
ided un algoritmo para generar objetos fractales que asemejaban 
la costa de una isla o de un continente. En su libro se pueden ver 
algunos de estos objetos generados por ordenador. 

Cabria decir que La geometria fractal de la naturaleza dis- 
gustaba por igual a fisicos y matematicos. Para un matematico, 
Mandelbrot no demostraba nada. Para un fisico, Mandelbrot no 
explicaba nada. La costa de Inglaterra es fractal, de acuerdo, 
pero. {por qué?, preguntaria un fisico. Incluso aunque Mandelbrot 
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Li jas imagenes'4, 5 y 6 nos 
ieritramos en la. Gola de uno de los 
cubaliitos de mar. La propia-esiructura 
(le la espiral se-multiplica ea las 
liquras 6, 7 y'8. En el centra de 

la figura 7, emerge de nuevo la 
silueta inconttindible del conjunto 

(le Mandelbrot 
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proporclonara algoriimos para simula, {donde eataban las le 
yes generaios, Jos primeros principlog’ Sin embargo, la moderna 
leoria del enos y th conplejidad io es coneebible sin los con: 
eoplos introducidos por Mandelbrot. Mandelbret no demostraba 
jeorenas, Clerto, pero dibujaba, calewlaba, construia. Tampoco 
explicaba nada, pero deseribia, sugeria, proponia, encontraba re- 
glas alli donde otros solo veian desorden. 

La geometria fractal de ta naturaleza tuvo, y tiene, defenso- 
res y detractores. No es un libro para especialistas, pero tam- 
poco para no iniciados. Tampoco se gana al lector alguien que 
desde el prefacio reclama para si el honor de haber «fundado 
una nueva ciencia», como hace Mandelbrot. Pero es un libro que 
sorprende al lector a cada pAgina y que le hace pensar y replan- 
(earse muchas ideas preconcebidas. Y ha sido, y sigue siendo, 
un superventas. Treinta afios después de su publicacién Benoit 
Mandelbrot ha recibido todos los honores y medallas imagina- 
bles. Y su obra es, sin duda, uno de los hitos de la literatura cien- 
lifica contemporanea. 


VUELVE LORENZ 


Como vimos en el capitulo anterior, el meteordélogo Edward Lo- 
renz habia ideado un sistema de doce ecuaciones diferenciales 
que simulaban una especie de atmésfera «de juguete». Recorde- 
mos que una ecuacion diferencial expresa matematicamente el 
modo en que un sistema cambia continuamente en el tiempo. 
En esa atmdsfera de juguete descubrié Lorenz la sensibilidad a 
las condiciones iniciales de los sistemas cadticos: dos estados 
inicialmente muy pr6éximos daban lugar, al cabo del tiempo, a 
situaciones muy diferentes. 

Lorenz entendié que habia algo en la estructura matematica 
de su sistema que merecia estudiarse mas a fondo. La atmés- 
fera es lo que se llama un sistema disipativo: la viscosidad del 
aire tiende a frenar el movimiento, disipando la energia cinética 
del viento y convirtiéndola en calor, que, a su vez, se dispersa 
gracias a la conduccion, la conveccién y la radiacién. Un siste- 
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min Hisipalive que aio o908 alipeninde por inp Men.e eer de 
tnerigin fende al eqiittbria. Por ejemplo, un eolunmpie eb que le 
Hemos dads sole un inpulso va lrenindose poco a poco, debi- 
do al rovaniento, haste que lerming parandose en el punto mas 
bajo. Podemos deseribir su estado mediante su posicion (con el 
Anpulo 0) y su velocidad » (figura 3). 

Si representamos la trayectoria del columpio en el espacio de 
las fases, esto es, representando su posicion en el eje de absci- 
§as y Su velocidad en el eje de ordenadas, obtenemos la imagen 
izquierda de la figura 4. La curva que describe el movimiento del 
columpio tiende a un punto, que representa al columpio parado 
@n su posicién de equilibrio. El columpio es un sistema dinami- 
co, y el punto mas bajo, su posicidn de equilibrio, es lo que se co- 
noce como un atractor porque, independientemente del impulso 
inicial que demos al columpio, todas las trayectorias que dibuje- 
mos en el espacio de las fases acaban en el punto de equilibrio. 

Supongamos ahora que decidimos impulsar el columpio re- 
pularmente, con la frecuencia adecuada para mantener una os- 
cilaci6n permanente. Si al inicio el columpio esta en reposo, la 
trayectoria seguida en el espacio de las fases es la que muestra 
la imagen derecha de la figura 4. Esta trayectoria representa el 
hecho de que el columpio oscila cada vez con mayor amplitud 
hasta alcanzar un estado estacionario, en el que oscila sin va- 
riar su altura. El impulso que damos en cada periodo es justo el 
que compensa la pérdida de energia por la friccién de la cuerda 
con el punto de sujecién, La elipse que representa el movimiento 
sostenido en el espacio de las fases es también un atractor; en el 
argot matematico se conoce como ciclo limite. 

En un caso genérico habra estados del sistema que se repi- 
tan, o que sean aproximadamente los mismos una y otra vez, 
acercandose mas y mas. Estos estados pertenecen a un conjunto 
restringido, el atractor. 

El columpio que impulsamos en cada oscilacion es un sistema 
disipativo, pero que constantemente recibe un aporte de energia 
desde el exterior. De nuevo los fisicos y matematicos tienen una 
expresién especifica para este tipo de sistemas: se trata de un 
sistema disipativo con un forzamiento externo. 
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‘Thiyortorla de un columpio en el espacio de las fases. Izquierda: sin forzamiento; derecha: impulsado con 
Ih Tecuencia adecuada. 
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angulo y la velocidad 
Ghee es = marcados). : 
t 


Loreng sabia die au abméefern de jugtiele era im later dist 
palivo foraznd, in el cave de Jo atmosfera, la energin externa la 
Sporta el sol. Localmente, el ol ealienta él suelo y las masas de 
aire Circundantes tienden a ascender. La diferencia de calenta- 
miento entre el ecuador y los polos alimenta también los siste- 
mas globales de circulacién atmosférica. Lorenz esperaba que la 
nimésfera, como sistema disipativo forzado que era, tendiera a 
un atractor, aun Conjunto restringido de estados. Pero sus resul- 
(aclos mostraban claramente que ese atractor no era ni un punto 
ni un ciclo limite, porque sus variables no tendian a valores es- 
(acionarios ni a oscilaciones regulares. Ni siquiera se trataba de 
ima combinacion de ciclos de diferentes frecuencias. Pero su 
atmosfera de doce variables era todavia demasiado complicada. 
Lorenz intent6 simplificar ain mas sus ecuaciones sin perder lo 
esencial: la variabilidad de las soluciones, que no se tratara de 
soluciones estacionarias, periddicas o facilmente predecibles. 


El primer atractor extrafo 


En 1963 Edward Lorenz publicé en Journal of Atmospheric 
Sciences un articulo titulado «Flujo no periddico determinista». 
En él presentaba un modelo de tres ecuaciones diferenciales con 
solo tres variables. Este medelo no era ya una atmésfera simpli- 
ficada, sino que representaba la convecci6n de Rayleigh-Bénard 
en un rollo convectivo similar a los que Libchaber estudiaba en 
su caja de helio (figura 5). 

Era una estructura convectiva espacialmente confinada, ya 
que el liquido no podia salir de la caja en la que estaba encerra- 
do, pero dependiente del tiempo, porque la intensidad y la di- 
recciOn del movimiento podian variar. De las tres variables, una 
daba cuenta de la velocidad maxima del liquido dentro del rollo; 
otra, de la diferencia de temperaturas entre la columna de fluido 
ascendente y la descendente, y la otra variable representaba el 
gradiente vertical de temperatura. E] sistema contenia tres para- 
metros relacionados con las propiedades del liquido, las dimen- 
siones de la caja y la diferencia de temperaturas entre los limites 
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apariecclt iapon, peer Senate ne Mibrinoa rc icbalsa nace 
sarios para ofrecer un reo compertamiento dinimico, Recorde- 
mos que en una funclon no teal la funeién de wna sume no es 
la suma de las funciones. Ei] ejemplo mas senciilo de funcién no 
lineal es x’, ya que (#,+2,,)" = a7*+2,". La presencia de términos 
no lineales es iriprescifdibla en el caos determinista, 

I} sistema de Lorenz reproducia cualitativamente el compor- 
tamiento experimental observado a medida que se aumenta la 
diferencia de temperatura. Cuando esta es baja, el liquido no se 
mueve y la situacion es de equitibrio estable. En el espacio de las 
fases el estado del sistema esta representado por un punto en el 
origen de coordenadas. Que se trate de un equilibrio estable sig- 
nifiea que ante cualquier perturbacion el sistema regresa a él. Es 
dec, si se produce externamente, por agitacién o por cualquier 
Olra causa, una pequefia escilacién de temperatura o de veloci- 
dad, el sistema vuelve a la situiacién de equilibrio, en este caso, al 
reposo. Al aumentar progresivamente la temperatura se produce 
ine primers bifureacion, un cambio cualitativo en el comporta- 
tiiento observado, y el liquido comienza 4 moverse formando un 
ralio: sube por un lade y baja por el otro. Si la temperatura no 
ea muy alta, la velocidad no varia en el tiempo y esta es también 
una situacién de equilibrio estable. Pero hay dos posibilidades 
pguivalentes: que el rollo gire en el sentido de las agujas del reloj 
© en sentido contrario. Esta situacién viene representada em el 
espacio de las fases por dos puntos simétricos. 

Lorenz comprebé gue, al aumentar suficientemente la dife- 

rencia de temperatura, el sisteraa de ecuaciones comenzaba a 
presentar la variabilidad que ya habfa observado en su anterior 
_ modelo atmosférico. Como ejemplo, la figura 6 muestra la evolt- 
eién en el tiempo de la. velocidad maxima dentro del rollo para 
determinade valor de la diferencia de las temperaturas superior 
e inferior. El signo positivo de la velocidad representa un rollo 
que se mueve en el sentido de las agujas del reloj. El signe nega- 
tivo representa un movimiento en e! sentido contrario. 
Observamos que la velocidad toma unas veces valores posi- 
tivos y otras valores negativos, cscilando un nimero cambiar- 
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Eequema ie un salle convective 
CoMmG #) alimulda por Lorenz 
On Sus @cuaclones, El liquide 
#e Callenta por abajo. Todo el 
dominic esta lleno de Ilquido, 
que asciende por un lado 

y desciende por el otro. 


Fuente de calor 
FIG. 6 


Velocidad 


Tiempo 


Evolucién temporal de la velocidad maxima del liquido en el sistema de Lorenz para un conjunto escogido de 
parametros. La velocidad oscila irregularmente alrededor de dos valores de signo contrario, El signo positivo 
significa que el fluido se mueve girando en el sentido de las agujas del reloj y el signo negativo, que se mueve 
en el sentido contrario. Las lineas discontinuas horizontales muestran los valores en torno a los cuales oscila 
la velocidad. El tiempo y la velocidad estan expresados en unidades arbitrarias. 
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te de veces en tome a dos valores, cercarog e #7, TP slavenie 
erocuerdae que exiatian esos puntos de equilibria, que repre: 
sentaban moviniientos estacionarios del (hide, y da vuellas al- 
rededor de ellos, pero sin Jlegar a esiabilizarse en ninguno. [1 
aumento de la diferencia de ternperaturas y, por tanto, el ritmo 
al que inyectamos energia en el sistema, es la causa de esta ines- 
tabilidad. 

La figura 7 es una representacién en tres dimensiones de una 
trayecloria en el espacio de las fases durante un gran interva- 
lo de tiempo. El sistema no acaba en un punto de equilibrio ni 
en un ciclo limite, sino que es atraido a una situacién en la que 
vaga indefinidamente en torno a una zona. limitada del espacio. 
Ha caido en un atractor extravio. Es un atractor porque todas 
las trayectorias iniciadas en sus cercanias acaban cayendo en 
él. Pero se le conoce con el apelativo de extraro porque no tie- 
ne ninguna de las formas geométricas convencionales: no es un 
punto ni una circunferencia ni una curva o superficie conocida. 

El atractor extrafio representa un estado estable del sistema 
en cuestion. Las trayectorias van hacia él y, una vez que llegan, 
permanecen alli. En un sistema cadtico que muestre un atractor 
extrafio, un estado inicial cualquiera no tiene por qué pertenecer 
al atractor, pero al evolucionar se le acerca cada vez mas y ter- 
mina en él. Una de las razones por las que Richardson fracaso en 
su primer intento de calcular numéricamente la evolucién de la 
atmésfera, como vimos en el capitulo anterior, es que el estado 
inicial que él escogié no estaba suficientemente caracterizado. 
La dinaémica atmosférica ya esta en un atractor, por muchas di- 
mensiones que tenga y por muy dificil que sea de caracterizar. 
Un conjunto cualquiera de presiones, temperaturas y demas va- 
riables no tiene por qué pertenecer a un estado posible de la 
atmosfera en ese atractor. Calcular la evoluci6n de la atmésfera 
a partir de unos datos cualesquiera lleva, en general, a una situa- 
cidn transitoria, que no se corresponde con la evolucion real de 
la atmoésfera. Si un monstruo gigantesco llegara a la Tierra y so- 
plara brutalmente sobre ella, sacaria a la atmosfera de su atrac- 
tor. Eso originaria un periodo transitorio mas o menos largo, con 
grandes oscilaciones de vientos, presiones y temperaturas, hasta 
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que la dindmica propia del sistema volviera a traer a la atmésfera 
a su régimen habitual. Richardson, sin querer, hizo el papel de 
monstruo, introduciendo unos valores iniciales que originaban 
fluctuaciones no realistas. Este es un problema que los meteor6- 
logos han sabido resolver. Los valores observados no son siem- 
pre suficientemente precisos y pertenecen a una red irregular de 
observatorios. Antes de introducir los datos observados en los 
modelos numéricos, los meteorélogos los someten a un proceso 
de «inicializacién» que produce unos valores iniciales que, por 
decirlo asi, estén ya dentro del atractor, con lo que se suprimen 
oscilaciones espurias. 

EJ atractor de Lorenz es el primer atractor extrafio encontra- 
do, pero el hecho de que Lorenz fuera meteordlogo y publicara 
sus resultados en una revista de meteorologia hizo que pasara 
desapercibido en otros campos durante muchos afios. No fue 
hasta mediados de la década de 1970 cuando el descubrimien- 
to de Lorenz se dio a conocer en otros ambitos y se xeconocié 
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su relovanelo, La oxprealon afradtor evtraie tie aeutiada por el 
bela David Tuelle ¥ él holandés Mors Tikens en 1971, Aunque 
(6 Un nombre «extrafoe pata iin objeto matemalies, e¢ incluso 
Jn denominacion ha side eriiicada por alunos matematicos que 
Plensan que este Lips de atractor resulla extratio solo a quienes 
no lo entienden, la denominacion se ha impuesto tanto en la lite- 
rilurn eapecializada como en los textos divulgativos. 


LA RECETA DEL CAOS 


iin un atractor extrafio las trayectorias tienden a juntarse en una 
direccién, porque el atractor atrae todas las trayectorias de 
una cierta regién del espacio de las fases, pero también tienden 
4 separarse, porque si no no observariamos la sensibilidad a las 
condiciones iniciales. En el atractor de Lorenz vemos estas dos 
caracteristicas. La forma de mariposa es aproximadamente una 
superficie de dos dimensiones, aunque tiene cierto espesor, y 
cualquier trayectoria cuyos valores iniciales estén lejos de esta 
superficie casi bidimensional tiende a ella. Esto representa la 
contraccién. Pero, a la vez, dos trayectorias inicialmente muy 


juntas pueden divergir rapidamente y dar vueltas cada una alre- 


dedor de uno de los dos centros del atractor. Ahi vemos la sensi- 
bilidad a las condiciones iniciales. Por esas dos razones, la rece- 
(al para construir un atractor extrafio es relativamente sencilla. 
Se trata de plegar y estirar, plegar y estirar. 

E] sistema de ecuaciones de Lorenz tiene tres variables. Es 
el nGmero minimo de variables que se necesita para producir 
caos-en un sistema que dependa del tiempo de forma continua, 
esto es, en un sistema de ecuaciones diferenciales. Este hecho 
fue comprendido ya por Henri Poincaré, que en sus estudios de 
la mecanica celeste vislumbr6 que los sistemas de ecuaciones 
diferenciales no lineales de tres o mas variables podfan llegar a 
tener soluciones muy complicadas, inabordables con los méto- 
dos de la época. En su camino por comprender los sistemas no 
lineales, Poincaré inicié el estudio cualitativo de las ecuaciones 
diferenciales. Renunciando a resolver con todo detalle las ecua- 
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ones, ded meétodos y hertamients pare entender la eatruetarn 
Wlobal de las soluciones, Se interas6 por cuestiones tales como 
CuAntos puntos de equilibrio habia, al eran estables 6 no, si habla 
soluciones periddicas, etcetera. 

La dificultad progresiva en abordar las ecuaciones diferen- 
ciales no lineales con los métodos del analisis mateméatico tra- 
dicional llevé a Poincaré del andlisis a la geometria, y de ahi a 
la topologfa, de la que es considerado uno de los creadores. La 
1opologia se interesa no por la forma concreta de los objetos 
feorétricos, sino por algunas caracteristicas cualitativas tales 
como el numero de piezas que los forman o el numero de agu- 
jeros que tienen. Dos figuras son topolégicamente equivalentes 
si podemos deformar una de ellas, sin cortar ni pegar trozos, de 
forma continua, hasta hacerla coincidir con la otra. Para un to- 
p6logo un tridngulo, un cuadrado y una circunferencia son lo 
mismo. 

El matematico estadounidense Stephen Smale hizo el camino 
inverso al recorrido por Poincaré. Inicié su carrera matematica 
trabajando en topologia. Aficionado a coleccionar minerales y 
a navegar a vela, Smale hizo sus grandes descubrimientos ins- 
pirado por el mar en Rio de Janeiro. Aunque, a decir verdad, el 
atractivo de Rio no era solo la playa: allf estaba el matemati- 
co brasilefio Mauricio Peixoto, autor de un importante teore- 
ma sobre la estabilidad de las érbitas de un sistema dindmico, 
con quien establecioé un fructifero intercambio de ideas. Smale 
resolvi6 parcialmente uno de los problemas que habia dejado 
abierto Poincaré: la conjetura de Poincaré. Poincaré habia con- 
jeturado que toda figura geométrica que no tuviera agujeros 
y que estuviera hecha de una sola pieza era topolégicamente 
equivalente a una esfera. En el espacio ordinario la conjetura 
de Poincaré es facilmente reconocible: cualquier poliedro, por 
ejemplo un cubo, puede deformarse continuamente hasta for- 
mar una esfera. Poincaré generalizaba este enunciado a obje- 
tos de cualquier dimension (aunque no podamos hacernos una 
imagen intuitiva de estos objetos, los matematicos idean esferas 
y poliedros en espacios de mas de tres dimensiones). Pero no 
fue capaz de demostrarlo, y lo dej6 como problema abierto en 
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uno de sus TILTON eseritos sobre topologin, Pues bien, en 161 
Sule consigiié demositw lo canjotira de Pomeant pari ole 
jetos de 7 6 mis dimensiones. Wae hallaggo le valid a Stephen 
Simle In medalla Melds, uno de los mas presliglosos galardo- 
noes que puede recibir un matemaAtico, que le fue entregada en el 
songreso Internacional de Matematicos en Moseu en 1966. La 
conjetura de Poincaré para dimensiones menores se mostré un 
hueso duro de roer, y no fue hasta los primeros aiios del siglo 
xx cuando el ruso Grigori Perelman consiguié demostrarla para 
tres dimensiones. Ello también le valié a Perelman la medalla 
Fields, a la que renunci6. 

Temiendo que su éxito con la conjetura de Poincaré agotara 
sus posibilidades como matematico dentro del campo de la to- 
pologia, Smale decidié cambiar de asunto y se interes6 por los 
sistemas dinamicos. Como hemos dicho, en un sistema dinami- 
co continuo, en el que el tiempo fluye continuamente y puede, 
por tanto, tomar cualquier valor, el nGmero minimo de variables 
necesarias para producir caos es tres, y es también el minimo 
necesario para producir un atractor extrafio. Recordemos que 
en un sistema dindmico discreto registramos la evolucién del 
sistema en intervalos fijos de tiempo: dias, afios, generaciones, 
etcétera. En el caso de los sistemas dinamicos discretos, como la 
aplicacion logistica estudiada por May y Feigenbaum, podemos 
conseguir crear un atractor extrafio con dos dimensiones. Y eso 
es lo que hizo Smale. En 1967 ideé un sistema dinamico discreto 
que contiene todos los ingredientes del caos. A este sistema se le 
conoce como la herradura de Smale. 


Smale idea una herradura 


La idea esencial de la herradura de Smale es tomar un rectaéngu- 
lo, estirarlo en una direcci6én, doblarlo hasta formar una especie 
de herradura y luego quedarse con la interseccién de la herradu- 
ra y el rectaéngulo original. Después se trata de repetir la opera- 
cién una y otra vez, estirando y doblando, estirando y doblando, 
tal y como haria un panadero con la masa del pan. 
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Sucesién de operaciones 
que conducen a la 
construccidn de la 
herradura de Smale. 


Vedmoslo graficamente en la figura 8. La herradura de Smale 
consta de dos transformaciones sucesivas. Consideremos la re- 
gidn de la figura 8a, que consta de un rectangulo en el que hemos 
marcado dos bandas sombreadas como referencia. En primer 
lugar, contraemos toda la regién en la direccién horizontal has- 
ta hacer que su anchura sea menos de la mitad de la anchura 
inicial. A la vez estiramos en la direccién vertical, de forma que 
el area total se mantenga. Con ello obtenemos una banda como 
la de la figura 8b. A continuacién, doblamos la banda obtenida 
para formar una especie de herradura, como muestra la figura 
8c. De la herradura asi obtenida nos quedamos solo con la parte 
que resulta de hacer la interseccién con el rectangulo original. El 
resultado de la operacién completa es que hemos llevado cada 
punto del rectangulo original a otra posicién dentro del mismo 
rectangulo. ; 

En la figura 8d vemos el resultado de realizar dos veces con- 
secutivas la misma operacidn. Las dos bandas horizontales ori- 
ginales se han transformado en cuatro bandas verticales. Si vol- 
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estrechas, Doblando y estirando @ hacer mucho énfasis 
t as veces obtenénios una gran 
-eantidad de bandas estrechas. Vamos a 
generando asi una estructura fractal: ™Matematicas. 

las bandas no llenan todo el rectangu- 

Jo, por lo que no forman un objeto de 

dimensién dos, pero son més que una linea, por lo que su dimen- 
sién puede ser mayor que uno. 

El efecto de estirar y doblar es semejante al que se produce 
cuando amasamos una masa de pizza. Si espolvoreamos de ha- 
rina la masa, doblar y estirar muchas veces va extendiendo la 
harina incorporada por toda la masa, de manera mas o menos 
homogénea. Particulas de harina que estaban muy cerca al espol- 
vorear pueden acabar en extremos diferentes de la masa: la ope- 
racion presenta sensibilidad a las condiciones iniciales. La herra- 
dura de Smale tiene pues dos de las sefias del caos: fractalidad y 
sensibilidad a las condiciones iniciales. 


El atractor de Hénon: la herradura codificada 


En 1976 el astrénomo francés Michel Hénon introdujo un sis- 
tema dinamico discreto bidimensional que contenfa en una 
formula sencilla todos los elementos de la herradura de Smale. 
Hénon habia llegado a interesarse por los sistemas dinamicos a 
partir de sus estudios sobre la formacidn y la estabilidad de los 
cimulos globulares, agrupaciones densas de miles de estrellas 
ligadas gravitatoriamente. En 1976 Hénon se incorporé al Ob- 
servatorio de Niza y allf, advertido de los estudios de Lorenz 
y otros sobre los atractores extrafios, se decidié a usar su ex- 
periencia para tratar de visualizar un atractor de la forma mas 
sencilla posible. 

La transformacién que Hénon idedé produce un atractor ex- 
trafio bidimensional que se conoce como atractor de Hénon. 
Map, El sistema dinamico de Hénon condensa en ecuaciones mate- 
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mAiions lo# elementos condeptunios de ln herradura de Smale; 
doblar y estirar, Dade win punto del plano de coordenadas (x, 
V,), las ecuaciones de Hénon proporcionan Ins coordenadas de 
aire punta (wy). A partir de abl, la aplieacion sucesiva de la 
fransformacion de Hénon produce una secuencia de puntos (2,, 
Uy) Ce) UV), py UW), My Uy). Que se pueden entender como una 
trayectoria en el plano, Después de varias iLeraciones, los suce- 
sivos puntos tienden a un atractor que no es ni un punto ni una 
curva, sino un alractor extrano. Este alraclor se representa en 
Ja figura 9. 

Ml atractor de Hénon es una figura fractal. Su dimensién es 
de 1,26, lo que quiere decir que no Hlega a rellenar un area com- 
pletamente, pero es algo mas que una simple linea que se curva 
sobre si misma. Si hacemos una ampliacién de cualquier zona 
del atractor se nos revelan nuevas estructuras, similares unas 
n olras a diferentes escalas: es la autosimilaridad. La figura 10 
muestra dos ampliaciones sucesivas. Se observa que, a medida 
que ampliamos la imagen, aparecen nuevas hebras que a una es- 
enla anterior no eran observables. Estas hebras se abren suce- 
sivamente, ampliacion tras ampliacién. En principio el proceso 
se repite hasta el infinito. En la practica, el nimero de hebras 
que consigamos revelar esta limitado por la precisién del or- 
denador utilizado y el tiempo de cdlculo. La autosimilaridad es 
una caracteristica tipica de los atractores extrafios y aparece 
¢n otros objetos que ya hemos descrito, como el diagrama de 
Feigenbaum o el triangulo de Sierpinski. 


El atractor de Rossier 


Para terminar la presentaci6én de los atractores extrafios, presen- 
taremos un atractor mas, obtenido esta vez en un sistema dina- 
mico continuo, esto es, descrito por ecuaciones diferenciales: el 
atractor de Réssler. Este atractor fue descubierto por el cientifico 
alemén Otto Réssler en el contexto de ciertas reacciones quimi- 
cas. Se obtiene de un sistema de tres ecuaciones diferenciales, y 
lo vemos representado en la figura 11. De nuevo tenemos los dos 
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El atractor de Hénon es un atractor extrafio. Se trata de una figura fractal compuesta por infinidad de fibras que 
se contienen unas a otras. 


; \ S : 
i [ - { An \ \ L\ 4 : x 


ul n 
0,3 0,35 0,4 0,45 0,5 0,55 0,6 0,65 


04 02 0 02 04 06 08 


Ampliaciones del atractor de Hénon. Notese la escala en los ejes y comparese con la figura 9. 
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Ingrediontes neceyariog: Gatlirimients, que xe produce a lo largo 
del ajee, ¥ plagada, que se produce en el plano ay, 


MIDIENDO EL CAOS 


iin ¢1 primer capitulo vimos que tanto el sistema solar como la 
nimésfera se Comportaban cadéticamente. Resulta contradicto- 
rio que podamos predecir un eclipse con siglos de antelacién 
y seamos incapaces de predecir el tiempo a dos semanas vista. 
Parece que un sistema es més cadtico que otro. gCémo podemos 
expresar esa diferencia de manera cuantitativa? Fisicos y mate- 
mAaticos han ideado varias herramientas que ayudan a cuantificar 
el caos en un sistema. En este epigrafe vamos a presentar tres 
de ellas: los exponentes de Liapunov, la informacion y las dimen- 
siones fractales. 

Los exponentes de Liapunov, llamados asi en honor al ma- 
fematico ruso Aleksandr Liapunov (1857-1918), cuantifican la 
skénsibilidad alas condiciones iniciales. Consideremos un siste- 
mia dinémico continuo y, en él, dos trayectorias que se inician 
muy Cerca una de otra, de tal forma que la distancia entre ellas 
es 6, (figura 12). Al evolucionar, la distancia d entre estas dos 
trayectorias variara con el tiempo. En los instantes iniciales la 
divergencia o convergencia de ambas trayectorias es exponen- 
cial, y el factor de crecimiento o decrecimiento de la distancia 
recibe el nombre de exponente de Liapunov, A. 

Si el exponente A es negativo, las trayectorias se acercan; si 
es positivo, se separan. Un exponente de Liapunov positivo es 
sintoma inequivoco de caos determinista, porque significa que 
puntos que inicialmente estan cerca se separan rapidamente: 
esto es la sensibilidad a las condiciones iniciales. E] inverso del 
exponente 1/A es un tiempo, que representa la escala tempo- 
ral a la que las divergencias entre las trayectorias se hacen no- 
tables. Para el sistema solar, ese valor es del orden de varios 
millones de afios (el valor concreto depende de qué parte del 
sistema analicemos en detalle). Para el modelo atmosférico de 
Lorenz que presentamos en el capitulo 1, es del orden de sema- 
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Atractor de Réssier. 


Exponentes de Liapunov. 8, es una distancia muy pequefia que inicialmente separa las dos trayectorias; des 
la separaci6n entre estas dos trayectorias al cabo de cierto tiempo. En este caso se ha representado un 
exponente de Liapunov positive, ya que las trayectorias se separan. 


ZQUE TIENE DE EXTRAS UN ATRACTOR EXTRANO? 


107 


108 


nie, PD) caon atmoaférico se reveli pal en une eaenla de tiempo 
Mee TA Paquelsa, 

Loa niraclores extranos tenon una Maselinante relacion con el 
canceplo de informacion, Una seial periodica es perfectamen- 
te predeeible, Si conoecemos un ciclo podemos predecir con 
exactilud cual sera el nivel de la sefial mucho tiempo después, 
ya que basta con repetir el ciclo tanlas veces Como sea necesa- 
rio, Pero nna senal caédtica ¢s impredecible a largo plazo, y el 
comportamiento del sistema durante un cierto lapso de tiempo 
no nos da informacion de lo que va a ocurrir después. A medida 
que el sistema evoluciona vamos, por lo tanto, perdiendo infor- 
macion. Esto ocurria, por ejemplo, con nuestra fotocopiadora 
reductora multiple: los rasgos de la imagen inicial iban difumi- 
nandose a medida que realizabamos las iteraciones, para des- 
vanecerse por completo en el infinito. Pero también hay cierta 
creaciOn de informacion en un atractor extrafio. Cualquier tra- 
yectoria dentro del atractor lo recorre una y otra vez, visitando 
continuamiente regiones cercanas pero siempre diferentes. Asi, 
a medida que seguimos una trayectoria cualquiera, se va reve- 
lando progresivamente la estructura del atractor. En este senti- 
(lo, el atractor puede ser visto como un generador continuo de 
informacion. 

Los atractores extrafios son objetos fractales y, como tales, 
se caracterizan por su dimension fractal. Pero hay una diferen- 
cia entre objetos fractales como el triangulo de Sierpinski y un 
atractor como el de Lorenz. El tridngulo de Sierpinski es un ob- 

jeto homogéneo. Todas sus partes son equivalentes, y el atractor 
tiene el mismo aspecto en cualquier region en la que nos fijemos. 
Pero un atractor extrafto es el resultado de un proceso dinami- 
co. Hay zonas que son visitadas con mas frecuencia que otras. 
En una simulacién por ordenador, las zonas més visitadas del 
atractor se resolveran con mayor precision que las menos fre- 
cuentadas, por lo que la determinacién de la dimension fractal 
puede resultar dificil. Estas circunstancias han levado a definir 
diferentes dimensiones, que se utilizan como herramientas que 
caracterizan no solo la geometria, sino también la dinamica del 
sistema. 
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EL EFECTO MARIPOSA 


«EM aleteo de una mariposa en Brasil, Zorigino un tornado en 
Texas?» Ese era e) titulo de una conferencia que Edward Lo- 
renz promuncié en Washington en 1972 durante una reunion de 
ln Asociacién Estadounidense para 


el Avance de la Ciencia. Lorenz esco- Puede ocurrir que pequenas 
pid el titulo de su conferencia con la diferencias en las condiciones 
intencién explicita de que resultara —jnjcigles engendren grandes 


provocativo, sabiendo lo improbable 
que era una respuesta afirmativa a 
su pregunta, y es seguro que no pre- 


diferencias en los fenodmenos 
finales. Un pequefio error en las 


vid —no podia preverlo— las conse- primeras produciria un enorme 


cuencias de su elecci6n. Porque el error en los ultimos. 
llamado efecto mariposa ha trascen- 
dido el mundo de Ja meteorologia y 
de la ciencia para pasar a ser un lu- 
gar comtn y dar nombre o aparecer en peliculas, canciones fe) 
series de television. Lo que Lorenz queria plantear, en un ambito 
puramente meteoroldgico, era hasta qué punto la sensibilidad 
a las condiciones iniciales de los sistemas cadticos dificulta la 
prediccién meteoroldégica, y ese era el tema de su charla. Con el 
aleteo de la mariposa queria ilustrar esa pequeha perturbacion 
que puede hacer que dos condiciones iniciales aparentemente 
iguales no lo sean en realidad. Si la atmésfera es extremada- 
mente sensible a las condiciones iniciales, la capacidad de pre- 
decir el tiempo sera siempre muy limitada. Desde que Lorenz 
pronuncié su conferencia, los sistemas de prediccion se han 
perfeccionado hasta puntos insospechados y hoy en dia las pre- 
dicciones a menos de una semana vista son muy fiables. Pero 
algunas de las cuestiones que Lorenz planteaba siguen vigentes 
y la misma dindmica cadtica de la atmésfera dificulta las pre- 
dicciones a largo plazo, por mas que se mejoren los sistemas de 
calculo y de obtencién de datos. 

Aunque una mariposa en Brasil no pueda provocar un tornado 
en Texas, el efecto mariposa ha quedado como metafora de la 
sensibilidad a las condiciones iniciales, una caracteristica que 
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PreHenlA en mayor o Menor medida todos 168 slstenmas cadtieos 
deterministas, 

Aunque ei nombre de Loreng ha quedado vineulado al efecto 
Mariposa y la senaibilidad a las condiciones iniciales, en reali- 
dad se trata de un redescubrimiento, porque ya Poincaré lo habia 
formulado, casi en los mismos términos, en su obra Clencia y 
metodo: 


4Por qué los meteordélogos tienen tantas dificultades en pre- 
decir el tiempo? 4Por qué la lluvia o las tempestades parecen 
llegar por azar, de forma que la gente encuentra natural rezar 
para que Nueva o haga buen tiempo, mientras que encontra- 
rian ridiculo rezar para que se produzca un eclipse? Vemos 
que las grandes perturbaciones se producen generalmente en 
las regiones donde la atmésfera esta en equilibrio inestable. 
Los meteorélogos ven bien que este equilibrio es inestable, 
que un cicl6on va a originarse en alguna parte. Pero ;dénde?, 
80n incapaces de saberlo. Una décima de grado de mas 0 de 
menos en un punto cualquiera y el ciclon estalla aqui y no alla 
y extiende su destruccién en zonas que de otra forma se hu- 
bieran salvado. [...] Aqui de nuevo encontramos el contraste 
entre una causa minima, inapreciable para el observador, y 
efectos considerables, que son a veces terribles desastres. 


Gracias al uso generalizado de los ordenadores, hoy tenemos 
multitud de ejemplos de esta sensibilidad a las condiciones ini- 
ciales. Los trabajos de Lorenz ilustran un comportamiento que 
Poincaré habia predicho mas de medio siglo antes. 
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Caos clasico y caos cuantico 


Los atractores extrafnos aparecen en sistemas 

en los que la energia no se conserva. Pero hay 
sistemas fisicos en los que la energia se conserva. 
Es el caso del sistema solar, pero también 

del mundo de los atomos. En ellos el caos 
determinista se presenta bajo otro aspecto. 


Ya sabemos que tanto el sistema solar como la atmésfera son 
sistemas cadticos. Sin embargo, hay una diferencia fundamental 
entre ambos: la atmésfera es, como hemos visto en el capitulo 
anterior, un sistema disipativo, mientras que el sistema solar no 
lo es, al menos en una primera instancia. La viscosidad del aire 
tiende a frenar los ciclones y torbellinos, disipando su energia 
mecanica en forma de calor. Por ello, la atmésfera necesita un 
aporte externo de energia, que recibe del Sol, para mantener su 
actividad. Pero el caso del movimiento de un planeta es distin- 
to. Un planeta orbita alrededor del Sol viajando por el espacio 
vacio. No hay ninguna fricci6n que frene su movimiento y el 
astro describe una y otra vez su 6rbita elfptica alrededor del Sol 
sin necesidad de aportes externos de energia. La velocidad del 
planeta no es constante a lo largo de su 6rbita, sino que varia 
un poco. Cerca del perihelio, cuando esta mas cerca del Sol, su 
velocidad es mayor. En el afelio, cuando esta mas alejado, su mo- 
vimiento se ralentiza. La energia mecanica del planeta consta 
de dos términos: la energia cinética, debida exclusivamente a 
su movimiento y proporcional al cuadrado de su velocidad, y 
la energia potencial, que es menor cuanto mas cerca esta del 
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Sol Aunque cada uni de estas energion yaa a lo large de In 
OFDIUG, au Sui, In energin mecyien (oll, ae maniene conslane 
te, Cuando un distemn meeinico se comport de esin forma se 
dee que es un s7s/emea conserpaliyo, iin {sien ledrica este Lipo 
de sistemas recibe el nombre (eenico de sistemas hamilionia- 
nos, en honor del hsico y malematico trlandés William Hamilton 
(1805-1865). 

fur la vida cotidiana no hay sistemas hamiltonianos estric- 
lamente hablando, La friecién y otros mecanismos disipativos 
son omumipresentes en la naturaleza, por lo que wna descripcién 
hamiltoniana de un sistema mecanico es siempre aproximada. 
Podemos describir relativamente bien, por ejemplo, el vuelo de 
una pelota de baloncesto hacia el aro como un sistema hamilto- 
nano, Pero esa descripcién se revela totalmente insatisfactoria 
si queremos explicar la trayectoria de un bal6n de futbol cuan- 
do, en un lanzamiento de falta, el jugador le da efecto para supe- 
rar la barrera. En este caso no se pueden despreciar las fuerzas 
vjercidas por el aire, entre las que se incluye la friccio6n debida 
a la viscosidad. Igualmente, podemos imaginar un péndulo que 
oscila libremente como un sistema hamiltoniano, pero la ener- 
pia se conserva solo aproximadamente, si consideramos unas 
pocas oscilaciones, porque la friccién del aire o el rozamiento 
de la cuerda con el punto sobre el que pivota lo detienen poco 
a poco. 

Pero los sistemas hamiltonianos son de especial interés tan- 
to en las pequerias como en las grandes escalas espaciales. Ya 
hemos visto que el movimiento de un planeta, y por tanto el del 
sistema solar como un todo, es conservativo. En realidad, esto 
tampoco es estrictamente cierto porque hay algunos mecanis- 
mos disipativos, como los efectos de marea, pero estos son des- 
preciables en muchos casos de interés, y sus consecuencias son 
perceptibles solo a muy largo plazo. En el otro extremo de la 
escala, en la escala atémica, la dinamica de los electrones dentro 
del atomo tiene lugar también en el espacio vacio y, segtin las 
leyes de la mecanica cuantica, en ausencia de mecanismos disi- 
pativos. Dentro del atomo, la energia del electrén es suma de la 
energia cinética y la energia potencial eléctrica, ya que la fuerza 
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de atraecion entre el ndeles alomleo y los electronen en eléetr 
ca, Por 16 tanto, los dlomos tambien a6n sistemas hamillonianos. 

In 1608 sistemas disipativos cl cuos déeterminista tiene su mejor 
expresién én la nocién de atractor extrafio, un objeto geomeé- 
irico singular, caracterizado por tener una estructura fractal y 
autosimilar, En los sistemas hamiltonianos el caos se manifiesta 
de una forma diferente. 


LECCIONES DE UN PENDULO 


E:1 paradigma de sistema hamiltoniano es el péndulo sin roza- 
miento. Consideremos un péndulo en el que una bola esta su- 
jeta al punto desde el que cuelga por una delgada barra rigida, 
un columpio idealizado. El espacio de las fases consiste en este 
caso en un diagrama donde anotamos en un eje la posicién y en 
el otro la velocidad. La posicién la da el angulo (6) que la bola 
forma con la vertical. 

Cuando la bola se encuentra en reposo en su punto mas bajo, 
la posicion en el diagrama de fases de las variables 0 y v corres- 
ponde al punto en que el Angulo y la velocidad valen cero: 0=0, 
v=0. Si desplazamos ligeramente el péndulo de su posicion de 
equilibrio, la bola ejecuta pequefias oscilaciones. Cuando el des- 
plazamiento es maximo @ tiene un valor determinado por el 
desplazamiento inicial y su velocidad es cero. Luego la velocidad 
va aumentando (en sentido negativo) a medida que el angulo dis- 
minuye. Cuando el Angulo se hace cero, la velocidad es maxima 
(de nuevo en sentido negativo). El 4ngulo cambia entonces de 
signo y la velocidad va disminuyendo (en valor absoluto). Cuan- 
do el péndulo se encuentra de nuevo en su separacion maxima 
pero a la izquierda, la velocidad vuelve a ser cero y de nuevo 
empieza a acelerarse (figura 1). 

En el diagrama de las variables angulo y velocidad (@,v) el 
movimiento de ida y vuelta del péndulo viene representado por 
una linea cerrada, una elipse. Podemos imaginar ahora que re- 
petimos la experiencia con distintos angulos iniciales del pén- 
dulo. Cada movimiento subsiguiente vendra representado por 
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ute lipase diferente, Ee poatble, por lo tante, representar en ol 
diagiama todie ln soluclones de ih eounoidn, que en este ene 
CoOrresponden 8 Carved cermndis, En pequenas oscilactones el 
periodo del pendulo cs siempre el rismo y sleue una ley que 
fue determinada en su dia por Galileo, Pero si desplazanios el 
péndiwo mas y mas cada vez, cl periodo de las oseilaciones se va 
‘largando, Llega un momento en que puede aleanzar el punto en 
el que la bola esta por encima y en la vertical del punto del que 
cuelga. Silo lanzamos con mucha fuerza, el péndulo empezara a 
rotar alrededor del eje, si el punto donde pivota lo permite. De 
un movimiento de oscilacién pasamos a uno de rotacién. Las 
irayecLorias que describen este movimiento en el espacio de las 
fases tienen entonces un aspecto diferente y son, ademas, abier- 
tas. La figura 2 resume estas posibilidades. 

Junto al punto de equilibrio inicial, en el que la velocidad y 
el angulo son iguales a cero, vemos que hay otro punto intere- 
sante en este diagrama en el que la velocidad también es cero 
pero el Angulo es nr. Ambos son puntos de equilibrio del sistema, 
porque si damos justo esos valores a las variables, el péndulo 
se queda ahi, sin moverse. Pero tienen un caracter diferente. El 
punto en que la velocidad y el Angulo tienen valor cero corres- 
ponde al péndulo en la posicién mas baja, y es lo que lamamos 
un eqguilibrio estable. Si desplazamos un poco el péndulo de esa 
posicion, tratara de volver a ella. Mientras tanto, cuando el 4n- 
gulo vale x y la velocidad es cero, lo que corresponde a tener Ja 
bola arriba del todo, el equilibrio es inestable: cualquier pequefia 
perturbacion hace que la bola se aleje de él, cayendo. Ademas 
hay dos tipos de curvas cualitativamente distintas en el diagra- 
ma: las curvas alrededor de @=0 y v=0 son cerradas, pero las 
curvas que representan movimiento de rotacién, cuando supe- 
ramos el punto superior de equilibrio, no lo son (aunque en la 
dinamica real también representan un movimiento periddico). 
En el diagrama los dos puntos de velocidad cero y angulo n o 
—n son diferentes, pero representan la misma situacién fisica: 
el péndulo en su posicién mas alta. La trayectoria que une es- 
tos dos puntos es una trayectoria especial y tiene reservado un 
papel particular en el caos determinista. El tiempo que tardaria 
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Punto mas alto: 
velocidad nula 
NE 
\ 
\ 


Punto mas bajo: —~ 
velocidad maxima 


Representacién en el espacio de tas fases del movimiento de un péndulo. 
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Rotacidn horaria 


> 
=> 
=> 
~»> 
ie 
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Oscilaciones Rotacion antihoraria 


Espacio de las fases del péndulo sin rozamiento. 
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ol péndulo en recerrer esta trayeetoria seria infinite, porque va 
Tiley 78 lenlo sein se aceres Al punto mas alto, Al vequeririun 
Hiempe infinite para desermbinla, eatn iayectoria os irrenlizable 
en la praetica. bs importante darse cuenta de que el periodo de 
tas oscilaciones depende de la. amplitud: crece y crece a medida 
que las oscilaciones son mas amplias. Esta variacién del periodo 
con la amplitud del movimiento es una consecuencia de que la 
écuacion diferencial que describe el movimiente del péndulo sea 
una ecuacion no lineal. De todas formas, el péndulo asi tratado 
é8 un sistema integrable: existe una expresidn matematica que 
predice la velocidad y el 4ngulo como funciones del tiempo, da- 
dos sus valores iniciales. Y ello es asi para todas las trayectorias 
posibles del espacio de las fases. 

Perturbemos ahora esta situacién. Sometamos nuestro pén- 
culo a una fuerza periddica, como si columpiaramos a un nifio 
con un periodo determinado. Y supongamos que el periodo de la 
fuerza es un poco mayor que el de las oscilaciones de pequefia 
amplitud. Teniendo en cuenta que para cada amplitud —cada al- 
{ura maxima del péndulo— el periodo de oscilacion es diferente, 
habré una amplitud para la que el periodo de la fuerza aplicada 
sea resonante, es decir, el periodo natural de oscilacién para esa 
iuuplitud sera el mismo que el de la fuerza aplicada. En esas con- 
(liciones, empujaremos al péndulo de manera sincronizada con 
su movimiento e ira adquiriendo mas y mas energia. El péndulo 
aumentaré su amplitud, pero al hacerlo aumentara también el pe- 
riodo natural de oscilacién, y pronto dejara de tener el mismo 
periodo que la fuerza aplicada, por lo que esta dejara de estar 
sincronizada con el péndulo y dejara de ser resonante. Todo ello 
produce una estructura particular en el espacio de las fases. 

Para revelar mejor los detalles de esa estructura se usa habi- 
tualmente una herramienta matematica conocida como aplica- 
cién o mapa de Poincaré, ideada por el matematico francés, En 
un mapa de Poincaré no se representa una trayectoria completa 
en el espacio de las fases, sino el punto de la trayectoria cada 
cierto tiempo. Para un péndulo sometido a una fuerza periddica, 
el mapa de Poincaré se construye de la siguiente forma: se toma 
un punto inicial cualquiera del espacio de las fases como valor 
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iniclal, se calenia su evalieién a is largo del tiempo, y se van 
mareando en el diagrania lis posltiones que ocupa la Lrayectoria 
6 intervalos de tleripe iguales al perlodo de la fuerza periddica. 
1) proceso se replte para distintos puntos iniciales, para explo- 
far asi todo el espacio de las fuses. 

Fijémonos en la figura 3 en el mapa de Poincaré del péndulo 
(odavia sin perturbar. Hemos tomado como intervalo de tiempo 
para hacer el mapa uno un poco mayor que el tercio del periodo 
de las trayectorias centrales. Se observa que hay curvas cerradas 
que estan completamente rellenas. Estas corresponden a trayec- 
torias cuyo periodo no es multiplo del periodo elegido, y por eso 
la trayectoria pasa cada vez por un punto diferente de la grafica. 
Pero observamos también trayectorias que solo presentan un 
nimero finito de puntos. En particular hay una trayectoria que 


Mapa de Poincaré de un péndulo sin perturbar. 
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solo tlene tres puntos en ete nape de Polnewrd, as la cure 
conlando desde el interior hacio tier, Matos tres puntos corres 
ponden a una trayectora eyo period ea juslo tres veees el que 
hemos elegido para construir el mapa. Al cabo de tres iteracio- 
nes, el punto escogido vuelve al mismo sitio. También es notable 
la trayectoria mas externa, muy cercana a un periodo que es seis 
veces el de referencia, y por eso presenta casi seis puntos. Entre 
ambas hay otra trayectoria con mayor numero de puntos. 

En general, las trayectorias con un periodo exactamente igual 
al que hemos usado de referencia aparecerian como un solo pun- 
to en el mapa de Poincaré, porque al cabo de un periodo se en- 
contrarian de nuevo en el] mismo sitio. Aquellas con un periodo 
igual a dos veces el de referencia aparecerian como dos puntos. 
En definitiva, las trayectorias cuyo periodo sea una fraccién del 
periodo de referencia se convierten en una serie discreta de pun- 
tos en el mapa de Poincaré. Las trayectorias cuyo periodo no 
esta en proporcién racional al de referencia se van rellenando 


’ completamente. El mapa de Poincaré sin perturbar consiste en- 


tonces en una mezcla de curvas casi completamente dibujadas y 
puntos discretos. 

Qué ocurre si aplicamos ahora una fuerza periddica aunque 
su intensidad sea pequefia? Pues sucede lo mismo que con los 
asteroides en el sistema solar. Las trayectorias que estan en re- 
sonancia con la fuerza aplicada se ven fuertemente perturbadas. 
Esta perturbacién da lugar a una curiosa estructura en el mapa 
de Poincaré, como se muestra en la figura 4. La zona en la que se 
encontraba la 6rbita cuyo periodo era tres veces el de la fuerza 
aplicada se ve fuertemente perturbada. Ahi se originan una es- 
pecie de islas. Al mismo tiempo, se ve que hay zonas en las que 
las trayectorias originales casi se mantienen. Estas trayectorias 
estan muy alejadas de resonancias importantes. 


La escuela soviética 


Durante la época comunista, la fisica y las matematicas encon- 
traron en la URSS mas reconocimiento que otras ramas de la 
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Mapa de Poincaré de un péndulo sometido a una fuerza periddica. 


ciencia. Grandes nombres de estas disciplinas estan asociados a 
la Unién Soviética, como es el caso de los fisicos Lev D. Landau 
y Andrei Sajarov, o el matematico Andrei Kolmogorov. Kolmo- 
gorov fue una de las grandes figuras mundiales de las matem4- 
ticas del siglo xx. Trabaj6 en 4reas muy diversas y se interesé 
por todos los grandes problemas de la época. Pero ademas, y 
no menos importante, Kolmogorov fue un gran profesor y bajo 
su direccién se formaron e iniciaron sus investigaciones otros 
grandes matematicos soviéticos. En lo que aqui nos interesa, el 
nombre de Kolmogorov esta asociado a uno de los teoremas mas 
importantes de la teoria del caos determinista en los sistemas 
hamiltonianos: el teorema KAM. 

Como vimos en el primer capitulo, las resonancias entre las 
distintas 6rbitas planetarias pueden desestabilizar las Orbitas 
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(a los planets, Lelie resorincias se producen eacdw vey que las 
Orbitas de dos planetas ienen periodos cuye cocionte es un ri 
mero ryaeional Come ya dijimos, Hen Poincaré demoastro que 
las resonancias hacian imposible encontrar una solucion anali- 
tiga que representara el movimiento planetario en un intervalo 
infinito de tiempo. En el epigrafe anterior hemos visto cémo Jas 
resonancias producen una estructura particular en el espacio de 
las fases del péndulo sometido a una fuerza periddica. 

In 1957 Kolmogorov dio una serie de lecciones sobre diversos 
temas relativos a los sistemas dinémicos. En la segunda parte 
del curso, expuso por primera vez el enunciado de un importan- 
te teorema: dado un sistema haruiltoniano integrable (como el 
péndulo sin perturbar, por ejemplo), hay una infinidad de trayec- 
torias que se mantienen estables ante una pequefia perturbacion. 
In la figura 4 vemos que, aunque hay zonas muy perturbadas por 
las resonancias, hay muchas trayectorias cercanas a sus corres- 
pondientes trayectorias sin perturbar. A medida que se aumen- 
ta la intensidad de la perturbaci6n, mas y mas trayectorias iran 
destruyéndose. 

Kolmogorov no dio una demostracién pormenorizada de su 
teorema, sino que esboz6 solo los pasos principales. La lista de 
problemas que interesaban a Kolmogorov era tan amplia que no 
escribia en detalle todo lo que hacia. Kolmogorov nunca public6é 
una demostracion detallada de su teorema y todo parece indicar 
que hay puntos que no resolvié del todo. Uno de sus colaborado- 
res, Vladimir Arnold, asumi6.la tarea de hacer una demostracién 
pormenorizada del teorema de Kolmogorov. Arnold habfa rea- 
lizado su tesis doctoral bajo la direccién del gran matematico. 
En particular, juntos habfan resuelto el conocido como décimo 
tercer problema de Hilbert, relativo a las raices de la ecuacién 
de séptimo grado. En 1963 Arnold hizo una exposicién detallada 
del problema de la perturbaci6én de un sistema hamiltoniano, de- 
mostrando la existencia de 6rbitas estables para pequefias per- 
turbaciones. La idea central de la demostracién esta relacionada 
con la teoria de numeros, una de las ramas mas dificiles de las 
matemAaticas, y radica en el hecho de que hay 6rbitas en las que 
el cociente entre su periodo y el periodo de la perturbacién es 
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lo sullclontomente eielonal Come pare que TINgun resonanein 
le afecte 0 linge plage, Comoe hemos explicada, in resonancia se 
produce cuando el cociente de los periodos en juego, por ejem- 
plo el del péndulo y el de la fuerza a la que esta sometido, es un 
numero racional. Luego si este cociente es irracional, la trayec- 
toria en cuesti6n nunca sera resonante. El matematico de origen 
aleman Jurgen Moser extendié y completd la demostracion de 
Arnold. El teorema se conoce hoy en dia por las iniciales de es- 
tos tres matematicos: teorema KAM. 

Si Kolmogorov creé una escuela de la que salieron grandes 
matematicos, como el mismo Arnold o Yakov Sinai, Arnold, a su 
vez, influy6 enormemente en toda una generacion de matemati- 
cos rusos. A partir de finales de la década de los sesenta, y hasta 
su muerte, Arnold organizaba todos los martes por la tarde un 
seminario en el que se reunian sus colaboradores. Al principio 
de cada semestre, Arnold se encargaba de presentar una lista de 
problemas interesantes, de todas las areas de las mateméaticas, 
para ser tratados y analizados. Cada uno escogia un problema 
y cada semana se debatian los progresos y las dificultades. De 
esta forma surgieron numerosos teoremas y nuevas teorias. Ar 
nold era un gran deportista y solfa esquiar grandes distancias 
en invierno. El deporte formaba parte de su método de trabajo. 
Seguin su propio relato, si un problema se le resistfa, se ponia los 
esquis y se iba al campo. Al cabo de 40 km vislumbraba, si no la 
solucién, al menos una idea de por dénde seguir. Y la idea solia 
ser correcta. 


UN CARO ERROR 


En la figura 4 vefamos la estructura del mapa de Poincaré de un 
péndulo forzado. Las trayectorias predichas por el teorema KAM 
son claramente distinguibles, como también lo son las islas. Con 
ellas conviven zonas caéticas. Henri Poincaré fue el primero en 
vislumbrar la existencia de estas trayectorias cadticas. La histo- 
ria de su descubrimiento esta unida a un error que le cost6 una 
buena suma de dinero. 
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un ndmero irracional, Cuanto mas grandes sean los nlimercs 4, @,..4,, mejor se aproximara el 
__ ntimero en cuestién por un numero racionall, EI teorema de KAM asegura que al perturbar un 
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_ Que estar suficientemente alejacdas de las resonancias, es decir, su period6 no debe estar en una 


cl6n positiva es: 


pects 
; 2 

sta es la famosa razon durea, la que existe entre la altura y a base del rectangulo dela fi igura, a 

ue considerada jpor los romanos como la proporcién. mas arménica posible. En el.con- 


del teorema.de KAM, aquella orbita cuyo periodo esté en proporcion aurea con ne alge . 
fuerza perurtatora sera la Ultima en desestabilizarse. 


Pero si se trata de un numero inracional (irracional es ith numero gue no’ se puede expresal 
come cociente de dos enteros), la fraccin continua sé extiende hasta et infinito. Cuaiguier nie 
mero irracional se puede aproximar por Su cortespondiente fraccion continua truncada en un 
momento dado. Por ejemplo, el numero x $e aproxima bastante bien por la fraceion: 


Sin ey una representacién nn ble de T requiere infritas fracciones, ya que se trata de 


sistema hamittoniano hay trayectorias que apenas cambian. Para ello, estas trayectorias tienen 


proporcién racional con el periodo de la perturoacion. Pues bien, hay un ntimero irracional que 
esta mas alejado que ningun otro de los numeros racionales. A este numero le corresponde la 
siguiente. fraccion continua: 


En el rectangulo aureo se da la proporcién ats y este cociente es la razon aurea. 
En cada paso, el valor de fos numeros a, a, a, es el menor posible, 1, y por ello ningun ra- arb a 


cional es una buena aproximacién. En la conatruccion de este numero, lo que hay debajo de 
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Uno se quoda impresionada por Ya comentinios en el primer en 
la complajidad de esta figura, que pitulo que el rey Oscar Tl de Suecia 
fo infento siquiera dibular, Nada eon sry 7 cai ii oes dala 

per su OO.” CUIMpVOAnos, Que TUVO Te 
08 mas aproplado pat a dar gar en enero de 1889, [1 premio es- 
una idea de la complejidad del 


taba dotado con una medalla de oro 

problema de los tres cuerpos. y 2600 coronas suecas, una cantidad 

Henai Poncang ada despreciable para la época. Al 

premio optaron varios de los gran- 

des matematicos de entonces y Henri 

Poincaré presenté a concurso una memoria de cerca de trescien- 

tas paginas sobre el problema de los tres cuerpos. El fondo del 

asunto era la estabilidad del sistema solar como un todo. Para 

ilustrar sus nuevos métodos y herramientas, Poincaré se centré 
en una version restringida del problema de los tres cuerpos. 

La memoria de Poincaré era de tal complejidad y extension 
que ni siquiera el jurado, formado por grandes mateméaticos 
como Mittag-Leffier, Charles Hermite y Karl Weierstrass, la 
pudo entender en todos sus detalles en una primera lectura. El 
premio llevaba aparejada la publicacién del ensayo ganador y 
Mittag-Leffler encarg6 a un joven matematico sueco, Lars Ed- 
vard Phragmén, la preparacion del trabajo de Poincaré para su 
edicién. Poincaré y Phragmén se cartearon con regularidad para 
depurar la memoria y aclarar los muchos puntos oscuros, Cuan- 
do ya estaba todo casi listo para su publicacion, Phragmén des- 
cubrié un razonamiento que no le parecia completo. Tenia que 


laquierda: trayectoria homoclina sin perturbar. Derecha: enredo homoclino. 


al punto de equilibrio la velocidad tiende a cero, y lo mismo 
ocurre cuando salimos de él. En una parte de su memoria, Poin- 
caré se interesaba por la estabilidad de las trayectorias homo- 
clinas en el contexto del movimiento planetario: ise mantenian 
si se consideraba la perturbacién creada, por ejemplo, por un 
planeta pequefio? 

Poincaré habia considerado previamente las trayectorias que 
entran o salen de un punto de equilibrio. En una homoclina estas 
dos curvas son la misma. En un principio, el matematico fran- 
cés crey6 demostrar que las dos curvas seguian siendo la misma 
ante pequefias perturbaciones: la homoclina se mantenia como 


ver con la estabilidad de las érbitas homoclinas. 
Las 6rbitas o trayectorias homoclinas son trayectorias espe- 
ciales que unen un punto de equilibrio consigo mismo. La figura 


5a muestra un ejemplo. 


Si la trayectoria une dos puntos de equilibrio distintos se de- 
nomina heteroclina. En la figura 2 (pag. 117), en la que hemos 
representado el espacio de las fases de un péndulo sin roza- 
miento, aparecen dos trayectorias heteroclinas, que unen los 
equilibrios en los que la velocidad es cero y el angulo x 0 —1. 
Estas trayectorias son especiales porque se tarda un tiempo in- 
finito en recorrerlas. En efecto, a medida que nos acercamos 
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tal. Pero un comentario de Phragmén le hizo darse cuenta de 
que estaba equivocado. No tard6é en encontrar la verdadera res- 
puesta: las dos curvas no coinciden, sino que se cortan infinidad 
de veces. La figura 5b muestra un esquema de algunos de es- 
tos primeros cortes. La parte de la homoclina perturbada que se 
acerca al punto de equilibrio se retuerce cortando una y otra vez 
a la otra curva que llega al equilibrio. Los cortes se producen a 
distancias cada vez mas cortas y la curva se alarga cada vez mas 
en uno y otro sentido, retorciéndose. Por otro lado, lo mismo 
ocurre con la otra rama, trazada con puntos discontinuos. Esta 
curva tiene un aspecto similar, es la imagen invertida de la curva 
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de eave continue, Tambien ve corn iafinihis yeeos y, pir hacer 
lo, se Glairga ani y otf vex enn sentido y Giro, 

Los Infinitos Cortes sucesivos de las dos ramas en que se ha div 
vidideo la antipua homocling crean una marafa que estira y estirn 
puntos proximos, plepandolos luego una y otra vez. Ein palabras 


FG. 6 


Mapa de Poincaré de un 
péndulo forzado. En él 
se ven zonas cadticas. 
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de Poineareé: 


|...| esta figura esta formada por dos curvas y sus intersec- 
ciones en nimero infinito [...] Estas intersecciones crean 
una especie de enrejado, de tela o malla infinitamente tupida; 
estas curvas no se cortan nunca consigo mismas, pero deben 
replegarse sobre si mismas de manera muy compleja, para 
volver a cortar una infinidad de veces las fibras de la malla. 


Dentro de un enredo homoclino como el de la figura 5b, los 
puntos sucesivos de una trayectoria se van estirando y plegan- 
do. Si, el lector ya lo habra adivinado, estos son los elementos 
esenciales del caos: estirar y plegar. E] enredo homoclino es al 
caos hamiltoniano lo que el atractor extrafio al caos en sistemas 
disipativos: la estructura geométrica subyacente. 

En la figura 6 mostramos el mapa de Poincaré del péndulo 
perturbado en toda su complejidad: hay islas, curvas KAM (muy 
poco perturbadas) y zonas ca6ticas, correspondientes a las he- 
teroclinas perturbadas. La zona cadtica es mas notable en la 
parte exterior y se manifiesta en puntos dispersos, Estos puntos 
corresponden a una trayectoria que vaga entre las cercanias de 
los puntos de equilibrio inestable iniciales, donde la velocidad 
es cero y el angulo z 0 -n. Esta figura también presenta otra de 
las caracteristicas del caos, la autosimilaridad. Las islas de las 
trayectorias resonantes son similares al mapa sin perturbar ini- 
cial, y en su interior se desarrollan otras islas y aparecen nuevas 
trayectorias heteroclinas. Al ir aumentando la intensidad de la 
perturbacion, estas heteroclinas se desestabilizan a su vez, dan- 
do lugar a mas caos. 

Cuando Poincaré comunicé por carta a Mittag-Leffler que ha- 
bia descubierto un error en la memoria original presentada para 
el premio del rey Oscar I, y que queria corregirlo antes de que se 
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publicara, ya era demasiado tarde. El articulo ya estaba impreso 
y habia empezado a distribuirse. Mittag-Leffler par6 la distribu- 
cidén, pero pidié a Poincaré que se hiciera cargo de los costes de 
reimpresion del articulo corregido. El coste total de la tirada era 
de 3500 coronas, mil mas de lo que Poincaré habia recibido por 
el premio. 


CAOS CUANTICO 


La mecAanica cuantica es la teoria fisica del mundo microscépico: 
las moléculas, los 4tomos y el nticleo atémico. Fue desarrollada 
en las primeras décadas del siglo xx por una generaci6n irrepeti- 
ble de fisicos. Una lista no exhaustiva incluye los nombres, ya le- 
gendarios, de Planck, Einstein, Bohr, Heisenberg y Schrédinger. 
La fisica cuantica se asienta sobre una serie de conceptos que 
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Ponipen COn Ta ecinion newtoniind Entre ellos eatin ta hipdte: 
ala cuiintion, In dunlidad ondmcarpiseuls, el principio de ineeri 
dumbre y el principle de correspondencia, Como necesliamos de 
ellos para entender el resto del capitulo, conviene que recorde 
mos, aunque sea brevernente, en qué consisten estos principios, 

La hipétesis cudntica fue formulada por el fisico aleman 
Max Planck a finales del aio 1900. Segtin esta, la energia de 
un sistema mecanico oscilante solo puede tener un conjunto 
disereto de valores, multiplos de una cantidad m{nima, el cuan- 
to, de valor F=hv, donde v es la frecuencia de oscilacién, y 
h, una constante, la constante de Planck. Esta serie de valo- 
res recibe el nombre de niveles de energta del sistema. En su 
formulacion inicial, la hipé6tesis de Planck se aplicaba a siste- 
mas oscilantes con una frecuencia fija: un muelle, por ejemplo. 
En su formulacidén final, la mecdnica cuantica proporciona las 
herramientas matemAaticas para calcular los niveles de energia 
permitidos en cualquier sistema mec4nico. En el caso de siste- 
mas ligados, por ejemplo, un electrén y un protén que forman 
un 4tomo de hidrégeno, estos niveles de energia permitidos 
forman un conjunto discreto. Heisenberg y Schrédinger, entre 
1924 y 1926, propusieron dos formulaciones diferentes, pero 
equivalentes, de la mecanica cudntica. Ambas formulaciones 
permiten, en principio, calcular los niveles de energia de un 
sistema mecanico cualquiera. 

En la figura 2 habiamos dibujado el espacio de las fases de 
un péndulo. Las trayectorias cerradas se corresponden con un 
movimiento oscilante, de vaivén, del péndulo. A cada una de es- 
tas trayectorias le corresponde una energia mecanica definida, 
suma de la energia cinética y la energia potencial, y esta energia 
es tanto mayor cuanto mas amplias son las oscilaciones. En la 
fisica clasica todas las energias son posibles: los valores de ener- 
gia permitidos forman un continuo. La hipotesis cudéntica implica 
que, en un sistema cudntico, no todas las energias son posibles, 
solo lo son un conjunto discreto de ellas. Antes de que Heisen- 
berg y Schrédinger hicieran sus respectivas contribuciones, el 
aleman Arnold Sommerfeld dio una regla de cuantificaci6n que 
permitia calcular los niveles de energia de ciertos sistemas sim- 


CAOS CLASICO Y CAOS CUANTICO 


plow, cane el péndtils, La regia de Sommerfeld suponia que dn 
Jongitud total de las trayeetoriia admitidia en 6] 6gpacio de lus 
fuses debia ser un mlliplo de la conslante de Planck, Bata regla 
permitia calewar con éxito los niveles de energia del oscilador 
arrndénico de Planck y del Atomo de hidrégeno, por ejemplo. 

La dualidad onda-corpusculo fue intuida por Einstein, en pri- 
mer lugar, y explotada en profundidad por Louis de Broglie y 
frwin Schrédinger. Seguin De Broglie, toda particula cuantica 
lleva asociada una onda cuya longitud de onda es: 


La masa de la particula esta representada por ™, y v es su velo- 
cidad. El producto mv es la cantidad de movimiento, p, de la par- 
ticula. En la formulacién de Schrodinger de la mecanica cuantica, 
conocida como mecdnica ondulatoria, toda particula leva aso- 
ciada una funci6n ‘¥(x, y, 2), la funcién de onda, cuyo cuadrado 
da la probabilidad de encontrar a la particula en un punto deter- 
minado del espacio. Schrédinger dedujo una ecuacion, que leva 
su nombre, cuya resolucién proporciona la funcion ‘P(x, y, 2). En 
la figura 7 hemos representado el cuadrado de la funcién de onda 
para una particula sometida a una fuerza andloga a la de un muelle 
cldsico en el estado mAs bajo de energia. Desde una perspectiva 
clasica, la particula ejecutarfa un movimiento oscilatorio alrede- 
dor del punto de equilibrio (marcado en la figura con la linea ver- 
tical discontinua). Desde una perspectiva cuantica, la funcién de 
onda nos da la probabilidad de encontrar la particula en una zona 
del espacio. La funcioén es maxima en el punto de equilibrio, lo que 
significa que es ahi donde es mas probable encontrar la particula. 
A medida que nos alejamos del equilibrio la probabilidad disminu- 
ye, y es practicamente cero para puntos alejados. 

El comportamiento de una pequefia particula, como un elec- 
trén en un atomo, se parece mas al de una onda en la superficie 
de un lago que a la traza que deja un boligrafo sobre una hoja de 
papel. Este tipo de movimiento ondulatorio es el que describe 
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bubien fenGrienios core th propagicion de in tie o el sonido 
en tode ipo de medion, Matis ondas son la solver de eeineic 
Nee eA enente sionitires ala de SehrOdinger, euyas so 
lieiones suelen lener in dapecto suave y continuo, Las particu. 
laa cufinticns comparten asi propiedades Mindamentales de las 
ondas clisicas, fn partieular, cuando dos ondas se superponen 
se producen interforencias; los maximos se refuerzan y los mi 
mmos y maximos se cancelan, El fenédmeno se puede observar, 
por ejemplo, si tiramos simullaneamente dos piedras a un lago; 
en la zona donde las dos ondas coinciden se produce un patron 
de interferencias caracteristico. La interferencia de las ondas de 
materia es uno de los fen6menos cuanticos mas notables y sera 
um punto importante en nuestra discusién del caos cuantico. 

Un atomo tiene ciertas similitudes con un sistema planetario. 
In el Alomo de hidrégeno un electrén esta sometido a la fuerza 
éléectrica atractiva del nicleo, que consiste en un solo protén. 
La fuerza eléctrica, que sigue la ley de Coulomb, tiene la mis- 
tna. forma que la atraccién gravitatoria: su intensidad varia como 
el inverso de la distancia al cuadrado. En el] caso de un cuerpo 
celeste que orbite alrededor del Sol se pueden distinguir dos ti- 
pos de trayectorias. Por un lado estan las trayectorias cerradas, 
como las de los planetas. Estas ocupan una region determinada 
del espacio, de la que el planeta no sale nunca, y son 6rbitas 
elipticas. Pero existen otras trayectorias posibles en las que el 
cuerpo viene de muy lejos, pasa cerca del Sol solo una vez y 
luego se aleja indefinidamente. Estas trayectorias son hipérbolas 
y las describen algunos cometas. En un diagrama en el espa- 
cio de las fases, las érbitas cerradas tienen un aspecto similar a 
las trayectorias interiores de la figura 2, mientras que las 6rbitas 
abiertas, hiperbdlicas, presentan un aspecto semejante al de las 
trayectorias externas a las heteroclinas de la misma figura. 

En el caso del atomo de hidrdégeno, el electrén también puede 
hallarse en dos tipos de estados. Hay estados ligados, los que 
corresponderian a las orbitas elipticas clasicas, y estados no li- 
gados, en los que el electrén no esta preso y puede abandonar el 
atomo. Para que los estados ligados sean estacionarios, la fun- 
cion de onda de Schrédinger debe ser una onda estacionaria, Si- 
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Poncl6n de onda para la 

\ version cuantica de una 

! particula sujeta a un muelle. 
Se ha representado la funcion 

correspondiente al nivel mas 

j bajo de energia. La linea 

| discontinua representa 

\ el punto de equilibrio. 


milar alas ondas que se producen en las cuerdas de una guitarra. 
En una onda viajera los méximos y minimos van cambiando de 
posicion a medida que la onda se propaga. En una onda estacio- 
naria los nodos (puntos donde la onda se anula) y los vientres 
(puntos donde las oscilaciones alcanzan mayor amplitud) ocu- 
pan siempre la misma posici6n. En los sistemas ligados ello im- 
plica ciertas restricciones sobre la forma de la onda, que solo se 
cumplen para algunos valores de la energia. El conjunto de esta- 
dos ligados es, por lo tanto, un conjunto discreto, al que corres- 
ponde un conjunto discreto de niveles de energia. Los estados 
no ligados vienen representados por ondas viajeras, que pueden 
abandonar el atomo, y no tienen ninguna restriccion, por lo que 
los valores posibles de la energia forman un continuo. 

El principio de incertidumbre o de indeterminacion fue for- 
mulado por Heisenberg en 1927, afio en el que se encontraba 
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en Copentiague trabajando estrechanente con Niels Bohr E1 
poncipio de indetermiiacion se enunela met ln precision Av econ 
la que se puede medir Ja posicion ¢ de una particula y la preci- 
sién Ap con la qué se puede medir su cantidad de movimiento 
y no son independientes, sino que deben cumplir la siguiente 
desigualdad: 


h 
ApAx 2— - 
P An 


La consecuencia de ello es que no podemos determinar a la vez 
la posicién y la cantidad de movimiento de una particula con ab- 
soluta precision. Si conocemos exactamente su posicién no hay 
incertidumbre en las coordenadas espaciales de la particula y Ar 
sera cero. La desigualdad obliga entonces a que Ap sea infinita, 
por lo que nuestra indeterminaci6n sobre la cantidad de movi- 
miento sera total y esta podra tomar cualquier valor. Nada sabre- 
mos sobre la cantidad de movimiento. Igualmente, si conocemos 
con total precisién la cantidad de movimiento de una particula, 
Ap sera cero y Ax infinito, con lo que la particula puede estar en 
cualquier lugar del espacio. Como la cantidad de movimiento de 
una particula es igual a la masa por su velocidad, lo que decimos 
de la primera se aplica también a la segunda: no podemos conocer 
simulténeamente con toda precision la posici6n y la velocidad de 
una partfcula. Con ello, la mecAnica cudntica despoja de sentido 
al concepto de trayectoria de un cuerpo, ya que si medimos la po- 
sicién de una particula en un instante dado, al no saber cual es su 
velocidad, seré imposible determinar dénde se encontrara en el 
instante siguiente. De nuevo Ja figura 2 nos sirve de ejemplo para 
ilustrar este concepto. En fisica cuantica el diagrama de fases de 
la figura 2 no tiene sentido, ya que en él representamos, en cada 
instante de tiempo, la posicién y la velocidad del péndulo. 

En este rapido repaso a la fisica cudntica, nos queda hablar 
del principio de correspondencia. Enunciado por Niels Bohr 
en torno a 1918, sirvié de guia en la elaboracién de la mecanica 
cuantica. Seguin este principio, los resultados obtenidos por la 
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mecinion cuintion deben eainedir con low de ln fiview clisien 
Si hacemos tender © cere la constinte de Planek. Ora forma de 
verlo, 8 nos fijmnos en In relneian entre longitud de onda y can- 
lidad de movimiento de De Broglie, es decir que cuando 1a longi- 
tud de onda de la onda asociada sea pequernia, Compurada con las 
dimensiones fisicas del sistema, la particula debe comportarse 
clasicamente, ya que los efectos ondulatorios no seran notables. 

El principio de correspondencia recuerda a la relacién que 
existe en 6ptica entre las ondas de luz y los rayos asociados. 
La luz es un fenédmeno ondulatorio y su propagacién viene 
descrita por frentes de onda que marcan los maximos o mini- 
mos de la onda. A una sucesién de frentes de onda le podemos 
asociar un Conjunto de rayos si trazamos lineas perpendicu- 
lares a los frentes en cada punto. En la figura 8 se muestra un 


FG. 8 


Frentes de onda 

y rayos asociados. 

Las circunferencias 
concéntricas representan 
lugares donde la onda 
tiene un maximo y estan 
separadas por una longitud 


de onda. 
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dJemplo en dl que lenemos unin onda que parte de un punta: 
low Trentes de onda son ealéreos y los Tayos Asoclados parten 
Tadivimente, Cuando no hay fenémenos de interferencia, loa 
rayos proporcionan tune informacion completa de la direecion 
dé propagacién de la onda y basia con trazar los rayos para 
describlr los fenémenos dopticos, Nsto es lo que hace la éptica 
peométrica, que es la base del disefho de diversos aparatos ép- 
ticos como prismaticos, telescopios, etcétera. Pero si las on- 
das interfieren, los rayos no proporcionan toda la informacién, 
es necesario saber si la onda se encuentra en un m4ximo o en 
un minimo para conocer el resultado de la interferencia: la 
descripcién de la onda completa, gracias a los frentes de onda, 
es ineludible. 

La funcién de ondas de Schrédinger contiene toda la infor- 
maci6n que la mecdnica cudntica puede darnos de una particu- 
la. Cuando el movimiento de la particula es tal que su longitud 
de onda es pequefia comparada con las dimensiones del siste- 
ma del que forma parte, basta trazar los rayos perpendiculares 
a los frentes de onda para tener una descripcién completa del 
movimiento, Estos rayos representan la trayectoria clasica de 
la particula y, segiin el principio de correspondencia, deben ser 
calculables a partir de las leyes de la mecadnica clasica. 

En lo que al caos determinista se refiere, no cabe esperar en 
el mundo cuantico el comportamiento que hemos encontrado 
en los sistemas cldsicos. La raz6én esta en el principio de incer- 
tidumbre. La nocidn de trayectoria de una particula carece de 
sentido en el mundo cuantico y, por tanto, tampoco se puede 
dar la sensibilidad a las condiciones iniciales, una de las se- 
fias de identidad del caos determinista. ;Cémo se reconcilia el 
aspecto extremadamente irregular de las trayectorias cadticas 
clasicas con las caracteristicas suaves de las ondas propias de 
la escala atémica? ,Existe el caos determinista en el mundo 
cuantico? Los sistemas cuénticos no presentan un comporta- 
miento cadtico al modo clasico, pero se usa la expresién caos 
cudntico para designar las caracteristicas de aquellos sistemas 
cuanticos cuyo sistema clasico asociado, aquel que les asigna- 
rfa el principio de correspondencia, es cadtico, Algunos, como 
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1 finies britinies Michael Berry, prefleren usar eaologia odin 
Hea en lugar de ¢aos cndnlice, 


El billar cuantico 


Un billar es un buen ejemplo para ilustrar las diferencias entre 
los comportamientos clésico y cuantico de un sistema cadtico. 
Imaginemos una particula obligada a moverse en una superficie 
bidimensional cerrada, como una bola en un billar. Pero, en lu- 
gar de un billar rectangular, empecemos por un billar un poco 
raro, uno cuyas paredes sean una circunferencia, un billar circu- 
lar. Supondremos también que nuestra bola de billar no sufre ro- 
zamiento con la superficie por la que se mueve y que las paredes 
son totalmente rigidas, por lo que no se pierde energia en cada 
choque. Se trata, pues, de un sistema hamiltoniano, en el que la 
energia de la bola se conserva. 

Consideremos primero que la bola se mueve segun las leyes de 
la mecanica cldsica. Una vez que le demos impulso, rebotara una 
y otra vez con las paredes sin perder nunca la energia inicial que 
le hemos dado. Cada vez que la bola choca con la pared se refleja 
siguiendo las leyes de la mecanica de Newton, cuya aplicacién, 
en este caso, conduce a una regla simple: el Angulo de incidencia 
es igual al de reflexion. La figura 9 muestra una trayectoria tipica 
de una bola en este billar circular. La trayectoria no se repite 
nunca, pero no es caédtica, se trata de lo que Ilamamos una drbita 
cuasi periddica. Como el contorno es una circunferencia, cada 
vez que la particula rebota «siente» una fuerza que esta dirigida 
hacia el centro de la circunferencia, una fuerza central. Por ello, 
ademas de la energia también se conserva el momento angular, y 
la trayectoria esta perfectamente determinada conocidos la po- 
sicién y la velocidad iniciales. Se trata de un sistema integrable. 
Ademas, hay toda una zona, un circulo interior, que la particula 
no visita nunca. Este sistema no presenta sensibilidad a las con- 
diciones iniciales, como muestra la figura 10, en la que hemos 
representado dos trayectorias que parten del mismo punto pero 
con 4ngulos ligeramente diferentes. Aunque las trayectorias se 
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Ven separnnide, iran cuatro (rielones aljuen relalivaamente cer 
va hun de la obra, 

Sein la teorie eudnties, Ia eeuacion de ondas de Schrédinger 
proporciona los posibles estados cuanticos. Como la bola esta 
obligada a moverse en una zona fintta del espacio, se trata de un 
sistema ligado y el conjunto de estados cuanticos es discreto. 
1] cuadrado de la funcién de ondas representa la probabilidad 
de encontrar la particula en un punto dado del espacio. A cada 
funcién de onda y, por tanto, a cada estado cuantico, le corres- 
ponde un nivel de energia determinado. La figura 11 muestra el 
cuadrado de la funcion de onda del billar circular de Ja figura 9. 
Se muestra la funcidén que corresponde al nivel de energia nime- 
ro 95. Hemos elegido un nimero tan alto para estar en condicio- 
nes en las que el principio de correspondencia sea aplicable, ya 
que a los niveles mas altos les corresponden longitudes de onda 
cada vez mas pequefias. Vemos que la trayectoria de la figura 9 se 
corresponde con la probabilidad que refieja la figura 11: la zona 
cuanticamente poco probable esta en la zona en la que clasica- 
mente no hay trayectorias, y la zona mas probable cuanticamen- 
te coincide con la circunferencia interior, que es tangente a un 
gran numero de trayectorias clasicas. 

Cambiemos ahora de billar. En la figura 12 tenemos un billar 
en forma de estadio. En él hemos representado una trayectoria 
lo bastante larga como para vislumbrar las caracteristicas ge- 
nerales del movimiento. El estadio no tiene simetria central y 
en é] no se conserva el momento angular. Es un sistema no in- 
tegrable: no es posible dar una expresié6n matematica general 
para la trayectoria que sigue la bola conociendo su posicidén y 
su velocidad iniciales. E] aspecto de la trayectoria es ahora mu- 
cho mas irregular. En realidad, si siguiéramos la trayectoria de 
la bola durante un gran nimero de rebotes, veriamos cémo las 
sucesivas trazas van llenando toda el area disponible, hasta re- 
llenarla completamente. Ahora si hay sensibilidad a las condicio- 
nes iniciales. Para comprobarlo, basta fijarse en dos trayectorias 
como las representadas en la figura 13. Tras cinco reflexiones en 
la pared, pasan por los puntos 1 y 2, en extremos opuestos del 
billar. Esta situacion difiere claramente de la del billar circular, 


Ln qrifien dle ta figura 9 muestra la trayectoria clasica de una particula en un billar circular. El punto marca 
In ponlelon Iniclal. En ta flgura 10 se han representado dos trayectorias que parten de la misma posicién, 
Hareuda Lon UN punto, pero con angulos tigeramente distintos. Tras cuatro reflexiones se llega a los 
fuilos marcados con los ndmeros 1 y 2. 


maou 


Tundtado de la funcion de onda para el nivel de energia numero 95 del bitlar circular. Las zonas mds oscuras 
Jopreeentan poca probabilidad de encontrar a la particula, las zonas claras son las mds probables. 
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como se aprecié al compararta con las fyuras Py 10, 1 pillar de 
ostadio en un sister eHdtleo, 

in Ja figura 14 ae representa ¢] cuadrado de la funcién de 
onda para él misma billar dé estadio. De nuevo hemos escogido 
el nivel 95 de energia. Sorprendentemente, la distribucién de 
probabilidad de hallar la particula no es uniforme, como seria 
de esperar considerando el hecho de que las trayectorias cla- 
sicas recorren con igual probabilidad toda el area del estadio. 
Esta circunstancia se explica por el papel que tienen las inter- 
ferencias entre ondas. En el billar de estadio clasico hay algu- 
nas trayectorias periddicas. Por ejemplo, si lanzamos la bola a 
lo largo de la linea central horizontal rebotara una y otra vez 
sin salir de la horizontal. También es periddica una orbita que 
atraviese el centro del semicirculo izquierdo e incida justo en 
el centro del lado recto superior, ya que al rebotar pasara por el 


Cuadrado de la funcién de onda para el nivel de energia numero 95 del billar con forma 
de estadio. Las trayectorias periddicas subyacentes crean un patron de interferencia que 
induce cierto orden en Ja cistribucién, 
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Niveles de energia del billar circular. 

Se han representado 40 niveles, entre el 81 

y el 120. La diferencia entre niveles, AE, llega 
a tomar valores muy pequefios. 
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Niveles de energia del billar con forma de 
estadio. Se han representado también los 


Estos tienden a distribuirse de manera mas 
uniforme que en el caso del billar circular. 


40 niveles comprendidos entre el 81 y el 120. 


felaciona con las simetrias. En el billar circular o en el 
je simetria de revolucidn, porque es {a misma para todos los 
a del centro del sistema. Ello hace que varias funciones de 
ismo valor ‘de. la -energia. En el argot ouéntico: 86 dice ¢ que son See 


ido iN campo exietno o en una molécula no simétrica, se rorrie la degeneracion yl fos’ nat 
a3 5 eee ai Los sistemas ia cuyo sistema clasico ascoiedo.. : 


Eugene Paul Wigner recibié el premio Nobel de Fisica en 1963. En la fotografia aparece contemplando ef diploma 
acreditativo de tan preciado galardon. 
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centro del semletreulo derecho y se reflajard por la miaraa lines 
que Hevaba, Ambas posibilidades catdn trazadas en la figura 165, 

Las Srbitas perlédicas en é) billar de estadio clisieo son muy 
excepcionales. De hecho, si lanzamos la bola al azar, la probabi- 
lidad de que siga una trayectoria periddica es cero. Sin embargo, 
las é6rbitas periédicas son muy relevantes en el billar cudntico. 
Ello se debe a que las ondas asociadas a las érbitas caéticas in- 
terfieren de manera incoherente. Las multiples reflexiones hacen 
que maéximos y minimos interfieran aleatoriamente, no dejando 
un patrén claro. Por el contrario, las ondas asociadas a las 6rbitas 
peridédicas interfieren siempre de la misma forma, con coheren- 
cia, y producen un patrén definido que perdura. E] billar cuantico 
es, asi, «menos cadtico» que su andlogo clasico. Este papel de- 
terminante de las érbitas periddicas en el comportamiento de los 
sistemas cuanticos fue desvelado por Martin Gutzwiller. 

Los billares cudnticos pueden parecer sistemas demasiado 
idealizados como para tener interés practico, pero lo cierto es que 


‘ hoy en dia se han conseguido realizar en el laboratorio gracias al 


avance de las técnicas de construcci6n de microdispositivos. 

Es posible construir cavidades microscépicas superficiales 
donde encerrar la luz. Estos dispositivos se asemejan a laseres 
bidimensionales y las cavidades se pueden construir con formas 


Dos trayectorlas periédicas en el billar con forma de estadio. 
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muy diversas, Las cavidades con for Pequefas maléculas como el 


ma de estado, de dimensiones del acetileno tienen decenas de miies 
orden de 40 micras de largo por 20 de de ggtados cudnticos para nuestro 


ancho, presentan espectros de emi- 
sién que estén directamente relacio- 
nados con los estados posibles, simi- 
lares a los calculados para un billar 
cuadntico con Ja misma forma. Se comprueba que las trayectorias 
periddicas del billar determinan el espectro, El que se trate de 
luz y no de particulas es irrelevante, porque la forma de las on- 
das estacionarias dentro de la cavidad es la misma en ambos 
casos. También es posible construir trampas para electrones con 
forma de estadio y los niveles de energia se distribuyen de la 
manera esperada. 

En un 4tomo de hidrégeno, la fuerza a la que esta sometido el 
electrén depende solo de la distancia al nucleo. Es una situacién 
altamente simétrica y el sistema clasico asociado es integrable y 
tiene una solucién conocida: las 6rbitas de Kepler. Por lo tanto, 
se trata de un sistema no cadtico. Pero si sometemos el 4tomo 
de hidrégeno a un campo eléctrico 0 a un campo magnético ex- 
terno, la situacidn se asemeja a la del péndulo perturbado y el 
sistema cldsico correspondiente se vuelve cadtico. Estos siste- 
mas han sido estudiados experimentalmente, y presentan carac- 
teristicas especiales en la distribucidn de los niveles de energia, 
similares ala de los billares cuanticos cadticos. 
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entretenimiento y edificacién. 
Martin GutzwiLLer 


Desde la aparicion de los primeros estudios sobre 
el caos determinista, la teoria del caos ha atraido 
la atencién de cientificos de todas las disciplinas. 
4Cdémo ha influido esta teoria en areas tan 
dispares como la ingenieria, la biologia, 

la economia, la quimica e, incluso, la filosofia? 


Si nos asomamos a contemplar un rio desde un puente ubicado 
en su curso alto, cuando la corriente baja con fuerza observa- 
remos que se forman remolinos en la estela de los pilares del 
puente. El] movimiento de los remolinos tiene el mismo poder 
hipnotico que el fuego en una chimenea: es siempre el mismo y 
siempre distinto. La estela tiene una estructura a gran escala que 
se abre a medida que la corriente se aleja del puente. Dentro de 
la estela se produce una sucesién de torbellinos, més o menos 
regular dependiendo de la velocidad del agua. Si el rio baja muy 
rapido, las estelas de los pilares llegan a ser muy turbulentas 
y los detalles nunca se repiten. Se pueden observar entonces 
torbellinos de todos los tamafos: desde algunos que tienen las 
dimensiones del pilar del puente hasta otros muy pequefios. Se 
forman debido a la presencia del puente y a la viscosidad del 
agua. La viscosidad es la resistencia que presentan los fluidos al 
movimiento. Cuando el agua pasa por debajo del puente, los pi- 
lares frenan la que pasa pegada a ellos, mientras que el agua que 
fluye en el espacio medio entre dos pilares lo hace a gran veloci- 
dad. La viscosidad hace que unas capas de fluido vayan frenando 
alas capas vecinas, transmitiendo la presencia del puente a todo 
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ef jG, 21 frenads viseouo induce la formaelén de low vortices, 
del mismo mode que ln frieeiGn con el suelo hace fiver une pelo- 
la Cuando In lanzninos, o que wn fulbollate en earrera cae pirando 
sobre at rmiamo Cuando un defensa lo agarra del brazo. Las carac- 
loristicas generales del movimiento turbulento son la formacion 
de torbellinos © vortices y la irregularidad del movimni ento, tanto 
espacial como temporal. 

La turbulencia aparece en rultitud de situactones: en el humo 
que asciende violentamente por la chimenea de la fabricas, en 
la estela que dejan en el aire los automéviles, en los chorros de 
agua, que salen de la depuradora de una piscina, etcétera. El que 
un flujo sea o no turbulento no depende solo de la velocidad 
del agua. Hacia 1880, el fisico inglés Osborne Reynolds estudié 
e] flujo en un tubo cilfndrico de gran longitud. Reynolds descu- 
brié que el hecho de que el flujo fuera o no turbulento dependia 
no solo de la velocidad del liquido en cuestidn, sino también de 
su viscosidad y del didmetro de la tuberia. También comprobé 
que, en lo que se refiere a la aparicién de la turbulencia, duplicar 
Ja velocidad tenia el mismo efecto que duplicar el didmetro del 
tubo o que dividir por dos la viscosidad del I{quido. Todo esto 
queda condensado en un ntimero que hoy en dia llamamos niti- 
mero de Reynolds, Re, que viene dado por el siguiente cociente: 


donde p es la densidad del liquido, » su viscosidad, U la velo- 
cidad del fiujo y d una distancia caracterfstica (el didmetro del 
tubo en el caso estudiado por Reynolds, el didmetro del pilar en 
el caso de un puente). Flujos que tienen el mismo numero de 
Reynolds son similares, aunque los liquidos 0 las escalas sean 
distintos. Si el ntimero de Reynolds es bajo, los flujos son siem- 
pre laminares, esto es, el fluido fluye con regularidad, como un 
rio cerca de la desembocadura. Si el ntimero de Reynolds es su- 
ficientemente alto, el flujo es turbulento, tanto mas cuanto més 
alto sea el ntimero de Reynolds. 
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A MEZCLAN 


Una propiedad esencin de ln turbulencia desde el punto de vis- 
ta Industrial 68 su capacidad de mezcla. En un flujo turbulen- 
to las distintas partes del fluido se mezclan eficazmente. Este 
hecho lo comprobamos, por ejemplo, cuando agitamos el café 
con una cucharilla para disolver el azcar. Si echéramos el azu- 
car en el café y no lo agitaramos, tardarfa mucho en disolverse. 
La raz6n es que la disolucién es un proceso que tiene lugar a 
escala molecular: las moléculas de azucar se difunden entre 
las moléculas de agua. La escala de tiempo a la que se da este 
proceso es proporcional al cuadrado del tamario del sistema, 
por ejemplo, la altura de la taza, e inversamente proporcional 
al coeficiente de difusién. Esto se expresa matematicamente 
como d@?/D, donde d es la escala espacial, y D, el coeficiente 
de difusion. El coeficiente de difusi6n da idea de la facilidad 
que tienen las moléculas de la sustancia que se difunde de pe- 
netrar entre las moléculas de la sustancia en que se introduce, 
y suele ser un valor muy pequeno. En el caso del azticar en el 
café, D~10* m*/s, por lo que si tomamos la férmula anterior 
para d~5 cm, obtenemos un valor cercano a una hora. En con- 
clusién, la difusién por si sola es un proceso muy lento para 
obtener una mezcla homogénea. ,Cémo lo remediamos? Remo- 
viendo el café con la cucharilla. Al hacerlo generamos un flujo 
turbulento que pone en contacto el azticar con todo el café de 
la taza y los mezcla més eficazmente. 

El matematico ruso Andrei Kolmogorov propuso, en los afios 
cuarenta del siglo pasado, una teoria de la turbulencia segiin la 
cual la energia asociada al movimiento turbillonario pasaba de 
los torbellinos de gran escala a los mas pequefios. E! mecanismo 
fisico por el que se produce esta cascada de Kolmogorov, como 
se le conoce, es el estiramiento que el propio flujo produce en los 
torbellinos grandes. Un torbellino es una estructura tridimensio- 
nal, como la que observamos en el desagiie del fregadero. Al ser 
arrastrado por el movimiento del fluido en un flujo turbulento, 
el torbellino se estira en la direccién de su eje y se encoge en el 
plano perpendicular a él. Al encogerse se va acelerando, como un 
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patinador eens plhoge low brasoe, y eventaliienio ae rormpe on 
varios dorbellinos nas pequetos y cvle nipides, Ml remultade [yal 
(le todo este procosG ca inn transfervencia de energin de los tor 
bellinos pvandes a los torbellinos pequetios, Kolmogorov evaluo 
CUaniLalivamente eate proceso y obtuvoe una expression mateniy 
tlea para la densidad de energia cinética que cabia esperar a cada, 
eseala, es decir, la forma en que se reparte la energia (otal que 
(rinsporta él flujo entre los torbellinos existentes, atendiendo a 
su tamano, Cuando los torbellinos aleanzan un tamaro suficien- 
Lemente pequend, la viscosidad los {rena y los destruye. La ener- 
fit, que es inyectada a gran escala en las grandes estructuras, 
68 disipada en las pequetias escalas. La turbulencia es eficaz a la 
hora de mezclar porque los grandes torbellinos van dando lugar 
a olros de menor tamario cada vez. A las escalas mAs pequefias 
actua eficazmente la difusién, culminando el proceso de mezcla. 

[in las tiltimas décadas se han desarrollado nuevas técnicas de 
manipulacién de liquidos a escala muy reducida, que se conocen 
con el nombre genérico de microflutdica. Este término englo- 
ba los flujos que se producen en microdispositivos, fabricados 
#iracias a las técnicas de la microelectrénica. Las aplicaciones 
tipicas de la microfiuidica son el andalisis quimico de pequefias 
muestras de materiales, la clasificaci6n de particulas 0 macro- 
moléculas y la preparacién de productos farmacéuticos a partir 
de cantidades infimas de principio activo. 

Un problema que aparece con frecuencia en microfluidica es 
el de la mezcla. Por ejemplo, puede ser necesario mezclar la sus- 
tancia que se quiere analizar con un determinado reactivo. Pero 
la escala tipica de un dispositivo microfluidico es del orden de 
una micra, por lo que el nimero de Reynolds es siempre muy 
pequefio. Por tanto, todos los flujos microfiuidicos son laminares 
y la tactica de remover el fluido con una cucharilla microscépica 
no es efectiva, ya que no se produce nunca un flujo turbulento. 

El caos determinista ha venido a resolver este problema. Se 
trata de inducir en el microdispositivo un flujo que, aunque la- 
minax, sea cadtico. El sencillo dispositivo de la figura 1, ideado 
por Hassan Aref, permite demostrar la capacidad de mezcla de 
los sistemas cadticos. El dispositivo consiste en una cavidad 
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Flujo generado por vortices intermitentes (grafico de Emmanuelle Gouillart). 


cilindrica llena de liquido y en la que se encuentran dos barras 
que pueden girar sobre su eje dispuestas como se muestra en la 
figura. Cada barra por separado generaria un flujo en forma de 
vértice cuyas lineas se muestran en la figura la. Para generar un 
flujo ca6tico, las dos barras se mueven alternativamente, primero 
una, luego la otra, y asi sucesivamente. E] resultado es que las 
capas de fluido se estiran y pliegan, dando lugar a un movimiento 
caético. La capacidad de mezcla de este dispositivo se ilustra en 
la figura 1b, donde se muestra cémo se ha dispersado una gota 
de tinta después de varias rotaciones alternativas de las barras. 
En los dispositivos microfluidicos los ingenieros disefan una 
serie de acanaladuras en forma de zigzag que obligan al liquido a 
estirarse y plegarse —la receta del caos— a medida que avanza. 
De esta forma se consigue un mezclado eficaz de las sustancias 
que entran en el dispositivo, aunque el flujo sea siempre laminar. 


EL CADS Y LOS GIRGUITOS ELECTRONICOS 


La teorfa del caos determinista lleg6 a la electrénica de la mano 
del ingeniero electrénico Leon Chua. Nacido en Filipinas en el 
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seno de une farnila de orgen ening, Chum estudiGé an Batidos 
Unidos, donde obtuve un maéster en ingenlerfa electronica en 
él Inatliute Tecnologico dé Massachusetis, Mn 1983 disefé un 
Clreulto Glectrénico que exhibia un comportamiento cadtico. 
El cireuito consistia esencialmente en un oscilador eléetrico 
acoplado a una resistencia no lineal. Un oscilador eléctrico es un 
circuilo en el que la corriente eléctrica circula alternativamente 
én uno y otro sentido a una frecuencia determinada. El circuito 
oscilante més sencillo se compone de un condensador y una in- 
ductancia. En el condensador se acumula carga eléctrica. En la 
inductancia el cable se enrolla sobre si mismo un gran numero 
de veces, formando una bobina, lo que produce en su interior un 
campo magnético al circular la corriente. Al poner en marcha 
un circuito formado por un condensador y una inductancia, la 
corriente circula alternativamente del condensador a la bobina y 
de la bobina al condensador. La frecuencia de la oscilacién vie- 
ne determinada por las caracteristicas fisicas del condensador 
y de la bobina. El otro elemento del circuito de Chua era una 
resistencia no lineal. En una resistencia estandar la diferencia 
de potencial eléctrico entre sus extremos es proporcional a la 
corriente que circula por ella, lo que corresponde a un compor- 
tamiento lineal. En una resistencia no lineal la diferencia de po- 
tencial no es proporcional a la corriente. 

Con la adecuada eleccién de los parametros de los elementos 
que componian el circuito de Chua, la tensién eléctrica entre 
dos puntos iniciaba oscilaciones erraticas de manera esponta- 
nea. Con la ayuda de un osciloscopio, un aparato que permite 
visualizar en una pantalla cOmo varia en el tiempo una sefial 
eléctrica, Chua podia visualizar el atractor correspondiente a la 
variacion cadtica de las tensiones dentro del circuito. La figura 
2 muestra este atractor, obtenido aqui mediante una simulacién 
por ordenador. Su estructura recuerda a la del atractor de Lo- 
renz. Las dos hojas en que se divide son conocidas en el mundo 
de la electrénica por la expresion en inglés double-scroll, que 
podria traducirse como «doble rollo». En los afios ochenta, esta 
figura se convirtié en un simbolo de la electrénica no lineal, en 
la que se hace un uso profuso de elementos circuitales no li- 
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neales. El circuito de Chua abrié todo un campo de la electré- 
nica moderna, que, a partir de entonces, incorporé de manera 
natural las técnicas matemAaticas de los sistemas dinémicos no 
lineales. 


CAOS EN BIOLOGIA 


En el segundo capitulo vimos que uno de los sistemas dinami- 
cos por excelencia estaba directamente relacionado con la bio- 
logia a través de la dinamica de poblaciones: el mapa logistico, 
que Robert May estudi6 en profundidad. Desde ese primer en- 
cuentro entre biologia y caos, los bidlogos se fueron interesan- 
do cada vez mas por el caos y los sistemas no lineales. Ese 
interés ha dado sus frutos en diversas 4reas y ha estimulado el 
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Infente de modelivar matemitloamente ayninos proceioe bie- 
logicos, 

Dos 6rganos de} ser humane han reelbidd especial alencién 
por los expertos en caos: el cerebro y el corazén, in ambos ca» 
sos es posible registrar una serie temporal de datos a partir de 
las sefiales eléctricas generadas por el cuerpo, lo que se con- 
sigue mediante electroencefalogramas y electrocardiogramas. 
La pregunta en ambos casos es la misma: ges posible atribuir el 
origen de estas sefiales a un sistema dinamico con pocos grados 
de libertad? 

El caso de los electroencefalogramas es especialmente com- 
plejo. El cerebro no esta sincronizado como un todo y funciona 
de manera modular, por lo que no se trata solo de la evolucién 
temporal de una sefial eléctrica, sino también de su distribucion 
espacial. El impulso inicial en el estudio del cerebro como un 
sistema dindmico ha llevado a nuevas ideas que se adentran en 
los sistemas complejos y se alejan del tema de este libro. Mas 


, éxito han tenido, desde el punto de vista de la teoria del caos 


determinista, los estudios sobre el comportamiento de una sola 
neurona. Es posible medir la actividad de una neurona registran- 
do el potencial de membrana del axon. El ax6én es la prolonga- 
cidn larga de las neuronas, que conduce el impulso nervioso, y 
el potencial de membrana es el potencial eléctrico que existe 
entre los dos lados, interior y exterior, de la membrana. Cuando 
una neurona es excitada, esa excitacién se propaga como una 
sefial eléctrica que modifica el potencial de membrana a lo largo 
del ax6n. Un filtrado adecuado de la sefial permite determinar 
el momento en que la neurona se excita, y obtener la secuencia 
temporal de los instantes en que se producen las excitaciones, 
De esta forma se han registrado comportamientos caéticos en 
algunas neuronas que ejercen funciones de marcapasos o que 
tienen un papel determinado en ciertas zonas funcionales, como, 
por ejemplo, las neuronas de los llamados generadores centrales 
de patrones, que son redes neuronales que generan conductas 
motoras ritmicas, como la respiraci6n o la locomocién. Un com- 
portamiento cadtico determinista puede ser el régimen normal 
de funcionamiento de estas neuronas. El! estudio de la senal ge- 
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‘ticos tanto temporales como espa- 


O O'all 
ee La doble flecha. expresa. el 
i ue tode i 
de reactivos | y productos que ter gamos & 1A 
leyes dol -equiliprio Uimico. Lo que’ nos interesa aqui es que cualquier reaccion quimica - 
puede contemplarse como un sistema dinamico. En. efecto, si ponemos en contacto dos 


reactivos, A y B, sus concentraciones fan disminuyendo a medida que la reaccion: ‘tiene 
lugar, mientras que las de los productos, C y 0, iran creciendo. Las leyes de ia quimica i 


permiten expresar la evolucion temporal de las: -concentraciones de reactivos y productos nig 


como un sistema de ecuaciones diferenciales. La situacion mas conocicda es aquella en la que- . 
el equiliorio esta muy desplazado hacia la derecha. Porejemplo, en una combustion se mezcla 
el combustible con. oxigeno y se produce agua y didxido de carbono, Esta reaccion tiene lugar: : 
con enorme rapicez, a veces de manera explosiva, por lo-que, vista corno un sistema dinémico, ha 
no tiene mayor interés; pronto se alcanza un estado estacionario en el que los reactivos se han 


~ agotado y la mezcla de productos es inerte. Pero hay situaciones mas complejas desde el a 


punto de vista dindmico. El esis) mas. conocido esla reaocion, de Belousov: Zhabotinsky. 


El descubrimiento de la reaccion BZ ‘ 
En torno a 1950 el quimico ruso Boris Belousay observe, por éssuiatidadl ‘osdlacianes en la E 


reaccion de oxidacién de acido cltrico por bromato de poiasic) catalizada por iones de cerio. _ . 


Belousov noté que [a solucién en ta que se producia la reaccién pasaba de amarilia a incolora 
sucesivamente con regularidad. Una década eecia otro Sie de la Union Sovietica, 
Anatol Zhabotinsky, dedicd su tesis 
doctoral a estudiar en detalle este 
fenémeno. Desde entonces esta 
teaccion ha sido objeto de muchos | 
estudios. Dependiendo ce las con- 
diciones particulares en las que se 
lleve a cabo, se puede observar 
toda una variedad de fenémenos, 
desde las oscilaciones originales 
hasta estructuras espaciales dife- 
renciadas, en las que se alternan 
colores diferentes, regimenes cad- 


ciales, etcetera. 


Patrones espaciales en la reaccion 
Belousov-Zhabotinsky. 
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SI alcanzas ef eguillbrio en blologia =~ nernda puede ayudar a comprender 
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él efecto de drojas o estimulos, Tam: 
blén se han realizado estudios én los 
que una neurona es infegrada en un 
clreuito electronico externo. De esa 
forma ha sido posible generar diversos comportamientos dind- 
micos que no aparecen en la neurona aislada. 

En el caso del corazén, los estudios realizados en el marco 
del caos determinista son mas concliyentes. El ritmo cardiaco 
de un coraz6n sano tiene siempre un componente cadtico, mien- 
tras que los rilmos excesivamente regulares estén asociados a 
patologias. La razén de este hecho parece estar en la robustez de 
los regimenes cadticos. Un ritmo estrictamente periédico no es 
robusto, pequefias perturbaciones pueden desestabilizarlo. Sin 
embargo un régimen dindémico caédtico bien establecido conlleva 
un atractor extrafio, que es un elemento estructuralmente esta- 
ble. El régimen cadtico soporta pequefias perturbaciones y per- 
mite pequefios ajustes sin salir de él. El ritmo cardiaco muestra 
la posibilidad de estabilizar un régimen alrededor de una perio- 
dicidad aproximada por medio de un régimen ca6ético. Este es 
mas resistente a los cambios del entorno. 


ES CAOTICA LA ECONOMIA? 


Durante varias décadas del siglo xx las teorias econémicas han 
estado dominadas por el concepto de equilibrio. Supongamos 
que tenemos dos individuos en disposicién de comerciar y que 
son seres completamente racionales en sus transacciones. Cada 
uno produce un bien diferente, por ejemplo, uno produce vino y 
el otro, aceite de oliva. La produccién anual de ambos es mayor 
de lo que consumen y pueden utilizar el excedente para adquirir 
el producto del otro. Pues bien, segtin la teorfa del equilibrio, el 
que produce vino vendera su excedente al otro a cambio de acei- 
te hasta que lo adquirido iguale en valor a lo vendido. A partir 
de ahi, no hay mas interés en vender y se alcanza un punto de 
equilibrio. Si las condiciones no cambian, los productores inter- 
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Inn esta teoria del equilibnioe es relativamente facil penerali- 
zar Ja sitmacion a muchos agentes. [1 precio de un determinado 
bien queda determinado por el equilibrio entre la oferta y la 
demanda dentro del sistema. Las desviaciones del equilibrio se 
tratan linealmente. Si por alguna razon, por ejemplo, una mala 
cosecha, los precios de un determinado producto se desvian 
del equilibrio, el propio sistema genera fuerzas recuperadoras, 
como lo hace un muelle mecanico, que devuelven la situacién 
al equilibrio. En este esquema, las desviaciones pequenas tie- 
nen consecuencias pequefias. También son posibles situacio- 
nes en las que los precios oscilan periddicamente, en ciclos de 
duracion fija, como los de un péndulo. Las grandes depresio- 
nes, las quiebras bursatiles y otras contingencias se atribuyen 
a efectos externos 0 excepcionales. 

Légicamente, el mundo no es perfecto, ni los agentes econémi- 
cos operan siempre de manera completamente racional. Debido 
a ello, los precios nunca estan en el equilibrio, sino que fluctuan. 
Durante buena parte del siglo xx se pensaba que estas fluctua- 
ciones eran aleatorias y su origen realmente indeterminable. Por 
esta raz6n, los economistas consideraban modelos gaussianos 
de las fluctuaciones econémicas. En un modelo gaussiano las 
fluctuaciones se distribuyen siguiendo la conocida campana de 
Gauss en torno al valor medio (figura 3, izquierda). Cuanto mas 
grande sea una fluctuacién mas improbable es, mientras que las 
fluctuaciones pequefias, por encima o por debajo del valor me- 
dio, son las que ocurren con mayor probabilidad. En la primera 
grafica de la figura 3 se ha considerado que ek valor medio es 
cero, y P(#) es la probabilidad de una fluctuacién de valor x. 
En la grafica de la derecha hemos representado el histograma 
de los precios diarios del barril de petréleo Brent entre enero 
de 2013 y junio de 2014. El histograma se construye anotando 
en el eje de ordenadas el numero de veces que el precio toma 
un valor dentro del rango expresado por las divisiones del eje de 
abscisas. Por ejemplo, durante ese periodo de tiempo el precio 
del barril estuvo entre 110 y 111 ddlares por barril 60 veces, y 
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entre 96 y 97 délares una sola vez. Si los precios del petréleo 
fluctuaran siguiendo una distribucién gaussiana, el histograma 
deberfa parecerse a la curva de la izquierda. Cabria pensar que 
la diferencia se debe a un nimero insuficiente de datos y que, 
tal vez, si tomaramos una serie de varios afios en lugar de los 18 
meses escogidos, el histograma se acercaria cada vez mas a la 
campana de Gauss. Pero no es asi. 

Varios afios antes de acufiar el término fractal, Benoit Mandel- 
brot comprendié que los precios de los productos no fluctuaban 
siguiendo una campana de Gauss. En los afios sesenta del si- 
glo pasado, Mandelbrot analizé la evolucién de los precios del 
algodén en Estados Unidos durante varias décadas y demostré 
que las fluctuaciones no se ajustaban a la campana de Gauss, 
como muchos esperaban. En particular, las fluctuaciones gran- 
des se presentaban con mas frecuencia de lo que predecirfa una 
campana de Gauss. Mandelbrot descubrio que las fluctuaciones 
seguian una ley que presentaba la misma forma a distintas esca- 
las, como sus fractales afios mds tarde. Los precios seguian un 
mismo patrén ya se analizaran por dia, por semana, por mes 0 
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por ahd, Me se correspande con iia ley en be qe las Mietinelo: 
nes grandes, sunque menos probables, slempre son poalbles, Ms 
una ley similar & 1a que alguen los terremotos y otros fenémenos 
catastréficos (los grandes terremotos son mucho menos proba- 
bles que los poco intensos, pero ocurren, luego su probabilidad 
no es nua). 

La irrupeién de la teorfa del caos determinista en diversos 
campos a finales de los afios ochenta del siglo pasado Ilev6 a al- 
gunos economistas, y también a fisicos y matemAticos, a probar 
las nuevas herramientas matematicas en el campo de la econo- 
mia. Tal vez las fluctuaciones de los precios de los productos, 
los tipos de interés o el valor de los indices bursdtiles no eran el 
resultado de fluctuaciones aleatorias de causa desconocida sino 
resultado de una dindmica interna. Si esto era asi, tal vez seria 
posible determinar la dimension del atractor fractal y construir 
modelos sencillos, de solo unas pocas variables, que dieran 
cuenta de toda la dinaémica observada. 

En los casi treinta afios pasados desde la aparicidn del caos 
determinista, los métodos de andlisis han experimentado un 
enorme desarrollo. Cientos de series de datos se han sometido 
a escrutinio, buscando exponentes de Liapunov, dimensiones 
fractales varias, informacion y otras cantidades relacionadas 
con el caos determinista. Pero los resultados son poco conclu- 
yentes. Muchas veces las series de datos no son suficientemente 
largas como para obtener resultados precisos. Otras veces estan 
sujetas a errores diversos, ya sea por la propia medida de la va- 
riable econémica considerada (por ejemplo, el producto interior 
bruto de un pais puede depender de los criterios seguidos para 
calcularlo, y estos pueden diferir entre paises o entre épocas) o 
por falta de fiabilidad de los agentes que las proporcionan. 

Lo cierto es que en ningin caso se puede afirmar con total 
certeza que determinado sistema econdmico sigue las leyes de 
un sistema cadtico determinista. La situacién difiere aqui de la 
que ocurre en otras dreas, como la fisica o la quimica, donde ex- 
perimentos bien controlados en el laboratorio han demostrado 
sin ningtin género de dudas la existencia de caos determinista en 
muy diversos sistemas. 
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Pere a) bien él andlisia de datos reales no ba arrojido con 
clusiones significativas, el cao® determinista he aporinde A ln 
economia ideas y conceplos que Hegaron para quedire, HL mnt 
importante de todos es el de la no linealidad, Sia mediados de} 
siglo xx muchos economistas consideraban que los modelos li 
neales podian explicar los aspectos generales fundamentales de 
la dindmica econémica, hoy en dia existe un consenso general 
acerca de la importancia de las interacciones no lineales entre 
los agentes econdémicos. 


EL PROBLEMA DE LA LIBERTAD 


Aunque histéricamente el problema de la libertad del ser huma- 
no y el libre albedrio se ha considerado una cuestién filos6fi- 
ca, el avance de las ciencias experimentales durante el siglo xx 
ha hecho que se convierta en tema de interés para la ciencia, 
especialmente para la psicologia y, mas recientemente, la neu- 
rofisiologia. Los filésofos de todas las épocas se han apoyado 
en la fisica y las matemdaticas para fundamentar sus posturas al 
respecto. Asi, durante los siglos del determinismo newtoniano, 
los conceptos de espacio y tiempo absolutos y el mecanicismo 
influyeron en los filésofos afines, sobre todo en aquellos que de- 
fendian la predeterminacion. La posibilidad de que el ser huma- 
no sea libre y tenga o no libre albedrio esta relacionada con el 
determinismo. Si la naturaleza es o no determinista y en qué gra- 
do, es una cuestién que no se puede ignorar a la hora de hablar 
sobre la libertad humana. 

El determinismo absoluto de la mecanica clasica, tal y como 
lo expres6 Laplace en su tratado sobre las probabilidades, no es 
compatible con la libertad del ser humano, ni con la capacidad 
del libre albedrio. Segtin decia Laplace: 


Una inteligencia que, en un instante dado, conociera todas 
Jas fuerzas de que la naturaleza esta animada, y la situacion 
respectiva de los entes que la componen, si ademas esa in- 
teligencia fuera lo suficientemente vasta como para analizar 
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locos paoe datos y aber en ine misina (6pm los movi 
nienios de los Herpes nila feandes del univeria asi como 
del tomo mas lifero, ntonees rida seria incierto para ella, 
y tanto el futuro como el pasado se harian presentes ante 
SUS OJOS. 


Laplace contemplaba asi la posibilidad, compatible con la me- 
cénica newtoniana, de calcular con absoluta precisién el movi- 
miento futuro de un cuerpo, conocidas las fuerzas a las que esta 
sometido y su posicién y velocidad iniciales. Hay que decir que, 
aunque esta sentencia ha quedado para la historia como la defi- 
nici6n del determinismo absoluto, el objetivo de Laplace al po- 
ner este ejemplo era precisamente el de ilustrar su imposibilidad 
practica y la necesidad de acudir al calculo de probabilidades en 
muchas facetas de la vida real. 

Ese determinismo absoluto, asi expresado, hace inttil cual- 
quier intento por nuestra parte de escapar a nuestro destino. Si 
las cosas fueran realmente asi, todo estaria escrito y nada de lo 
que nosotros hiciéramos podria cambiar el futuro. Incluso nues- 
tros intentos por cambiarlo no harian sino propiciar lo progra- 
mado, de modo parecido a lo que ocurre con los protagonistas 
de la pelicula Terminator: sus actos no hacen sino provocar lo 
que quieren evitar. 

Pero desde finales del siglo xx la fisica ha ido alejandose 
cada vez mas de esa posicién. Primero con Boltzmann, que in- 
troduce el concepto de probabilidad para explicar las magni- 
tudes termodinamicas, como la temperatura, pero sobre todo 
con la aparicién de la mecénica cuantica, ep la que la probabi- 
lidad figura en el mismo nucleo de sus fandamentos (edricos. 

La mecanica cuantica calcwa probabilidades, Por ejemplo, no 
predice en qué instante va a desintegrarse un nuclea inestable, 
sino que proporciona la probabilidad de que eso ocurra, A la 
postre, y aunque parezca contradictorio, jas predieciones de 
la fisica cudntica se encuentran entre las mas precisas que ha 
producido nunca la fisica, ya que lo normal no es observar un 
evento aislado, como la desintegracion de un atomo, sino una 
gran cantidad de eventos, por lo que normalmente se miden 
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El determinismo no implica que 
cualquier cosa que hagamos no 
podria haber sido hecha de otro 
modo, que todo suceso tenga una 
Causa 0 que nuestra forma de ser 


esté fijada. 


promediog, Hn ewnlquier case, JA meeinicn Midiities no eA una 
teorla determinista, 

Pero no solo In fishen culintica exeluye el determiniamo ab- 
soluto de Laplace. Como aprecié el fil6sofo Karl Popper, la mis- 
ma teoria especial de la relatividad lo excluye. En la teoria de 
la relatividad, la velocidad de la luz es un Iimite maximo a la 
rapidez con que se puede transmitir 
la informacién de un punto a otro. 
Por eso, la inteligencia imaginaria 
de Laplace no podria conocer en un 
instante dado el estado de todos los 
cuerpos que pudieran influir en ella 
en un futuro, por lo que no podria co- 
nocer las condiciones iniciales para 
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adivinarlo. 

Por todas estas razones, la mayo- 
ria de los fisicos de mediados del siglo xx exclufan el determinis- 
mo absoluto como categoria filoséfica. A pesar de ello, la fisica 
cuantica se aplica al mundo microscopico, y la teoria de la rela- 
tividad, a los dos extremos de la escala del mundo conocido: el 
microcosmos y las grandes estructuras del universo. La fisica 
del mundo macroscépico, el de la vida cotidiana, seguia siendo 
determinista, 

Pero, si un mundo absolutamente determinado excluye la po- 
sibilidad de la libertad humana, un mundo completamente alea- 
torio también lo haria. En efecto, si los acontecimientos se su- 
cedieran sin ninguna ley y fueran completamente impredecibles, 
nuestras acciones no servirian de nada. Elegir llevar a cabo o no 
una determinada accion no tendria las consecuencias previstas 
y nos dejaria igual que si no hubiéramos decidido. 

En lo que a la fisica se refiere, el caos determinista viene a 
ofrecer el marco justo que hace posible la libertad: no todo 
esta escrito, pero hay reglas, estructuras dinémicas, que ha- 
cen que el mundo no sea aleatorio. En estas circunstancias 
es posible prever lo que puede ocurrir, calcular o estimar la 
probabilidad de que ocurra, y elegir si se trata de evitar o no. 
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Y ean cleceion londiin efecto en el Matar, dé sederd6 con lis 
propias leyea de In nataralesa, 81 oan libertad existe de verdad 
oO no, si e8 O NO uns ilueidn, Como creen algunos eminentes 
neurofisiGlogos, es una diseusién que excede los propdésitos 
de esta obra. 
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INDIGE 


Val 


La teoria del caos 


Eri las primeras décadas del siglo xx la fisica vivid dos grandes 
revoluciones: la teoria de la relatividad y la mecanica cuantica. 
En las ultimas décadas del siglo tuvo lugar otra revolucién: la 
teoria del caos determinista. Los cientificos de varias disciplinas 
descubrieron que sistemas relativamente simples cuya evolu- 
cidn estaba determinada por leyes precisas podian mostrar un 
comportamiento irregular 0 cadtico. Eritre las caracteristicas de 
los sistemas cadticos estan la sensibilidad a las condiciones ini- 
ciales, e! conocido como «efecto mariposa», y la aparicién de 
estructuras fractales, fascinantes objetos matematicos. 
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